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SPACECRAFT  AND  STELLAR  OCCULTATIONS  BY  TURBULENT  PLANETARY 
ATMOSPHERES 

A  theoretical  investigation  of  various  wave  propagation  effects  and  their  impact  on 
derived  profiles  of  refractivity ,  temperature  and  pressure 


ABSTRACT 

the  long  propagation  paths  involved  in  radio  and  stellar  occultations  by  turbulent 
planetary  atmospheres  require  that  the  classical,  weak  scattering  scintillation  theory 
be  expanded  to  account  for  the  inhomogeneous  ambient  atmosphere  upon  which  the 
turbulence  is  superimposed.  Such  coupling  between  the  turbulent  and  the  ambient 
components  of  refractivity  reduces  the  scintillation  index  by  less  than  a  factor  of  two 
in  shallow  radio  occultations.  For  stellar  occultations  the  reduction  varies  between 
this  value  for  very  small  projected  stellar  radii,  to  approximately  a  factor  of  3.6  when 
the  projected  stellar  radius  above  the  planetary  limb  greatly  exceeds  the  radius  of  the 
free -space  Fresnel  zone  at  this  distance.  More  profound  changes  are  found  in  the 
scintillation  power  spectrum,  the  shape  of  which  depends  strongly  on  both  occulta- 
tion  depth  and  geometry  when  coupling  to  the  inhomogeneous  background  is  pro¬ 
perly  accounted  for. 

Second-order,  systematic  propagation  effects  calculated  from  both  geometrical  optics 
and  a  wave-optical  formulation  show  that  the  average  phase  velocity  is  increased  in 
the  presence  of  turbulence.  The  finite  (X’l/6)  wavelength  dependence  of  the  phase 
path  bias  implies  that  an  initially  non-dispersive  medium  becomes  slightly  dispersive 
by  the  addition  of  turbulence. 

Atmospheric  profiles  derived  from  occultations  by  turbulent  planetary  atmospheres 
differ  only  slightly  from  those  of  the  corresponding  non-turbulent  atmosphere  when 
the  weak  scattering  condition  is  satisfied.  Even  so,  there  is  strong  indication  that 
profiles  obtained  from  radio  Doppler  measurements  are  much  less  affected  by  turbu¬ 
lence  than  those  from  radio  or  stellar  intensity  measurements.  Systematic  errors  in 
refractivity.  temperature  and  pressure  arc  all  small  for  cither  kind  of  measurement  in 
the  limit  of  weak  scattering. 

Radio  and  stellar  o<<  ultations,  when  properly  conducted  and  interpreted,  arc  a  poten¬ 
tial  source  of  information  also  on  turbulence  strength,  bulk  flow  scales,  planetary 
rotation  and  atmospheric  winds. 


1  INTRODUCTION 
1.1  Historical  perspective 

Theoretical  interpretation  of  the  occultation  of  stars  by  planetary  atmospheres 
appears  to  have  originated  with  Pannekoek  (1903),  who  observed  the  occultation  of  a 
star  by  Jupiter  on  September  19  of  that  year.  He  recommended  that  such  observa¬ 
tions  be  made  carefully:  ’’Because  they  can  give  a  determination  of  the  temperature 
distribution  in  the  outermost  layers  of  the  planetary  atmosphere,  or  more  directly, 
the  way  in  which  the  horizontal  refraction  varies  with  height”.  The  theoretical  prob¬ 
lem  of  deriving  a  refractivity  profile  from  a  luminosity  record  of  the  occulted  star 
was  later  reviewed  by  Fabry  (1929),  by  Baum  and  Code  (1953),  and  finally  by  Link 
(1969),  who  was  the  first  to  employ  numerical  inversion  techniques  to  the  optical 
data. 

To  date  observations  of  stellar  occultations  have  yielded  information  on  the  atmo¬ 
sphere  of  Jupiter  through  the  occultation  of  a  Arictis  (Baum  and  Code  1953)  and  the 
(3-Scorpii  system  (eg  Hubbard  ct  at  1972,  Vcverka  ct  at  19741)).  For  a  review  of  stellar 
occultations  by  Jupiter,  sec  Hunten  and  Veverka  (1976).  Structural  information  on 
the  upper  atmosphere  of  Neptune  was  similarly  obtained  by  the  occultation  of 
BI)-17°4388  on  7  April  1968  (Osawa  ct  at  1968,  Kovalevsky  and  Link  1968,  Freeman 
and  I.ynga  1970,  Rages  ct  at  1974,  Vcvcrka  ct  at  1974a).  Stellar  occultations  by  Mars 
(Tcxas-Arizona  Occultation  Group  1977,  Elliot  ct  at  1977b,  Young  1977),  by  Venus 
(dc  Vaucouleurs  and  Men/cl  1960,  Vcverka  and  Wasscrman  1974)  and  by  Uranus 
(Elliot  ct  at  1977a,  Millis  ct  at  1977,  Hubbard  ct  at  1977)  have  also  been  recorded. 


With  the  advent  of  spacecraft  missions  to  the  planets,  occultation  of  the  spacecraft 
itself  has  become  a  very  important  and  complementary  source  of  information  of  the 
structure  of  planetars  atmospheres.  While  stellar  occultations  are  limited  to  fairly  tenu¬ 
ous  atmospheric  regions,  spacecraft  occultations,  by  virtue  of  their  much  higher  signal-to- 
noisc  ratio,  smaller  planet-to-spacecraft  distance,  and  the  Doppler  frequency  as  a  new 
data  source,  have  opened  up  for  study  much  deeper  atmospheric  regions  not  acces¬ 
sible  to  earth-based  optical  observations.  This  subject  was  one  of  several  related  bista¬ 
tic  radar  experiments  introduced  by  V  R  Kshleman  at  the  1962  Review  of  Space 
Research  conducted  at  the  State  University  of  Iowa  by  the  US  National  Academy  of 
Sciences  (Publication  1079),  and  carle  reports  by  Fjcldbo  (1904),  Kshleman  (1964), 
and  Kliore  ct  al  ( 1  964)  describe  details  of  its  scientific  and  engineering  aspects.  Since  the 
first  encounter  of  Mars  by  Mariner  4  on  If)  July  196")  (Kliore  ct  al  196f).  Kjcldbo  d 
al  1966a,  Kjcldbo  and  Kshleman  1968),  the  subsequent  occultations  of  Mariner  a  and 
10  bv  Venus  (<\g  Kliore  ft  ul  1967,  Mariner  Stanford  Group  1967.  Kjcldbo  and 
Kshleman  1969,  Kjcldbo  ct  al  1971,  Howard  ct  al  1974,  Kjcldbo  dal  197a),  of  Pio¬ 
neer  10  and  11  by  Jupiter  (eg  Kliore  dal  1974,  Kliore  dal  197a,  Kliore  and 
Woiceshyn  1976,  Kjcldbo  d  al  i97a),  of  Pioneer  1  1  at  Saturn  (Kliore  d  al  1980)  and 
the  occultations  of  Vovagcr  1  and  2  bv  Jupiter  ( Kshleman  d  al  1979b)  have  yielded 
important  new  results  on  the  structure  of  both  the  atmosphere  and  ionosphere  of 
these  very  planets.  Additional  in  situ  information  on  temperature,  pressure  and  chemi¬ 
cal  composition  have  been  provided  by  the  Venera  landers  descending  into  the  atmo¬ 
sphere  of  Venus  (e  g  Vakhnin  1968,  Vinogradov  d  al  1968)  and  recently  also  by  the 
Pioneer  Venus  entry  probes  ( Seif f  ct  al  1979). 

As  the  potential  viability  of  the  occultation  technique  has  been  demonstrated  in  a 
number  of  both  optical  and  radio  occultations,  a  question  of  utmost  importance 
concerns  the  accuracy  of  the  information  on  atmospheric  structure  obtained  in  this 
way.  As  it  turns  out,  the  answer  to  this  question  is  far  from  simple:  sources  of  widely 
different  character  and  magnitude  contribute  to  the  ultimate  error.  It  is  remarkable 
that  not  until  the  early  analyses  of  the  Pioneer  10  and  11  occultations  by  Jupiter 
produced  grossly  unrealistic  results,  was  the  question  of  accuracy  raised  on  a  serious 
and  systematic  basis  I  eg  Kshleman  197"),  Hubbard  ct  al  197fi). 

Prior  to  these  difficulties,  to  which  1  shall  return  shortly,  stellar  occultation  profiles 
ad  long  been  suspected  of  containing  substantial  errors.  An  outstanding  and  contri¬ 
buting  factor  has  here  been  the  large  sf’ikcs  in  intensity  that  appear  to  be  a  regular 
feature  of  stellar  occultation  light  curses.  Without  doubt,  these  features,  caused  bv 
some  as  yet  unknown  vertical  microstructure  in  the  atmosphere,  have  hampered  the 
derivation  oi  a  refractivity  profile  from  the  occultation  data,  in  particular  in  the 
earlier  days  before  the  numerical  inversion  technique  was  introduced  by  Kovalevsky 
and  Kink  (1968)  to  facilitate  the  analysis  of  the  Bl)-17°4888  occultation  by  Neptune. 
Hut  it  is  also  clear  that  other  major  sources  of  error  are  present.  Perhaps  the  most 
dramatic  example  of  the  difficulties  associated  with  this  technique  is  provided  by  the 
discrepant  refractivity  profiles  obtained  by  several  research  groups  for  the  (3-Scorpii 
occultation  by  Jupiter.  ()1  the  six  recorded  profiles  only  three  seem  to  approach  the 
factual  temperature  distribution  in  the  upper  Jovian  atmosphere  (Vapillon  dal  1978, 
Hubbard  ct  al  19)2.  Veverka  ct  al  1974b).  Fur  a  review  of  the  (5-Scorpii  occultations, 
see  Hunten  anil  Veverka  (1976).  It  is  dear  that  meaningful  results  demand  extreme 
(  are  in  both  the  experimental  design  phase  and  also  during  the  recording  of  data.  It  is 
also  essential  that  due  regard  be  paid  to  the  different  potential  sources  of  error  that 
may  arise  in  the  subsequent  data  reduction  to  obtain  profiles  of  refractivity,  tempera¬ 
ture  and  pressure  (sec  eg  Wasscrman  and  Veverka  1978,  French  dal  1978). 

When  the  first  radio  occultation  profiles  from  the  Mariner  4  Mars  encounter  appeared 
(Kliore  dal  196),  Kjcidho  dal  1966a, b,  Kjcldbo  and  Kshleman  1968),  there  was  no 
reason  to  doubt  their  general  validity,  except  that  the  upper  lew  scale  heights  of 
these  profiles  were  somewhat  indeterminate  because  of  partially  unknown  boundary 
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conditions,  a  problem  common  to  both  radio  and  optical  occultations.  The  sub¬ 
sequent  occultations  of  Mariner  5  and  10  by  Venus  on  19  October  1967  and  5  Tcb- 
ruary  1974,  respectively,  also  yielded  results  that  were  in  all  major  respects  reasonable 
(Kliore  et  al  1967,  Mariner  Stanford  Group  1967,  Howard  et  al  1974).  Indeed,  com¬ 
parisons  between  the  Mariner  5  profiles  and  that  of  the  Venera  4  landing  module 
(Kshleman  ct  al  1968,  Eshleman  1970),  and  between  the  Mariner  5  and  10  profiles 
and  the  Venera  8,  9  and  10  data  (Nicholson  and  Muhleman  1978),  display  a  substan¬ 
tial  degree  of  consistence,  both  between  the  individual  radio  occultation  profiles  and 
between  these  profiles  and  those  of  the  landers  where  a  comparison  could  be  made. 
On  a  finer  level  it  is  nevertheless  clear  that  the  small  differences  that  still  do  exist 
between  the  different  Mariner  10  profiles  are  not  due  to  either  statistical  fluctuations 
in  the  data  processing,  or  in  differing  assumptions  regarding  the  initial  conditions  for 
the  profile  computations  (Lipa  and  Tyler  1979).  Differences  by  as  much  as  10  K 
between  the  profile  computations  of  Nicholson  and  Muhleman  (1978)  and  Lipa  and 
Tyler  (1979)  and  the  initial  profile  of  Howard  et  al  (1974)  probably  require  an 
explanation  in  terms  of  the  difference  in  computational  procedures  used  to  obtain 
these  profiles  (Lipa  and  Tyler  1979). 

The  high  confidence  in  the  radio  occultation  technique,  as  a  tool  for  probing  plane¬ 
tary  ionospheres  and  atmospheres  to  depths  not  accessible  to  optical  occultations,  was 
seriously  shaken  bv  the  publication  of  the  preliminary  temperature  -pressure  profiles 
obtained  from  the  Pioneer  10  and  11  encounters  with  Jupiter  on  3  December  1973 
and  4  December  1974,  respectively  (Kliore  et  al  1974,  Kliore  et  al  1975).  Analysis  of 
the  ionospheric  part  was  hampered  by  multipath  propagation  until  the  several  simul¬ 
taneous  signals  were  laboriously  sorted  out  by  hand  ( Tjeldbo  et  al  1975).  In  the  lower 
neutral  atmosphere,  temperatures  were  far  in  excess  of  those  obtained  from  radiative- 
convective  balance  calculations  and  on-board  radiometric  observations  (Trafton  1973, 
Wallace  ct  al  1974,  Orton  and  Ingersoli  1976).  The  radiometric  observations  fixed  the 
effective  atmospheric  temperature  at  125±2K  at  the  500  mb  level,  far  below  the 
5  400  K  computed  from  the  first  analyses  of  the  Pioneer  1  0  and  1  1  occultation  data. 
This  circumstance,  further  aggrevated  by  Gulkis’  (Bcrge  and  Gulkis  1976)  argument 
that  the  13  cm  S-band  radio  signal  could  not  propagate  below  the  280  K  level  corres¬ 
ponding  to  the  13  cm  atmospheric  brightness  temperature  of  Jupiter,  led  one  to 
conclude  that  the  Pioneer  10  and  11  radio  occultation  profiles  were  in  all  probability 
wrong. 

This  unexpected  situation  spurred  hectic  efforts  to  uncover  the  problem  with  the 
Pioneer  occultation  profiles.  A  major  source  of  error,  the  neglect  of  planetary  oblate- 
ness  in  the  inversion  calculation  to  obtain  a  refractivity  profile  from  the  atmospheric- 
phase  perturbations,  was  identified  by  Hubbard  ct  al  (1975);  sec  also  Hubbard  (1976). 
Almost  simultaneously  Kshleman  (1975)  performed  a  general  parametric  sensitivity 
analysis  of  radio  occultation  experiments,  showing  that  small  errors  in  cither  Doppler 
rate,  spacecraft-limb  distance  or  in  spacecraft  velocity,  were  magnified  by  a  poten¬ 
tially  large  factor  in  such  experiments.  In  the  Pioneer  10  entry  occulation  this  factor 
became  as  large  as  2  70.  Kshleman  also  showed  that  there  was  more  than  one  source 
of  major  error  in  the  original  profiles,  although  correction  for  the  oblatcncss  or 
velocity  factor  did  make  it  possible  to  produce  several  temperature-pressure  profiles 
that  are  both  mutually  consistent  and  in  reasonable  agreement  with  the  Pioneer 
radiometer  observations,  in  particular  in  deeper  regions  of  the  atmosphere,  below 
about  10  mb,  where  the  effect  of  partially  unknown  boundary  conditions  is  rapidly 
decreasing  (Kliore  ct  al  1976,  Kliore  and  Woiccshyn  1976). 

Simultaneously,  and  independently  of  the  events  described  above,  a  different  line  of 
development,  from  which  also  the  present  work  arose,  sought  the  explanation  for  the 
erroneous  Pioneer  profiles  in  terms  of  smaller  scale  refractive  anomalies,  such  as 
turbulence,  in  the  Jovian  atmosphere.  While  layers,  vertically  propagating  internal 
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gravity  waves  or  turbulence  had  long  been  suspected  of  causing  the  large  spikes  seen 
in  stellar  occultation  light  curves,  virtually  no  systematic  effort  was  expended  up  to 
the  year  of  1976  to  explore  these  possibilities.  For  radio  occultation  a  first  quantita¬ 
tive  step  was  taken  by  Hubbard  and  Jokipii  (1975),  who  calculated  the  second-order, 
systematic  effect  of  turbulence  on  radio  occultation  profiles  derived  from  Doppler 
measurements.  Hubbard  and  Jokipii  showed  that  if  exponential  turbulence  was 
imbedded  in  an  otherwise  isothermal  atmosphere,  the  average  radio  phase  would  be 
retarded,  progressively  with  increasing  atmospheric  depth,  and  yield  a  resulting  tem¬ 
perature  profile  that  is  slightly  offset  from  the  true  profile  of  the  quiescent  atmo¬ 
sphere.  The  Hubbard-Jokipii  calculation  was  later  critisized  by  Eshleman  and  Haug- 
stad  (1977,  1978)  and  Haugstad  (1976),  who  showed  that  the  average  phase  bias 
caused  by  the  turbulence  was  in  the  opposite  direction  (i  c,  the  average  phase  was 
advanced),  and  furthermore  that  the  analysis,  contrary  to  their  claim,  was  valid  only 
in  the  limit  of  geometrical  optics  which  is  never  realized  in  radio  occultation  experi¬ 
ments.  It  was  also  demonstrated  that  taking  proper  account  of  diffraction  effects 
both  yielded  a  wavelength  dependence  of  the  phase  path  bias  and  reduced  its  magni¬ 
tude  far  below  the  value  calculated  from  geometrical  optics  (Eshleman  and  Haugstad 
1978,  Haugstad  1978b, c,d,  Haugstad  and  Eshleman  1979).  While  the  precise  form  of  this 
systematic  effect  of  turbulence,  its  magnitude  and  the  question  of  wavelength  depen¬ 
dence,  is  still  under  debate,  both  parties  appear  to  agree  that  this  propagation  effect 
is  too  small  to  be  of  any  practical  concern  (Hubbard  and  Jokipii  1977a,b,  Haugstad 
and  Eshleman  1979,  Hubbard  and  Jokipii  1979,  Hubbard  i979,  Eshleman  and  Haug¬ 
stad  1979). 

Another,  and  in  retrospect  more  important  step  to  clarify  the  effect  of  turbulence, 
was  taken  by  A  T  Young  (1976).  Young  considered  the  larger,  first-order  random 
effects  of  turbulence  and  argued  forcefully  that  the  intensity  spikes  typical  of  stellar 
occultation  luminosity  curves  were  uncorrelated  on  a  planetary  scale  and  were  in  all 
probability  manifestations  of  turbulence  rather  than  layers  or  other  horizontally 
extended  structures.  Developing  a  hybrid  geometrical  optics  —  wave  optics  scintilla¬ 
tion  theory  that  accounted  for  "coupling”  between  the  turbulence  and  the  inhomo¬ 
geneous  background  upon  which  it  is  superimposed.  Young  argued  that  also  radio 
occultation  experiments  should  be  severely  contaminated  by  turbulence  effects. 
Although  Young’s  heuristic  approach  did  not  yield  the  correct  form  of  this  coupling 
effect,  which  appears  as  a  result  of  propagation  over  very  long  distances,  his  paper  did 
bring  into  focus  turbulence  as  a  major  concern  and  undoubtedly  stimulated  other, 
independent  attempts  to  develop  the  correct  theoretical  basis  for  the  first-order, 
random  effects  of  turbulence  (Ishimaru  1977,  Haugstad  1977,  Haugstad  1978b, c,d,e, 
Hubbard  ct  al  1978,  Haugstad  1979a,  Woo  et  al  1980).  These  analyses  showed  that 
Young’s  heuristic  model  seriously  overestimated  the  effect  of  turbulence  and  further¬ 
more  indicated  that  previous  spacecraft  occultation  profiles  may  not  have  been 
seriously  affected  by  turbulence  (Haugstad  1978d,e). 

Along  a  parallel  line  of  development,  interest  was  shifted  from  the  possible 
degrading  effect  of  turbulence  on  radio  occultation  profiles,  to  the  possibility  of 
extracting  from  the  radio  signal  characteristics  of  the  underlying  turbulence.  As 
stressed  by  Woo  et  al  (1974),  information  on  the  strength  and  spatial  extent  of 
turbulence  is  important  as  a  means  of  constraining  theoretical  circulation  models  of 
planetary  atmospheres.  Not  surprisingly,  developing  procedures  for  extracting  such 
information  turned  out  to  be  but  another  aspect  of  the  general  problem  of  electro¬ 
magnetic  wave  propagation  in  a  randomly  inhomogeneous  medium  with  inhomo¬ 
geneous  (and  possibly  anisotropic)  mean  characteristics.  The  pioneering  work  here  was 
done  by  R  Woo  and  A  Ishimaru,  who  applied  standard  weak  scintillation  theory  to 
conditions  typical  of  radio  occultations  of  planetary  spacecraft  (Woo  and  Ishimaru 
1973,  Woo  and  Ishimaru  1974).  The  theory,  not  yet  accounting  for  the  coupling 
effect  discussed  earlier,  was  subsequently  applied  to  the  Mariner  5  and  10  occultations 
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to  yield  rough  estimates  of  the  scintillation  index  and  structure  constant  in  regions 
where  turbulence  could  be  reliablv  identified  (Woo  et  al  1974,  Woo  1975a);  see  also 
Gurvich  ( 1969),  and  Golitsyn  and  Gurvich  ( 1 97 1 ).  The  theory  was  later  expanded  by 
Haugstad  ( 1978c, d,e),  Hubbard  et  al  ( 1978),  and  by  Woo  et  al  (1980)  to  include  effects 
of  coupling  to  the  inhomogeneous  ambient  atmosphere.  Power  spectra  of  phase 
(Haugstad  1979a)  and  intensity  scintillations  (Haugstad  1979a,  Woo  et  al  1980)  were 
calculated  for  a  range  of  different  occultation  geometries  and  discussed  in  terms  of 
their  potential  usefulness  for  deriving  turbulence  characteristics  from  the  radio  occul¬ 
tation  data. 

Even  though  the  past  5—6  years  have  seen  major  advances  in  the  theory  of  turbulence 
effects  in  radio  and  optical  occultation*,  fundamental  problems  still  remain  unsolved.  A 
common  feature  of  all  the  turbulence  theories  quoted  here,  for  example,  is  their 
limitation  to  weak  scattering  conditions.  While  strong  scattering  theories  suitable  for 
short-path,  terrestrial  propagation  conditions  have  been  successfully  formulated  during 
the  past  decade  (for  a  review,  see  eg  Fante  1975),  inclusion  of  the  inhomogeneous 
background  essential  to  radio  and  optical  occultation  conditions,  still  poses  an 
unsolved  problem  in  the  context  of  strong  scattering.  This  and  other  related  prob¬ 
lems,  to  be  considered  in  the  next  section,  undoubtedly  constitute  a  principal  chal¬ 
lenge  to  future  theoretical  work  in  this  field. 


1.2  Objective  and  outline 

The  general  problem  to  be  addressed  in  this  report  is  twofold: 

a)  To  establish,  by  theoretical  analysis,  the  effect  of  turbulence  on  radio  and  stellar 
occultation  temperature  and  pressure  profiles 

b)  To  provide  a  theoretical  basis  from  which  the  characteristics  of  turbulence  in 
planetary  atmospheres  can  lie  denied  from  occultation  measurements 


In  order  to  uncover  the  various  facets  of  this  extensive  problem  and,  at  the  same 
time,  provide  a  justification  for  the  particular  aspects  of  this  problem  that  will  be 
subjected  to  detailed  examination  here,  it  is  convenient  to  consult  the  schematic 
illustration  in  Figure  1.1.  In  order  to  solve  the  problem  specilied  above,  an  essential 


first  step  will  be  to  formulate  a  realistic 


Figure  1.1  Schematic  illustration  of  the  effect 
of  turbulence  on  atmospheric  pro¬ 
files  derived  from  occultation 
measurements 


model  of  both  the  quiescent,  ambient 
atmosphere  and  the  turbulence. 
Once  established,  the  problem  con¬ 
sists  of  calculating  its  effect  on  the 
probing  electromagnetic  signal.  The 
resulting  signal,  contaminated  by 
effects  due  to  turbulence,  is  then 
subjected  to  the  standard  inversion 
algorithm  (i  c,  to  its  mathematical 
representation),  by  which  a  profile 
of  rcfractivity  versus  height  is  con¬ 
structed  from  the  probing  signal. 
Knowing  the  effect  of  the  turbulence 
on  this  signal  allows  one,  at  least  in 
principle,  to  calculate  also  the  corres¬ 
ponding  effect  of  the  turbulence  on 
the  rcfractivity  profile.  If,  further¬ 
more.  the  chemical  composition  of 
the  atmosphere  is  known,  also  pro¬ 
files  of  temperature  and  pressure 
may  be  computed  and  separated,  at 
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least  in  principle  and  in  an  average  sense,  into  an  ambient  part  and  a  fictitious 
component  caused  by  turbulence. 

While  this  procedure  is  in  principle  complete  and  will  solve  the  stated  problem, 
difficulties  arise  on  different  levels  when  practical  solutions  are  sought.  For  instance, 
when  applying  the  result  of  an  analysis  of  this  sort  to  am  specific  occultat ion,  a 
principal  difficulty  arises  from  the  fact  that  the  characteristics  of  the  turbulence, 
which  will  affect  the  outcome  of  the  wave  propagation  analvsis,  is  usually  not  known 
i i  priori,  l  he  problem  is  therefore  iterative  in  nature:  from  a  crude  model  of  the 
turbulence  its  effect  on  the  probing  signal  may  be  computed  and  compared  with  the 
actually  observed  signal.  From  this  comparison  a  refined  turbulence  model,  closer  to 
the  factual  distribution  of  turbulence,  may  in  principle  be  constructed,  etc. 

but  also  major  technical  problems  arise  when  implementation  of  the  scheme  in  Figure 
1.1  is  attempted.  As  in  terrestrial  propagation  experiments,  the  effect  of  turbulence 
or.  electromagnetic  radiation  frequently  falls  into  one  of  the  two  limiting  categories 
of  weak  or  strong  scattering.  In  the  former  case  the  fluctuations  in  signal  strength  arc 
small  in  terms  of  the  mean  signal,  whereas  in  strong  scattering  the  root-mean-square 
lluctuations  are  comparable  to  the  mean  signal  strength  (saturated  scintillations). 
While  weak  scattering  theories  were  fully  developed  bv  the  mid-sixties,  onlv  more 
recently  base  strong  scattering  theories,  due  to  their  greater  mathematical  complexity, 
been  developed  into  explicit  forms  that  allow  confrontation  with  observations.  In 
radio  occultations  two  additional  factors  add  to  the  complexity  of  this  problem: 
ti  the  actual  inhomogeneous  nature  of  the  ambient  atmosphere  cannot  be  neglected, 
a-  is  routinely  done  in  propagation  theories  tailored  to  terrestrial  conditions  involving 
much  shorter  path  lengths,  and  (ii)  the  fluctuations  in  angle  of  arrival  of  the  wave¬ 
front,  a  key  quantity  for  assessing  the  ellect  of  turbulence  on  atmospheric  profiles 
derived  from  Doppler  measurements,  requires  also  the  precise  spatial  extent  of  the 
turbulence  to  be  stipulated  (local  properties  are  not  sufficient)  and  leads  to  a  very 
•  omplicated  analysis.  (For  stellar  occultations  only  the  former  of  these  complications 
applies,  i  We  also  note  that,  because  of  the  non-linear  relationship  between  the  turbu¬ 
lent  refractivity  lluctuations  and  the  amplitude  and  phase,  the  calculated  profiles  of 
temperature  and  pressure  contain  both  a  systematic  and  a  random  part.  During  weak 
-  titering  conditions  the  systematic  effect  will  be  of  second  order  in  small  quantities, 
and  hence  much  smaller  than  the  leading  first-order,  random  effect. 

Having  described  and  categorized  above  the  full  problem,  we  here  restrict  ourselves  to 
giving  a  detailed  treatment  of  certain,  however  important,  aspects  of  the  overall 
problem,  essentially  covering  and  expanding  on  the  material  developed  in  the  quoted 
publications  by  this  author  (Haugstad  1078b, c,d,c,  llaugstad  1070a,  Kshleman  and 
i  laugstad  1977.  F.shieman  and  Haugstad  1978,  Kshleman  and  Haugstad  1979.  Haug¬ 
stad  and  Kshleman  1979,  Haugstad  1980).  I  bus,  giving  an  expose  of  the  basic  theory 
of  radio  and  stellar  occultations  in  section  2.  we  proceed  to  develop  in  section  3  the 
theory  of  electromagnetic  wave  propagation  in  a  non-turbulcnt  atmosphere  as  the 
limiting  form  of  a  general  wave-optical  formulation.  From  this  approach  the  error 
incurred  by  the  routine  neglect  of  wave -optical  effects  in  the  conventional  occultation 
theories  is  readily  assessed.  In  section  4  the  lecak  scattering  theory  of  radiation  propa¬ 
gating  through  a  turbulent  planetary  atmosphere  is  developed  from  a  stationary  phase 
method  (Haugstad  1978c).  After  having  established  an  atmospheric  model,  the 
resulting  formulae  arc  used  to  derive  explicit  results  for  the  weak  scattering  scintilla¬ 
tion  index  in  both  radio  and  stellar  occultation  experiments.  We  establish  here  also 
the  form  of  the  power  spectra  of  the  fluctuations  in  phase  and  intensity,  and  com¬ 
pare  the  predicted  spectra  with  experimental  data  from  the  encounters  of  Mariner  f> 
and  10  with  Venus.  Section  4  concludes  with  a  general  discussion  of  the  validity  of 
the  weak  scattering  approximation,  which  is  a  fundamental  limitation  of  all  the  re¬ 
sults  derived  in  this  work.  In  section  .">  we  investigate  to  what  extent  turbulence,  again 
within  the  realm  of  weak  scattering,  affects  the  derived  profiles  of  temperature  and 
pressure.  Answering  this  question  turns  out  to  be  relatively  simple  for  atmospheric 
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profiles  derived  from  stellar  or  radio  occultation  intensity  measurements.  However, 
for  atmospheric  profiles  derived  from  radio  Doppler  measurements  the  technical  prob¬ 
lems  are  formidable,  as  previously  indicated.  For  this  case  a  physical,  semi-quantita¬ 
tive  approach  avoiding  the  core  of  the  mathematical  difficulties  has  been  adopted. 
The  section  concludes  with  a  discussion  of  the  small  systematic  effects  of  turbulence 
on  occultation  profiles.  In  section  6  wc  devote  attention  to  a  relatively  new  topic  in 
this  field:  the  effect  of  turbulence  in  deep  radio  occultation  measurements.  The  basic- 
theory  of  deep  radio  occultations  (Eshleman  et  al  1 979a)  is  combined  here  with  the 
scintillation  theory  developed  in  section  4  to  yield  results  for  the  scintillation  index 
anil  p>>u«i  spectrum  close  to  a  refractive  focus  (Haugstad  1981).  The  theory  is 
applied  to  deep  radio  occultations  o|  planetary  atmospheres  (the  Voyager  1  Jupiter 
encounter),  and  to  t he  suggested  use  of  the  gravitational  lens  of  the  sun  for  eaves¬ 
dropping  and  communication  over  interstellar  distances  (Eshleman  1979).  We  con¬ 
clude  the  present  work  in  section  7  by  discussing  a  general  procedure  by  which 
information  on  the  strength  and  spatial  distribution  of  turbulence  may  be  extracted 
from  occultation  measurements,  exploring  both  the  potential  feasibility  and  the  inhe¬ 
rent  limitations  of  this  scheme.  Section  8  summarizes  the  results  of  the  present  work, 
and  expands  the  perspective  by  recommending  where  future  work  in  this  field  should 
be  directed. 


2  RADIO  AND  STELLAR  OCCULTATIONS;  FORMAL  THEORY 
2.1  Profiles  of  refractivity,  temperature  and  pressure 

1  he  derivation  of  vertical  profiles  of  temperature  and  pressure  from  occultation 
measurements  proceeds  in  two  stages.  In  the  first  a  refractivity  profile  is  determined 
from  the  atmospheric  perturbations  of  either  phase  or  intensity  of  the  probing  signal. 
In  the  second  stage  the  refractivity  profile  is  used  together  with  possible  knowledge 
about  the  chemical  composition  of  the  atmosphere  to  calculate  profiles  of  temperature 
and  pressure. 

For  the  first  of  these  steps,  consider  Figure  2.1  which  illustrates  a  spacecraft  partially 
occulted  by  a  planetars'  atmosphere,  the  refractivity  of  which  is  assumed  to  be  a 


Figure  2  I  Ray  Rath  and  simplified  occultation  geometry'  defining  the  rax  h ending 
angle  a.  the  impact  parameter  a,  the  radial  distance  to  the  point  of  closest 
approach  (periapsis)  r0,  and  the  planet  to  spacecraft  distance  I)  as  measured 
from  the  planetary  centre 

The  spacecraft  velocity  is  resolved  into  a  component  parallel  to  the  incoming  rav 
asymptote,  and  a  component  vt  in  the  plane  of  the  sky. 

I 

j 
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function  of  the  radial  coordinate  only-  (For  stellar  occultations  the  earth  replaces  the 
spacecraft  and  1)  is  now  the  planet-to-earth  distance.)  It  is  assumed  throughout  that 
ro  >  r0c,  where  roc  is  the  level  of  superrefraction,  below  which  the  ray  encircles  the 
planet  and  fails  to  reach  the  spacecraft.  (For  a  discussion  of  this  interesting  possibi¬ 
lity,  see  Lipa  and  Croft  (1975),  who  also  calculate  the  refractive  attenuation  of  the 
critical  ray.) 

Let  the  refractive  index  n  =  1+-  v,  where  f  js  the  atmospheric  refractivity.  The  quanti¬ 
ties  n  and  a  are  related  by  an  Abelian  transform  pair  (eg.  Fjeldbo  et  al  1971, 
Phinncy  and  Anderson  1968) 


=  —2a  7  f  -tjAtt, 

a  d?  n(£  -  a^ ) 

(2.1a) 

7r  In  n  (r0)  =  —  /  ^7^  cosh'1  {-)  da' 
a  da  a 

(2.1b) 

where  £  =  nr,  and  where  r0  and  a  are  related  by  Bouguer's  rule  (Born  and  Wolf  1975) 
as  in 

a  =  n(r0)r0  (2.2) 

Let  k  (k  =  1  k  I  =  2tr/X)  be  the  vector  wavenumber  of  the  radio  wave  after  refraction. 
If  the  planet  is  assumed  stationary  with  respect  to  the  earth,  the  additional  Doppler 
frequency  f  caused  by  atmospheric  refraction  through  an  angle  cv  is 

Af  -  k  ■  v  -  v^lcosa  ll  *  v,cosd  sino  (2. 5a) 

where  k  r  k/k.  and  where  0  is  the  (positive)  angle  between  v,  and  the  plane  of 
refraction.  If  v  is  lying  entirely  in  the  plane  of  propagation  ( <7  =0 ).  the  simpler  first 
order  form  of  equation  (2.5a)  is 


Af  =  VjQ  (2.5b) 

From  extrapolation  of  radio  tracking  data  during  times  when  there  is  no  atmospheric 
perturbations  of  the  radio  wave,  the  position  and  velocity  of  the  spacecraft  are 
known.  According  to  equations  (2.5a)  or  (2.3b;.  a  measurement  of  the  atmospheric 
Doppler  f  yields  the  bending  angle  a  and  hence  the  impact  parameter  a.  From  corres¬ 
ponding  values  of  a  and  a,  the  differential  da/da  may  be  formed  and  hence,  upon 
performing  the  indicated  integration  in  equation  (2.1b).  the  refractive  index  correspon¬ 
ding  to  a  g  cen  impact  parameter  is  obtained.  Application  of  the  Bouguer-relation 
(2.2)  finally  yields  the  desired  n(r0 (  profile. 

Calculation  of  a  refractivity  profile  from  stellar  or  radio  occultation  intensity 
measurements  follows  an  essentially  similar  path,  lhe  fundamental  equation  is  still 
(2.1b),  but  now  da/da  is  expressed  in  terms  of  the  intensity  of  the  probing  signal. 
Indeed,  for  small  bending  angles  and  neglecting  focusing  of  the  wavefront  around  the 
curved  limb  {al)  <  R),  geometrical  optics  flux  conservation  yields  for  the  intensity  0 
of  the  near-limb  ray 

1  0'1  =  Dda/da  (2.4) 

when  0  is  normalized  to  unity  before  refraction.  Thus,  since  a(t)  is  known,  a  measure¬ 
ment  of  0  versus  time  t  establishes  da/da,  and  via  equations  (2.1b)  and  (2.2)  also  the 
refractivity  as  a  function  of  atmospheric  depth. 
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For  an  isothermal  atmosphere,  where  v  varies  exponentially  with  height  with  a  con¬ 
stant  scale  height  11,  equation  (2.4)  may  he  manipulated  to  yield 

^  =  (l-2)+ln(i-1)  (2.5) 

110  0 

Here  is  the  component  of  the  spacecraft  velocity  in  the  plane  of  refraction  that  is 
perpendicular  to  the  ray  before  refraction.  The  time  t  is  taken  to  be  zero  when 
0=  1/2,  i  r,  corresponding  to  the  nominal  occultation  level.  An  essentially  identical 
form  of  equation  (2.5)  was  first  derived  by  Baum  and  Code  (1953)  and  has  since 
been  used  extensively  in  the  interpretation  of  stellar  occultation  measurements.  Thus, 
lor  an  isothermal  atmosphere  the  intensity-time  profile  0(t)  yields  the  temperature  T 
directly  if  the  mean  molecular  mass  m  is  known,  since  T  =  mgH/K,  where  g  is  the 
acceleration  of  gravity  and  K.  is  Bolt/.man’s  constant. 

Equations  (2.4)  and  (2.5)  strictly  apply  only  to  a  point  source.  For  stellar  occulta- 
tions  the  intensity  0  should  be  convolved  with  the  brightness  distribution  of  the 
projected  stellar  image  in  the  occulting  atmosphere.  As  shown  by  P  J  Young  (1977), 
however,  the  resulting  corrections  for  the  finite  angular  width  of  the  star  are  normally 
altogether  negligible 

Once  a  refractivity  profile  has  been  obtained  from  either  Doppler  or  intensity 
measurements,  the  hydrostatic  equation  may  be  integrated  to  yield  the  temperature 
profile  T(r0),  and  by  the  perfect  gas  law  also  the  pressure  profile  p(r0).  In  a  well- 
mixed  atmosphere  of  non-polar  molecules  we  have  (eg,  Fjcldbo  and  Eshleman  1968) 

l'(r0 )  =  T(r0 , )  +  -  n\-  /  g(r’0 )  ^(r'0  )  drj,  (2.6) 

v(r0)  k^(r0)r0l 


p(r0  )  ~  ^  -r-->  f(r0 )  (2.7) 

tn 

In  these  equations  m  is  the  (constant)  mean  molecular  mass  of  the  atmospheric 
constituents,  v  is  the  mean  "refractivitv  volume”  (i  e.  the  refractivity  at  standard 
conditions  divided  by  Eoschmidt’s  number),  and  K  is  Boltzman’s  constant.  Equation 
(2.6)  assumes  a  boundary  condition  high  in  the  atmosphere  at  r0  =  r0,. 

Having  reviewed  the  formal  occultation  theory  it  is  appropriate  to  add  a  few  com¬ 
ments  relative  to  our  objective,  as  stated  in  section  1.2.  We  note  first  that  the 
procedure  described  above  assumes  the  atmospheric  refractivity  to  be  spherically 
stratified.  Turbulence,  representing  small  inhomogencities  in  v  over  a  continuum  of 
scale  sizes,  clearly  violates  this  provision.  One  might  think  that,  since  occultation 
experiments  arc  sensitive  to  atmospheric  properties  over  a  distance  C  =  (RH)12  in  the 
direction  of  the  ray,  the  calculated  refractivitv  profile  when  turbulence  is  present 
would  contain  an  additional  component  equal  to  the  average  turbulent  refractivity  over 
the  distance  C.  This  is  not  the  case,  however.  As  indicated  in  section  1.2,  turbulent 
refractivity  fluctuations  in  fact  produce  changes  in  signal  strength  that  grow'  with 
distance  and  may  eventually  become  comparable  to  the  mean  occulted  intensity,  even 
though  the  turbulent  refractivity  fluctuations  arc  a  small  fraction  only  of  the  ambient 
refractivity.  The  intensity  fluctuations  also  depend  on  the  radiation  wavelength  X, 
whereas  the  geometrical  optics  intensity  formulae  (2.4)  and  (2.5)  by  their  nature  arc 
independent  of  X.  We  thus  conclude  that  turbulent  refractivity  fluctuations  cause  the 
derived  refractivity  profile  to  depart  from  the  "true"  profile  of  the  ambient  atmo¬ 
sphere  in  ways  that  may  be  largely  independent  of  the  generic  fluctuations  in  refractivity 
at  the  ray  periapsis.  For  these  reasons,  exploring  the  elici  ts  of  turbulence  on  atmo¬ 
spheric  profiles  obtained  from  occultation  measurements  constitutes  both  a  legitimate 
and  important  scientific  endeavour. 
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2.2  Effects  of  boundary  conditions  and  various  measurement  errors 


For  practical  reasons  the  integration  in  equation  (2.1b)  has  to  be  truncated  at  some 
upper  limit  a’=ai,  say.  For  an  isothermal  atmosphere  the  resulting  fractional  error  in 
derived  refractivity,  Su,  is  approximately  (Wasserman  and  Vevcrka  1973) 

5„  =  1  -  erfl0(rol  -  r0)|  (2.8) 

where  |3'1  =  11  is  the  scale  height,  erf(-)  is  the  error  function,  and  r0!  is  related  to  U| 
via  Bouguer’s  rule.  The  form  of  equation  (2.8)  implies  that  the  refractivity  error 
decreases  rapidly  with  increasing  penetration  of  the  rav;  for  r0)  r0  =  211,  &„  -  3%. 
while  three  scale  heights  below  the  starting  point  of  the  integration  &v  s  1 .5%.  F'rom 
equations  (2.6)  and  (2.7)  it  is  found  that  also  the  associated  temperature  and  pressure 
errors  decrease  rapidly  below  the  first  few  scale  heights  of  these  profiles. 


An  additional  error  in  F  or  p  arises  from  the  essentially  unknown  boundary  condition 
in  temperature.  From  equations  (2.6)  and  (2.7)  it  is  seen  that  an  error  AT(r0i  )  in  the 
boundary  temperature  yields  errors 

AT(r0)  =  AT(r01)  ~~  (2.9) 

Kro ) 

and 

ApM  AT(r0l)  (2.10) 

m 

in  temperature  and  pressure,  respectively,  at  depth  r0.  While  the  temperature  error 
decreases  rapidly  with  depth  (eg,  exponentially  for  an  isothermal  atmosphere),  the 
absolute  pressure  error  is  constant  with  depth.  The  fractional  error  in  pressure 
decreases  rapidly  with  atmospheric  depth,  however. 


The  effect  of  different  types  of  measurement  errors  on  the  determination  of  p  and  T 
has  been  considered  by  Eshleman  ( 1 973),  and  in  part  also  by  Hubbard  ct  al  (1975). 
Following  Flshlcntan  (1973),  the  temperature  scale  height  of  an  assumed  isothermal 
atmosphere  may  be  expressed  in  terms  of  spacecraft  trajectory  data,  the  radio 
I  hippier  f,  and  its  time  derivative  df/dt  as  in 


MH(  v^2 - l)  -  l 

V1H  XD  |— | 

dt 

Alternatively  this  relation  may  be  written  as 


(2.11) 


M(S  -  1)  =  1  (2.12) 

where  the  dimensionless  magnitude  factor  M  and  the  sensitivity  factor  S  arc  defined 
by  comparison.  Since  M  for  planetary  occultations  may  reach  a  hundred  or  more.  S 
only  slightly  exceeds  unity  over  most  of  the  atmospheric  region  probed.  As  pointed 
out  by  Flshleman  (1975),  S  is  identically  unity  for  bending  mound  a  fixed  limb 
(11=0),  as  in  classical  knife-edge  diffraction.  All  information  about  atmospheric  struc¬ 
ture  is  therefore  contained  in  the  small  excess  of  S  over  unity. 

Referring  to  equations  (2.11)  and  (2.12).  absolute  errors  AS  in  S  and  A(aD)  = 
A(Xf/vj )  in  the  deflection  of  the  rav  produce  a  fractional  error  in  11  given  b\ 

H  OD  1 

Thus  fractional  errors  in  any  of  the  quantities  vf ,  I),  X  or  Idf/dtl  comprising  S  arc 
magnified  by  the  potentially  large  factor  M  =  ftl)/H,  while  a  fractional  error  in  a  (re, 
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in  f)  produces  only  a  comparable  error  in  H.  If  |df/dt|,  D  or  X  equals  (1— A)  times  its 
true  value,  or  if  the  derived  v j2  is  (1+A)  times  its  true  value,  then  to  first  order  in  A 
and  for  large  M  (Eshleman  1975) 


l>f 


=  1  +  MA 


T  3/2 


(2.14a) 


(2.14b) 


where  I  f  and  I't  denote  the  false  and  true  temperature  profiles,  respectively,  and  with 
a  similar  notation  for  p. 


Despite  the  potentially  large  error  in  derived  profiles  induced  by  small  errors  in  S, 
temperature/pressure  profiles  derived  from  Doppler  measurements  are  still  more  accu¬ 
rate  than  corresponding  profiles  derived  from  radio  or  stellar  occultation  intensity 
measurements,  provided  that  the  radio  Doppler  can  be  measured  at  least  M  times 
more  accurately  than  the  intensity,  which  is  normally  the  case.  From  equation  (2.5) 
it  is  indeed  found  that  errors  in  either  v^  or  0  yield  comparable  errors  in  the  scale 
height  determination,  and  hence  in  T.  This  is  also  the  reason  why  stellar  occultations 
are  routinely  terminated  when  the  residual  light  flux  is  comparable  to  the  expected 
error  in  baseline  determination,  which  is  typically  a  few  per  cent. 

Digitization  and  subsequent  data  processing  also  introduce  errors  in  the  final  T— p 
profiles.  Such  errors  arc  highly  dependent  on  signal  to  noise  ratio  during  data  trans¬ 
mission,  and  on  the  specifics  of  the  data  processing  algorithms  used,  as  discussed  in 
detail  by  l.ipa  and  I  yler  (1979). 

The  previous  analysis  apparently  allows  the  following  arguments  to  be  made  concern¬ 
ing  the  effect  of  turbulence  on  occultation  profiles.  Turbulence  causes  the  actually 
observed  bending  angle  a  to  deviate  slightly  from  the  "true”  value  of  the  quiescent, 
non-turbulent  atmosphere.  According  to  equation  (2.13).  a  change  in  a  produces  a 
comparable  (t  c.  small)  error  in  H  and  hence  in  T  if  these  quantities  are  obtained 
from  radio  Doppler  measurements.  Atmospheric  profiles  derived  from  intensity 
measurements,  on  the  other  hand,  contain  errors  that  are  comparable  to  those  in 
intensity  according  to  equation  (2.5).  Since  intensity  scintillations  due  to  turbulence 
may  become  comparable  to  0  itself,  it  follows  that  the  resulting  temperature  error 
may  also  be  appreciable,  unless  some  "smoothing”  procedure  is  applied  that  suppresses 
the  scintillations.  We  shall  see  from  the  more  stringent  analysis  in  section  5  that  these 
contentions  are  indeed  true. 


3  ELECTROMAGNETIC  WAVE  PROPAGATION  IN  A  NON-TURBULENT  PLANE¬ 
TARY  ATMOSPHERE 

3.1  General  formulation 

Standard  formulations  of  electromagnetic  wave  propagation  in  a  spherically  sym¬ 
metric,  non-turbulent  atmosphere  arc  based  on  geometrical  optics.  The  rationale  for 
this  is  the  tact  that  the  radius  of  the  first  Fresnel  zone  is  smaller  than  the  atmo¬ 
spheric  scale  height  by  typically  a  factor  on  the  order  of  10,  cf  Table  3.1.  We  depart 
here  from  the  common  geometrical  optics  approach  by  formulating  the  propagation 
problem  in  a  wave-optical  context,  from  which  the  proper  approximations  leading  to 
the  standard  geometrical  optics  results  will  appear  in  a  natural  way.  The  small  wave- 
optical  corrections  to  these  results  arc  also  readily  obtained  from  this  approach, 
which  also  provides  the  notational  framework  for  the  generalization  in  section  4  of 
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30 

6.5-10* 

0.57 

BD17°43 88/NEPTUNE 

60 

4.5  10y 

1,6 

Table  3.1  Representative  values  for  various  radio  and  stellar  occultations 

The  table  Rives  representative  values  for  the  scale  height  H,  approximate  plancl-to- 
spacecraft  (or  earth)  distance  I),  radius  of  the  free-spat  c  Fresnel  zone  up  n .  and 
projected  linear  radius  as  n  of  two  occulted  stars.  All  lengths  are  in  kilometers.  The 
two  stars  are  observed  at  A  =  O.'j/im.  while  X  =  0.1  3  m  (S -band)  has  been  assumed  in 
the  calculation  of  a J.  „  for  the  radio  occultations. 


PHASE  SCREEN 
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Figured.!  .1  phase  changing  screen  oj 
thickness  represents  the 

inhomogeneous  ambient  atmo¬ 
sphere 

Position  vectors  I,  anti  T2  refer  to  the  exit  fate  of 
the  screen  and  an  arbitrary  field  point  a  distance 
, r 2 1  2  1)  C  from  the  sireen.  respei  lively 


If  KIT, 


cxp|ik  |r: 
ir2  r,  | 


the  weak  scattering  scintillation  theory 
to  include  the  inhomogeneous  ambient 
atmosphere. 

We  shall  assume  throughout  that  the 
effect  of  the  actual  atmosphere  on 
electromagnetic  waves  may  be  repre¬ 
sented  by  a  thin  slab  of  thickness  C  <  I) 
tnd  refract ivity  (V  the  effect  of  which 
is  to  change  the  phase  of  the  wave  but 
leaving  its  amplitude  unchanged,  see 
Figure  3.1.  We  shall  find  that  this  thin 
screen  approximation  is  always  valid 
for  ray  bending  angles  a  <S  (H/R),/2 , 
where  II  is  a  typical  scale  of  variation 
of  the  spherically  stratified  atmosphere. 

Let  F.(r, )  be  the  complex  electric  field 
of  the  wave  in  the  plane  x  =  C/2  upon 
exit  from  the  phase  screen.  In  position 
r2  behind  the  screen  the  field  is  (eg, 
Goodman  1968) 


cos(n.c)  dr, 


where  dru  -  dy,d/.,  is  a  differential  area  on  the  exit  plane  x  =  C/2  of  the  screen.  The 
obliquity  factor,  cos(n,e),  accounts  for  the  angle  between  the  normal  n  to  the  phase 
screen  and  the  direction  e  =  (r2  r,  )/lr2  r,  I  to  the  observer.  If  the  spatial  extent  of 
the  refracting  medium  perpendicular  to  the  ravpath  is  small  compared  to  D.  the 
angle  Ln,c'  0  and  cosfn.e)  '  1,  while  lr2  r,  I  2  1).  This  paraxial  approximation  is 
assumed  throughout.  Consistent  with  the  thin  screen  approximation,  an  incident  plane 
wave  of  unit  amplitude  has  an  electric  field  at  exit  given  by  F.(r,  )  =exp(ikS0),  where 
S0  is  the  change  in  phase  path  through  the  screen.  We  shall  see  in  section  3.2  that, 
consistent  with  the  thin  screen  approximation,  the  phase  change  may  be  calculated  by 
integrating  un  along  a  straight  line  path  through  the  screen,  instead  of  following  the 
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actual  curved  raypath.  Thus  S0  =  /^0(x,y0,z0 )dx,  where  (y0,z0)  arc  the  lateral 
coordinates  of  the  ray  upon  entering  the  screen.  If  C  =  (27fRH)|/2  ,  it  also  follows  that 
S0  =  i>0C,  where  v0  is  now  referred  to  the  ray  periapsis  of  the  actual  exponential  atmo¬ 
sphere  of  scale  height  11. 

With  these  provisions  equation  (3.1 )  may  be  expressed  in  t he  form 
i  •>  •  —  i j/ 

E(r2)  =  -*-//e  20  dru  (3.2a) 

2ttiD  -« 


where 


4>  =  2D(L)+  S0)  +  (y,  y,)2  +  (z2-z,)2 


(3.2b) 


Following  Haugstad  (1978c)  (sec  also.  Hubbard  and  Jokipii  1 977b),  we  expand  \p 
around  the  point  of  stationary  phase  (y0 ,z0 )  of  the  integrand,  i  c 


4*  =  -  Cy(y i  y0)'(/.|  z0)’ 

i.j  =  o  ‘J 


(3.3a) 


where 

c  =  _L  v  _ 

1J  l!i!  d\Jd*i  . '>0 

*1  -'-Cl 

The  stationary  phase  [joint  is  defined  by  the  conditions  c0 1  (y  i  =>'o . 
cio(Vi=y0,  z,  =z0)  =  0,  yielding 

Vo  =  V:  -  l>0  .  z0  =  z- 2  -  Ho 


(3.3b) 
Z|=70)  = 


where  a  and  (3  are  the  ray  bending  angles  in  the  xz-  and  xy-planc,  respectively,  as 
given  by  ot  =  jV0  /  (x.y0.y0)dx  and  (3  =  ft'o.y  (x,y0,z0)dx.  (Subscripts  v,  z,  yy,  yz,  etc 
denote  partial  differentiation  in  the  usual  way.)  Successive  determination  of  the  re¬ 
maining  coefficients  Cy  yields,  through  second  order 

Coo  =  2I)2[1  +  i(a2  +  (32)]+  2DS0 

Cot  =  c,0  =  0  (3.4) 


c02  1  +  M  .  C] ,  -  0  ,  c20  -  1  +  N 


I  he  two  quantities  M  and  N  play  a  central  role  in  the  subsequent  analysis;  they  are 
defined  by 


\1  - /^o.zz(x’>'o >70 )dx  ,  N  -  /r'o,yy(x,y0,z0)dx  (3.5) 

where  the  integration  is  again  understood  to  be  through  the  screen  from  x  =  —C/2  to 
x  =  C/2. 

Equations  (3.2)  (3.5)  constitute  the  basic  mathematical  framework  for  the  sub¬ 
sequent  analysis  in  this  and  also  in  the  next  main  section. 


3.2  Geometrical  optics 

We  shall  assume,  tentatively,  that  the  stationary  phase  expansion  (3.3)  can  be  trun¬ 
cated  after  the  second-order  terms,  i  e,  we  define  Cy  ( i  +  j  >  2)  =  0.  In  this  case  equa¬ 
tion  (3.2a)  becomes 


I 
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E(r2) 


ke 


>k^-° 

2D 


2  ff  i  D  -~ 


//e 


^  I0y  (yi -yo)2  +  <t> i1  (zi  zo )2 1 


dr 


lx 


(3.6) 


where  we  have  defined 

c20  -  0y  =  1  +  N  ,  c02  =  *-'  =  1  +  M  (3.7) 

The  integrals  in  equation  (3.6)  are  of  the  Fresnel  type  when  evaluated  between  the 
limits  ±  00  as  indicated.  Apart  from  possible  contributions  from  the  neglected  terms  in 
the  stationary  phase  expansion,  the  very  presence  of  the  planet  itself  requires  that  a 
circular  region  of  radius  R,  centered  at  Z|  =  r0,  be  excluded  from  the  integration. 
However,  due  to  the  rapidly  decreasing  contribution  to  the  integral  from  wavelets 
originating  from  positions  |z(  —  z0l  >  (XD)1/2  0Z,/2  and  Ivi  y0l  >  (XD/2  0'^  from 
the  stationary  phase  point,  the  presence  of  the  planet  can  be  ignored.  Evaluating 
equation  (3.6)  for  the  intensity  0  =  IE|2  thus  yields 

0  =  0y0z  (3.8) 

This  result,  independent  of  k,  is  the  geometrical  optics  intensity  formula  for  a  spheri¬ 
cally  symmetric  refractivity  field.  It  can  be  derived  from  a  flux  conservation  argument 
directly,  or  from  a  geometrical  optics  formulation  based  on  the  Eikonal  equation 
(Haugstad  1978e).  Truncation  of  the  stationary  phase  expansion  after  the  second- 
order  terms  therefore  restricts  the  validity  of  the  resulting  intensity  formula  to  the 
realm  of  geometrical  optics.  We  calculate  in  the  next  section  the  lowest-order, 
k-dependent  correction  to  equation  (3.8),  and  thereby  also  establish  the  region  of 
validity  of  this  result. 

Equation  (3.6)  provides  an  interesting  wave-optical  interpretation  of  ray  bending  in 
an  inhomogeneous  medium.  Indeed,  the  phase-factor 

0  =  ^l0y  (y  t  ->’o)2  +  <t>z  (Z|  -Zo)2 ) 


has  the  constant  value  tr  on  the  circumference  of  an  ellipse  centered  at  (y0.Zo)  with 
semi-axes 


aF-h  =  0^(XD)'/2  ,  aFtV5  0«'2(XD),/2  (3.9) 

in  the  y-  and  z-directions,  respectively.  Refraction  thus  distorts  the  originally  circular 
shape  of  the  frce-space  Fresnel  zone  of  radius  aF  ()  =  ( Al>) 1/2  into  an  ellipse  with 
semi-axes  given  by  equation  (3.9).  The  property  that  the  intensity  is  proportional  to 
the  square  of  the  area  of  the  first  Fresnel  zone  is  preserved,  however,  since 
(n’aF,ha f.v/^F.O  )2  =  0y0z>  *n  agreement  with  equation  (3.8).  We  note  also  that 
diffraction  effects,  through  the  inclusion  of  higher-order  terms  in  the  stationary  phase 
expansion  and  in  0  above,  will  cause  the  atmospheric  Fresnel  zone  to  deviate  slightly 
from  its  strictly  elliptical  shape  in  the  geometrical  optics  limit. 

For  a  spherically  symmetric  and  exponential  atmosphere  the  definitions  (3.5)  may  be 
used  to  show  that 

M  =  aD/H  ,  N  =  -aD/R  (3.10) 

in  this  case.  The  corresponding  result  for  the  intensity  is 

0’ 1  =  (1  +  5(D)  (1  -5®) 

'  H  ' '  R 


(3.11) 
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which  reduces  to  equation  (2.4)  when 
focusing  is  negligible  (al)  R).  In  the 
present  case  refraction  expands  the 
frce-space  Fresnel  zone  by  a  factor 
( 1  -oD/R)'1'2  >  1  parallel  to  the  limb, 
while  at  the  same  time  compressing 
the  vertical  scale  by  the  factor 
(1  +  aD/H)"1'2  <  1,  see  Figure  3.2.  As 
aD  -*■  R,  the  spacecraft  approaches  the 
center  of  the  planet,  as  viewed  in  the 
plane  of  the  sky,  and  equation  (3. 1 1 ) 
predicts  Diffraction  effects 

limit,  however,  the  intensity  to  a  finite 
value  whose  magnitude  is  readily  ob¬ 
tained  from  the  physical  picture  de¬ 
veloped  above  (cf,  Eshleman  ct  al 
1979a). 

Fhc  phase  in  the  geometrical  optics 
approximation  is  obtained  from  equa¬ 
tion  (3.6)  by  writing  the  electric  field 
as  F.  =  A  exp  |i(kS  -o>t)|,  where  A  is 
the  amplitude  and  kS  is  the  phase.  Per¬ 
forming  the  integration  over  the  screen 
it  is  found  that 

kS  ~  k|l)  ‘  S0(y„ .•/.„)  t  4l)(a2  t  tf2)|  (3.12) 

1  he  sum  kl)|  1  *  l/2(or  t)l2)|  is  simple  the  increase  in  phase  by  following  the  geomet¬ 
rical,  refracted  path  to  the  receiver,  while  kS0  is  the  additional  atmospheric  phase 
perturbation  inside  the  refractive  medium. 

The  previous  developments  are  all  based  on  the  thin  screen  approximation,  as  defined 
in  section  3.1.  Equation  (3.8),  or  more  conveniently  equation  (3. 1  1 ),  may  now  be 
used  to  explore  the  validity  of  this  approximation.  For  this  approximation  to  be  valid 
it  is  sufficient  that  the  change  in  intensity  <!>  through  the  medium  be  small.  Thus, 
putting  1)  =  H  5  (RH)1  2  in  equation  (3.1  1 ),  we  find  that  the  thin  screen  approxima¬ 
tion  is  satisfied  if 

a  <  W/'i  '  ( 1 1/ R ) 1  2  (3.13) 

That  is,  the  deflection  of  the  ray  inside  the  refractive  medium  should  be  a  small 
fraction  of  a  scale  height  only.  The  integration  of  v0  along  the  actual  curved  raypath 
to  obtain  the  atmospheric  phase  perturbation  can  therefore  be  approximated  by  inte¬ 
gration  along  a  straight  line  path,  as  already  assumed. 

The  sceminglv  restrictive  condition  (3.13)  is  found  to  have  been  satisfied  for  most 
past  radio  occultations.  It  is  invariably  satisfied  in  stellar  occultations  where  the 
optical  rays  suffer  much  smaller  bending.  Using  a  direct  flux  conservation  argument 
not  restricted  by  the  thin  screen  approximation,  Eshleman  <7  al  ( I  979a)have  derived  a 
generalization  of  our  equation  (3.1 1)  that  reduces  to  this  result  when  the  condition 
(3.13)  is  satisfied. 


o 
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figure  3.2  Schematic  illustration  of  the 
changing  shape  of  the  atmo¬ 
spheric  Fresnel  zone  as  the 
observer  approaches  the  focal 
line  in  the  plane  of  the  sky 
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3.3  Wave-optical  effects 

Since  the  tree-space  Fresnel  zone  is  a  small  but  finite  fraction  of  atmospheric  scales 
of  variation,  the  true  phase  and  intensity  of  the  radio  wave  deviate  slightly  from  the 
geometrical  optics  results  of  the  previous  section.  Close  to  the  focus,  where  aFh 
increases  without  bound,  wave-optical  effects  eventually  become  dominant  for  any 
finite  value  of  the  wavenumber  k,  as  noted  previously.  VV’c  calculate  below  the  magni¬ 
tude  of  these  diffraction  effects  in  the  two  limits  where,  (i)  a  F  S  aj.  0,  but  with  no 
additional  restriction  on  aj.  v,  and  ( ii )  when  at  h  ie,  when  the  intensity  is 

critically  controlled  by  focusing  of  the  wave  around  the  curved  limb. 


The  former  of  these  cases  is  most  easily  treated  by  expanding  the  phase  factor  0 
beyond  the  quadratic  terms  and  retaining  the  lowcst-order  correction  to  the  electric 
field  E.  Writing 


0  -  [0y  (v ,  -y0)2  +  0zlU i  -z 0)2 


c0j(zi  -Zo)-1  +  c,2(z,  -z0)2(V| -y0)  + 


(3.14) 


it  is  found  that  the  leading,  finite  wavelength  correction  to  the  electric  field  involves 
terms  of  the  form  (z,  -z0)4  and  (z,  -z0)8.  Putting  K  =  E0(1+6E),  where  the  geomet¬ 
rical  optics  field  E0  =  0V,/20Z,/2  exp(i(kS — cot)]  and  6E  is  the  lowest-order  fractional 
correction  to  E,  we  find  after  considerable  algebra  that 


5E  =  -  -1-  <6(1  -0., 
32  z  z 


j32  D 
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(3.15) 


where  we  have  defined  j3  -  1/11  as  the  reciprocal  scale  height.  Thus  there  is  a  correc¬ 
tion  kSS  to  the  geometrical  optics  phase  by  an  amount 


d2I) 


k5S  =  “  32  0z(l  " ^  k 


(3.16) 


There  is  also  a  fractional  intensity  correction  given  by 
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9/J  k2 


(3.17) 


This  latter  result  is  functionally  similar  to  a  corresponding  result  by  Hubbard  (1979) 
for  propagation  in  two  dimensions. 


Since  the  atmospheric  phase  perturbation  kS0  =  kr’o^-  kilo,  it  follows  from  equa¬ 
tion  (3.16)  that  the  fractional  change  in  phase  due  to  diffraction  is 


k6S^  1  ,:daI>2 

kS0  32  z  k2 


(3.18) 


That  is,  there  arc  fractional  corrections  in  phase  and  intensity  that  arc  both  propor¬ 
tional  to  K^j34D2/k2  and  that  decrease  rapidly  with  occultation  depth.  Since 
x  spans  a  range  of  order  of  magnitude  10'4  10’8  for  past  radio  and  stellar  occulta- 

tions,  these  corrections  will  normally  be  altogether  negligible. 


Sufficiently  close  to  the  focal  line  all  =  R,  diffraction  effects  eventually  become 
dominant  for  any  Finite  value  of  k.  I.ctting  rx  denote  the  perpendicular  distance  of 
the  spacecraft  from  the  focal  line,  Hubbard  (1977)  has  shown  that  for  an  exponential 
atmosphere  the  intensity  close  to  the  foclis  is  given  by 
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(3.19) 
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where  J0  is  a  Bessel  function  of  zero  order.  The  focal  intensity  maximum  is  thus 

0max  =  47f2  RH/XD  (3-20) 

with  a  characteristic  width 

w  =  XD/R  (3.21) 

For  w  <<  rx  <  R,  the  average  intensity  reduces  to  the  geometrical  optics  result  (3.11) 
when  allowance  is  made  for  a  doubling  of  the  intensity  close  to  the  focal  line  where 
the  two  limbs  yield  nearly  identical  contributions  to  the  average  intensity.  Since  w  is 
typically  a  small  fraction  of  a  meter,  the  ray-optical  intensity  formula  (3.11)  is  valid 
extremely  close  to  the  focal  line.  As  shown  by  Eshleman  et  al  ( 1 979a),  the  results 
(3.20)  and  (3.21)  may  also  be  derived  from  simple  physical  arguments  relating  to  the 
limiting  distortion  of  the  atmospheric  Fresnel  zone  as  the  focal  line  is  approached. 


4  PROPAGATION  IN  A  TURBULENT  PLANETARY  ATMOSPHERE 


4.1  Phase  and  intensity  in  the  weak  scattering  approximation 


The  development  of  the  previous  section  provides  a  framework  for  calculating  the 
effect  of  turbulence  on  electromagnetic  wave  propagation.  In  the  thin  phase  screen 
approximation  turbulence  merely  adds  a  component 


Si  = 


«/•-> 


v\  (x.  y0.  z0)  dx 


(4.1) 


to  the  total  optical  path  length  at  the  exit  face  of  the  screen.  It  is  assumed  through¬ 
out  that  the  turbulent  component  of  refractivity,  t>, ,  has  zero  mean  and  satisfies 
02>1'2  <  tv  Equation  (3.1)  takes  the  form 


E(r2) 


// 


Jr 

2D 


2D2+(v2-yI)2  +  (z2-zI)2  +  2I)(S0+S1)l 


27ril)  -~ 


dr, 


(4.5 


In  the  weak  scattering  approximation  kiS ,  I  <  1.  and  the  integrand  in  equation  (4.2) 
may  be  expanded  around  the  point  of  stationary  phase  as  before.  Expanding  the 
electric  field  E  =  E0  +  E,  +  ....  where  F>0  =  (9(l  +  t'0).  E,  =Q  (v ,  ),  etc,  we  thus  obtain 


E,  (r2 )  =  k 0v0ze 


ik^ft/2 
2D  ' 


/  v\ (x,  v0,  Z0)dx  ffc 

-oo 


i2n,0v (y'  “y°)2+ (z* 


-ZoJ 


dr, 


(4.3) 


after  retaining  only  the  quadratic  terms  in  the  stationary  phase  expansion.  Let 


•  *.)=///  c 


i|uxx,  +  uv(y,  -y0)+uz(z,  -z0)J 


dZ(u) 


(4-4) 


be  the  Fourier-Stieltjcs  representation  of  ,  where  dZ(u)  =  Z'(u)du  (du  = 
duxduyduz),  with  ux,  Uy,  uz  being  the  spatial  wavenumber  components  in  the  x-,  v-, 
and  z-dircctions,  respectively.  The  Fourier  amplitude  Z'  of  i',  has  been  referenced  to 
the  stationary  phase  point  (y0,  z„).  For  a  homogeneous  and  isotropic  random  field 
<Z  (u,  )Z  ,l(u2)>  =  4>(U|  )  5(u,  u2 ),  where  ‘l’  is  the  spatial  power  spectrum  of  f,  and 
5(-)  is  the  Dirac  delta-function.  Combination  of  equations  (4.3)  and  (4.4)  yields  after 
algebraic  manipulations 

ik  C  °-  ^ 

Ei  (r2 )  k0yd/e  21)  )  Z'(u)  il{uL)  <iu  (4~r>) 
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Mere  sinc(-)  is  the  sine-function,  ux  =  (uy,  uz)  anti  we  have  defined 

'*  2k  l0yuy  f  0zuz  1  f  ‘  (0v‘  y<:)  +  0Z  z" )  -  '  ™l0v  y  I2  +  *i'  Z12 1 

“  _  "c  ~  "  ■  f  f  e  ~  dy'i  dzj 


(4.6) 


The  first-t>rder  perturbation  in  intensity  normalized  by  the  average  intensity,  0|  ,  and 
the  first-order  perturbation  in  phase  path,  S, ,  are  related  to  Kt  as  in 

0i  ~  2  Re[K0‘  K|  |,  S,  =  k'1  ItnlKj,1  K,  | 

respectively.  Since  K0  =  0y  0Z  exp(ik  c0o/-D),  we  obtain  after  integration  over  the 
phase  screen  in  equation  (4.6) 


~  «  u  V  , 

0,  =  2kV  /// sinc(  x)  Z (u)  sin  (P)  du 

2  tr 


(4.7] 


u  t  , 

S(  =  C  J/f  sinc(  x  )  /  (u)  cos  (P)  du 
2  tr 


(4.8) 


where  we  have  defined 

P  =  ft  (0yuy+0z.uz)  (4'9) 

for  the  mean  square  of  these  quantities  we  obtain,  assuming  1.0  <  V.  (1.0  is  the  outer 
scale  of  turbulence) 

<0?>  =  87Ttk*  //  <l»(ux)  sin2  (P)  du^  (4.10) 

<S?>  =  2ttC  /j«h(ux)  cos2(P)  dux  (4.11) 

I  he  quantity  <0j>  is  called  the  scintillation  index,  the  square  root  of  which  is  often 
termed  the  modulation  index. 


liquations  (4.10)  and  (4.11)  are  the  appropriate  generalizations  of  the  weak  scatter¬ 
ing,  mean  square  fluctuations  in  intensity  and  phase  path  to  turbulence  superimposed 
on  an  inhomogeneous  background.  (For  stellar  oooultations  the  intensity  should  be 
convolved  with  the  projected  stellar  brightness  distribution  in  the  occulting  atmo¬ 
sphere,  sec  section  4.8.)  The  effect  of  the  ambient  atmosphere  is  seen  to  appear  onlv 
through  the  quantities  0V  and  0Z,  whose  interpretation  in  terms  of  the  mean  signal 
intensity  was  established  In  the  previous  section.  A  sufficient  condition  for  neglecting 
coupling  it)  the  inhomogeneous  background  atmosphere  is  therefore  that  0V,  0Z  both 
be  essentially  unity,  or  equivalently,  that  the  change  in  intensity  caused  by  the 
ambient  atmosphere  alone  be  negligible.  1'his  condition  is  not  satisfied  in  cither  radio 
or  stellar  occultation  experiments,  however,  implying  that  interpretation  of  weak 
scintillation  data  in  such  experiments  should  proceed  from  the  generalized  formulas 
derived  here.  Surveillance  of  the  earth's  upper  atmosphere  using  a  pair  of  satellites 
may  also  involve  sufficient  ray  bending  for  weak  scintillation  data  to  be  analyzed  in 
terms  of  the  results  given  here,  especially  if  the  raypath  traverses  regions  of  thermal 
inversion  where  the  local  value  of  i)‘  v0/dz? ,  contained  in  the  definition  of  0Z>  may 
attain  appreciable  values.  It  may  be  noted  that  although  the  previous  derivation  has 
assumed  e,  to  be  a  strictly  homogeneous  anti  isotropic  random  variable,  the  results 
(4.10)  and  (4.11)  arc  still  valid  if  n,  is  inhomogeneous,  provided  that  the  scales  on 
which  the  inhomogcneitics  occur  arc  large  compared  to  any  other  scale  of  interest, 
i  e,  to  1 , o  and  the  Fresnel  scales  a(,  h,  a,  . 

Having  established  the  first-order,  weak  scattering  results  for  phase  and  intensity 
appropriate  to  radio  and  stellar  occupations,  it  is  important  to  explore  their  validity. 
For  the  ambient  atmosphere  the  thin  screen  approximation  was  found  to  require 


Z3 


a«(H/R)12.  Since  the  first-order  intensity  and  phase  perturbations  are  both  func¬ 
tions  of  0V<  0Z,  the  condition  on  the  smallness  of  a  also  applies  in  this  case.  1  he  thin 
phase  screen  approximation  also  requires  that  the  scintillations  in  intensity  be  small 
inside  the  screen.  This  condition  is  automatically  satisfied,  however,  since  we  shall 
prove  in  section  4. 3  that  <0J>  is  small  compared  to  unit)  at  any  distance  from  the 
screen  when  kJS1(r,)|<  1,  i  e,  when  the  turbulent  phase  perturbation  at  exit  from 
the  screen  is  small.  W'e  observe,  finally,  that  since  C  <  1)  the  quantities  0V  and  <t>z  do 
not  change  appreciably  through  the  screen  and  may  hence  be  associated  with  a  fixed 
spacecraft-limb  distance  x  =  1). 


The  stationary  phase  expansion  also  imposes  certain  restrictions  on  the  validity  of  the 
previous  results.  These  limitations  are  most  easily  appreciated  by  noting  from  the 
derivation  of  equation  (4.5)  that  if  Ay  =  y,  y0  and  A/.;  z,  -z0,  a  Fourier  compo¬ 
nent  of  wavenumber  ux  =  (uv,uz)  gives  a  contribution  proportional  to 


Jy>  l<\'  I  Ay  t  ^  0yuy)2  +  <t>:)  (Az  +  ^  0zuz)2  ] 


iAv  dAz 


(4.12) 


to  the  perturbed  electric  field  at  the  receiver.  Defining  uv  =  2tr/tv  and  uy  ~  it  is 

seen  front  this  expression  that  an  eddy  of  wavenumber  components  uv,  uz  gives  a 
contribution  to  tin.  perturbed  electric  field  that  originates  from  a  position 
( i A\ ,  0vM)Tv,  ,Az!  =  <p7  Xl)/t.  )  relative  to  the  stationary  phase  point.  Such  ele¬ 
mentary  wave  contributions'  are  tints  scattered  through  angles  0V  -  0vX/Cv  and  d7  = 
Qy  A,fz  in  the  \\  ami  xz  planes,  respectively.  When  0V  -  0Z  =  1,  this  is  essentially 
the  Bragg  condition  for  diffraction  by  a  sinusoidal  spatial  diffraction  grating.  Now. 
for  Kolmogorov  turbulence  <S;>  and  <02>  may  be  shown  to  derive  their  principal 
contributions  from  eddies  of  dimension  (f'v,  tz)  ~  L0  and  (Cv,  Sz)~  (ap  j,,  ap.v),  res- 
pcctively.  f  he  quadratic  approximation  <if  the  stationary  phase  expansion  should 
therefore  be  validated  to  distances  Ay  ~0vXI)/I,o.  ,A/,i  ~  0Z  Xl)/L0  and  |Ay|  ~ 
,q  h  ,  i A/. i  ,,  for  these  two  quantities,  respectively.  F. valuation  of  the  next-order 

terms  in  the  stationary  phase  expansion  indeed  turn  out  to  be  small  at  these  distances 
if  H  I,„  a j.  v  ,  and  provided  also  that  0V  <  <  I-0 /a  j.- _() )  (L0  R/a f.-  ) 1 7  in  the  evalua¬ 

tion  of  <Sf  -.and  <  (Rill1  2  /'  a  j.  (l  in  the  ev.duation  of  <0i>.  In  shallow  occulta- 
lions  i  hese  constraints  are  alwav  s  met .  but  very  close  to  the  focal  line  both  conditions  arc 
violated  and  neither  ^$2>  nor  <S2>  will  be  given  by  simple  formulas  of  the  type 
14.10)  and  (4.1  1). 

We  are  now  prepared  to  disc  uss  another  quantity  of  principal  interest  in  radio  occul- 
tation  experiments:  the  ray  bending  angle.  It  is  readily  shown  from  the  previous 
equations  that  a  major  component  of  the  mean  square  perturbation  in  bending  angle 
involves  an  integral  of  the  form 

f  f  4>l  iq  )  u2  cos2  ( I’)  duz  (4.13) 

For  Kolmogorov  turbulence  this  integral  diverges  as  a  function  of  the  upper  limit,  and 
is  thus  controlled  bv  the  smallest  eddies  of  dimension  -  li0.  the  inner  scale  of  the 
turbulence.  Bv  relation  (4.12),  eddies  of  this  size  yield  scattered  waves  that  originate 
from  distances  d  --  XD/t0  from  the  stationary  phase  point.  With  C0  1  cm.  d  is 
commensurable  with  the  dimensions  of  the  planet.  The  implication  is  therefore  that  a 
calculation  of  the  perturbation  in  ray  bending  angle  requires  the  complete  spatial 
extent  of  the  scattering  medium  to  be  stipulated,  and  the  stationary  phase  expansion 
cannot  be  truncated  after  the  second  order  terms.  The  problem  is  compounded  by 
the  fact  that  in  going  from  equation  (3.1)  to  (3.21,  me  exact  kernel 


exp(|k|r?  r,ij 

1*2  *1  I 
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of  equation  (3.1)  was  approximated  by  the  Fresnel  kernel 


—  exp( i  —  p2 
D  1  2D 


which  requires  that  x  k  p"* / O'1  1,  where  p  -3  [(y,  y,  ):  *  (/.,  I  ]1  ;  is  the  lateral 

separation  of  the  receiver  and  an  arbitrary  point  on  the  exit  face  of  the  phase  screen 
( Tatarskii  1971,  pp  223).  For  the  perturbations  in  intensity  and  phase  this  condition 
is  strongly  satisfied,  since  for  these  quantities  p  (XI))1  2  and  <  Al)/L0,  respectively. 
However,  the  fluctuations  it;  ray  bending  angle  derive  their  principal  contribution 
from  waves  scattered  from  distances  f  from  the  stationary  phase  point,  where  fc  is  the 
effective  lateral  extent  of  the  turbulent  medium.  For  turbulence  that  varies  exponen¬ 
tially  with  altitude  on  a  scale  ilt  (cf  section  4.2),  the  medium  extends  effectively  to  a 
distance  p  =  t  a  (RHt)‘  2  parallel  to  the  limb.  Assuming  Ilt  -  II  =  10  km.  we  find  that 
for  X-  and  S-baml  occulta! ions  of  f'S  spacecraft  typical  values  of  x  an  such  that 
x  S*  1 .  Hence  a  numerical  approach  is  apparently  required  to  obtain  the  mean  square 
fluctuations  in  bending  angle 


Most  previous  analyses  of  scintillation  data  from  radio  and  stellar  occult  ition  experi¬ 
ments  have  ignored  the  coupling  effect  between  the  turbulence  and  the  ambient 
background  described  and  calculated  here  (Gurvich  1 969,  Golitsyn  and  Gurvich  1971, 
Woo  and  Ishimaru  1973,  Woo  ct  it/  1974.  Woo  and  Ishimaru  1974,  Woo  1973a, 
Jokipii  and  Hubbard  19771.  The  first  attempt  to  account  for  the  coupling  effect  on 

the  scintillation  index  was  apparently  made  by  Woo  eta/  (1973).  Later  and  indepen¬ 
dently  the  problem  was  treated  by  Haugstad  (1977,  1978c, e).  by  Hubbard  ct  al 

(1978)  and  recently  also  by  Woo  ct  al  (1980).  An  interesting  heuristic  approach  to 

the  problem  was  taken  by  A  Young  (1976).  I  lie  problem  of  calculating  the  mean 
square  fluctuations  in  phase  and  bending  angle  does  not  appear  to  have  been  addres¬ 
sed  by  other  authors. 


4.2  An  atmospheric  model 

F.xplicit  evaluation  of  the  intensity  and  phase  scintillations  requires  a  model  to  be 
stipulated  for  both  u0  and  u,  .  1  lie  former  quantity  has  already  been  assumed  to 
correspond  to  a  spherically  symmetric  and  exponential  atmosphere,  for  which 

=  Gilr,, )  e\p|  p'lr  r„  t|  (4.14) 

where  ^3' 1  =  II  is  the  scale  height,  and  r  r  r  is  measured  from  the  center  of  the  planet. 
For  the  turbulence  component  we  assume 

»q  (r)  =  <t-,;tr)>1 2t'?fr)  (4.13a) 

<v  j  (r)>‘ 2  -  <tq  (r0  )> 1  *  exp|  d,(r  fall  (4.13b) 


I  hat  is,  tq  is  the  produc  t  of  a  spherically  symmetric  component  that  varies  exponen¬ 
tially  with  height  on  a  scale  dj1  ■  Ht.  and  a  strictly  homogeneous  and  isotropic 
random  component  u*  of  zero  mean  and  unit  variance.  1  he  power  spectrum  ‘1>  of  tq 
will  be  assumed  to  be  of  the  Yon  Harman  type,  l  e. 


■I>(u)  - 


F(p;2i  I  .a  exp(  f,2,u‘  I 
TT1-!'!  \  |  1  •  I  f,  U2  I1' 


(4.16) 


Here  u  =  I  uj.  *-  u2 >  u^t 1  2  is  the  three-dimensional  sc  alar  w  ave  number.  I.u  and  f!„  are 
the  outer  and  inner  scale  of  the  turbulence,  respectively,  and  p  is  a  parameter  control¬ 
ling  the  slope  of  <I»  m  the  inertial  subrange,  where  I u  •'  1  .  for  Kolmogorov 

turbulence  p  ~  I  1  / 3. 
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For  purposes  for  numerical  evaluation  also  the  magnitude  of  <v]>  should  he  stipula¬ 
ted.  Turbulent  fluctuations  in  refractivity  may  be  related  to  those  in  temperature  and 
pressure  using  the  relation  (e  g,  Tatarskii  1971) 

v  =  7.9  •  1  O'5  p/T  (4.17) 


which  assumes  that  the  refractivity  from  any  vapors  present  can  be  neglected.  Small 
changes  in  refractivity,  V\  ,  are  thus  related  to  similar  changes  p,  in  pressure  and  I  , 
in  temperature  according  to 

m  =  7.9  ■  1  O'*  |  P-4  P;  T.  |  (4.18) 

i  {  i :  1 

For  the  centimeter  to  kilometer  eddy  sizes  of  interest  here  the  temperature  relaxation 
is  slow  enough  to  suggest  an  adiabatic  relationship  between  pi  and  lj  (Lawrence  and 
Strohbchn  1970),  so  that  p,/p  =  7/(7  - 1  )•' T,  / T,  where  7  is  the  ratio  of  specific  heats 
at  constant  pressure  and  volume.  With  this  assumption  equations  (4.1  7)  and  (4.18) 
combine  to  give 


<v.> 


(7-  1) 


•  2  <’Ij_> 

i  o 


(4.19) 


For  an  isothermal  ambient  atmosphere  the  fractional  mean  square  fluctuations  in 
refractivity  thus  vary  with  altitude  as  the  variance  of  the  temperature  fluctuations.  An 
estimate  of  the  latter  quantity  may  be  obtained  from  the  temperature  structure 
function 


1),  ( Ar)  -  C-j- lArl2'*3  (4.20) 

valid  for  Kolmogorov  turbulence  and  separations  K0  ^  I  An  <  L0.  An  approximate 
value  of  <  1;>  is  obtained  from  this  relation  by  putting  Ar  =  l.„  and  noting  from  the 
definition  of  I.„  that  l)t  (1,0)  2  2  <  1]>.  Thus* 

i>,(i.o)  =  hi  i -o'  '  2<t-;> 

The  fractional  fluctuations  in  refractivity  are  thus  approximately 


Willies  for  the  outer  scale  may  be  obtained  by  direct  measurements,  or  its  value  may 
be  inferred  from  the  bulk  flov\  properties.  The  temperature  structure  constant,  CJ  , 
may  be  estimated  from  the  mean  How  characteristics  for  various  types  of  turbulence 
(see,  e  g,  Tatarskii  1971,  pp  72  74).  Alternatively  one  might  use  the  estimate 

2  2- 1  O' 1  K;  tn'2  1  corresponding  to  an  altitude  level  of  '  17)  km  in  the  earth’s  atmo¬ 
sphere  (button  etui  1972).  Assuming  I.„  HI1  m  and  T„  -  200  K  we  obtain 

’  '  "I 

This  estimate  may  be  appropriate  for  Jupiter.  For  the  atmosphere  of  Venus  the 

higher  ambient  temperature  (and  perhaps  also  the  smaller  value  of  L0 )  probably 
makes  equation  (4.22)  an  over-estimate  of  the  turbulence  strength  in  this  case.  As  we 
shall  see  m  set  (ion  7,  the  best  wav  to  estimate  •-  v\  >  is  probably  bv  direct  analysis  of 
radio  or  stellar  scintillation  measurements.  Lacking  such  direct  information  on  the 
turbulence  strength,  we  shall  explore  in  a  forthcoming  section  the  implications  of 
turbulence  of  the  relative  strength  given  bv  relation  (4.22). 


*  tor  Kolmogorov  lurbiiloii  i-  I  In-  more  ■  om<  I  rrl.itionvhip  i\  given  bv  '  !.'U<I  j  - 


(4.22) 


4.3  Radio  and  stellar  intensity  scintillations 


The  turbulence  model  described  in  the  previous  section  may  now  be  combined  with 
the  scintillation  formula  (4. 10)  to  yield  an  explicit  result  for  the  radio  scintillation 
index,  [-'.valuation  of  this  equation  by  transforming  to  polar  coordinates  yields  for 
3  <  p  <  6  (Haugstad  197 8c) 


4  'a  a~-3\VR(0v,0z;p)<^>  (4.23a) 

^  JF.()  JK,0  K 


where 


l> 


7T  —  -  1 

Wr(0\.0z;p)  =  /  (0ycos20  +  0zsin2O)“  d0 

-7T 


(4.23b) 


In  deriving  this  result  we  have  assumed  that  L0  >  a^  i,  >aj.-  v,  implying  that  the 
major  contribution  to  the  scintillations  is  from  turbulent  eddy  sizes  in  the  inertial 
subrange.  In  its  present  form  this  result  is  therefore  not  applicable  to  very  deep  occulta- 
tions,  where  a  j.  h  may  exceed  l.().  We  shall  consider  this  case  separately  in  section  6. 

In  ecpiation  (4.23)  coupling  between  the  turbulence  and  the  inhomogeneous  back¬ 
ground  appears  through  WR  only.  For  a  strictly  homogeneous  background,  or  when 
the  differential  ray  bending  is  sufficiently  small  that  0V  0Z  S  1.  \\'  =  2tr  for  anv 
value  of  the  spectral  slope  parameter  p.  In  this  case,  alter  manipulations  invoking  the 
gamma  function,  equation  (4.23a)  may  be  shown  to  coincide  with  the  classic  Rvtov 
result  in  the  thin  screen  approximation. 

If  we  now  consider  occultations  for  which  the  approximation  0V  ?  1  can  be  made, 
but  with  no  restriction  on  <py,  evaluation  of  equation  (4.23b)  gives 

r  2tr,  <pz  =  1 

'VR(1.0z;p)=  J  „  ,  (4.24) 

1  r-s  ( L’.-.I ) 


Since  for  p  =  11/3  the  lower  va'ue  corresponding  to  0y<  1  is  about  0.6  times  the  upper 
value,  it  follows  that,  (i)  coupling  to  the  inhomogeneous  background  reduces  the 
scintillation  index  by  less  than  a  factor  of  two,  and  (ii)  this  reduction  is  essentially 
independent  of  0  s  0Z  whenever  defocusing  reduces  the  average  intensity  by  a  few  dB  or 
more.  I  bis  result  disagrees  with  a  prediction  by  A  Young  (1976),  who  concluded  from  a 
heuristic  calculation  that  the  scintillation  index  should  increase  by  a  factor  0' 23  as 
deeper  atmospheric  regions  were  probed.  Although  neither  the  works  of  Hubbard  ct  al 
(1978)  or  Woo  et  al  (1980)  gave  explicit  results  for  the  coupling  effect,  their  expres¬ 
sions  for  the  scintillation  index  may  be  manipulated  to  yield  a  similar  result  to  that 
derived  here.  We  shall  see  in  the  next  section  that  the  rather  small  effect  of  the 
inhomogeneous  background  on  the  scintillation  index  does  result,  however,  in  sub¬ 
stantial  changes  in  the  shape  of  the  radio  and  stellar  scintillation  power  spectra. 

lor  stellar  occultations  the  point  source  intensity  correlation  function  should  be 
convolved  with  the  normalized  stellar  brightness  distribution  as  projected  above  the 
limb  of  the  planet.  In  the  spectral  domain  this  operation  yields  the  general  expression 
(c  f,  Hubbard  et  al  1978) 


/?  'IMUj.)  i0(UjJiduj_ 

0i 

for  the  scintillation  index  of  a  stellar  source.  In  this  expression  (ujJ  is  the  spatial 
power  spectrum  of  <$j>  (i  c,  the  integrand  of  equation  (4.10)),  and 
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when-  ;tr)  is  the  normalized  brightness  distribution  of  the  projected  stellar  disk.  For  a 
uniformh  illuminaied  stellar  disk  of  projected  radius  as  one  obtains 


ivMuj!2 


2J.(  asV,J 

Vl 


where  |,  is  a  Bessel  function  of  order  one.  Instead  of  the  point  source  result  (4.10) 
we  thus  obtain 


-  di'>  -  HrrVk sin2  (I* 


a,  u 


sul 


(4.23) 


lhe  additional  factor  in  the  integrand  is  due  to  averaging  over  the  stellar  disk  and 
acts  to  suppress  contributions  to  the  scintillation  index  from  eddy  si/es  smaller  than 
the  projected  stellar  disk.  A  closelv  similar  result  to  equation  (4.23)  has  been  obtained 
bv  Hubbard  et  <il  1 19781. 


If  the  double  inequality  l.„  '■  as  aj.  is  satisfied,  evaluation  of  equation  (4.25) 
yields  for  3  <.  p  <.  fi 


l  (  — — *-  )  .  -J> 

i  2-  I  °)  ^-Ws(0v,0z;p)<rj>  (4.26a) 

|',P  =  )  *'s  as 


where  in  this  case 

U’s(0v,(>/;p)  /  (0vcos2d  *  0., sin 2 0 ) 2 <10 

•  77 


(4.26b) 


For  intermediate  and  small  values  of  the  ratio  as/ay  ,  numerical  integration  of  equa¬ 
tion  (4.23)  is  required.  For  shallow  occultations  equation  (4.26b)  gives 


2tr,  i,f  1 

.,s  l':(^/2l 
1(3)  '  z 


(4.27) 


Since  the  lower  value  for  0y  <-  1  is  "  1.76,  we  conclude  that  coupling  to  the  inhomo¬ 
geneous  background  reduces  the  stellar  scintillation  index  by  at  most  a  factor  s  8.6 
under  the  previous  assumptions.  A  Young’s  prediction  for  this  case  was  an  increase  in 
the  scintillation  index  by  a  factor  Oj,4  ',  attending  a  compression  of  the  stellar  image 
by  a  factor  0/  perpendicular  to  the  planetary  limb. 


It  is  instructive  to  express  the  scintillation  index  as  a  function  of  occultation  depth 
under  the  assumption  that  <f \>/Vq  is  constant.  To  this  end  we  first  note  that 


IF’ m2  <»>i> 

2trKl)2  v„ 


We  have  here  used  the  fact  t hat  a  -  (2jrH/l I)1  2 v0  for  bending  angles  tv  <  (H/R) 1  2 
and  M  «1)/II  as  before.  Fliminating  <tq>  between  this  relation  and  cither  of  equa¬ 
tions  (4.23a)  or  (4.26a)  yields  for  ~  I 


'  <>\  > 


K(i,)l,aw)'"t’(^)' 


V  <<’i> 
R  i’u 


V\V(l.0;p) 

<> 


(4.28) 


For  radio  occultations  W  \VRand  a  r  a,.  while  for  stellar  occultations  \V  =  \VS  anil 
a  av  l  he  numerii  al  factor  K  is  the  co-factor  of  equations  (4.23a)  and  (4.26a),  rrs- 
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Figure  4.1  Variation  with  occultation 
depth  <f>  f  of  the  t/uantities 
FI  and  FI  fur  radio  occul¬ 
ta  t  ions 


Figure  4.2  Variation  with  occultation 
depth  (f)'1  of  the  ipiantities 
S2  and  ft  for  stellar  occul- 
t  at  ions 


1  he  solid  curves  display  £2  and  FI  when  coupling  to  the  inhomogeneous  background  is  properly 
accounted  for,  while  the  broken  curves  are  calculated  for  a  model  where  this  effect  is  not  included.  1  he 
curves  are  normalized  such  that  FI  in  the  latter  case  approaches  unity  as  0 1  -* 


pectivcly.  If  the  spacecraft  (or  the  earth  for  stellar  occultations)  is  moving  perpendi¬ 
cular  to  the  incoming  ray  asymptote,  1)  is  approximately  constant  for  small  a,  and  all 
variation  of  <0:,>  with  depth  is  contained  in  the  two  last  factors  of  equation  (4.28). 
Figures  4.1  and  4.2  show  the  variation  with  occultation  depth  of  the  quantity 


ft  =  (  1-0 W(  1,0:1 1/3) 

0  2n 


(4.29) 


for  both  radio  and  stellar  occultations.  along  with  the  result  that  would  have  been 
obtained  if  coupling  to  the  inhomogeneous  background  were  falsely  neglected.  Note 
that  for  constant  1)  and  11,  M  is  proportional  to  atmospheric  density  or  pressure  at  the 
ray  periapsis.  Also  shown  in  the  figures  is  the  corresponding  variation  of  the  absolute 
mean  square  fluctuations  in  intensity,  which  are  proportional  to  FL  ~  <f>:  Si.  In  con¬ 
trast  to  the  scintillation  index  which  increases  rapidly  with  depth,  the  absolute  inten- 
sity  scintillations,  <0J>  =0'  <0}>,  approaches  a  constant  determined  by  the  asymp¬ 
totic  value  of  W  as  0  -*■  0. 


We  conclude  this  section  by  proving  that  in  the  weak  scattering  regime  <$]  >  <  1  also 
when  coupling  to  the  inhomogeneous  background  is  accounted  for.  For  turbulence 
superimposed  on  a  homogeneous  background  this  result  has  previously  been  proved 
by  Salpeter  (19ti7)  and  others.  Consider  the  mean  square  fluctuations  in  optical  path 
at  exit  from  the  phase  screen.  This  quantity  may  be  obtained  from  equation  (4.1 1 ) 
by  putting  P  =  0  and  evaluating  the  double  integral  over  the  power  spectrum,  yielding 


<S?  ( r |  )> 


l’(|,,)  )(p  2) 


(4.30) 


Eliminating  <fj>  between  this  result  and  equation  (4.23a)  gives 

<fi>  =  K(p)  (••K1..)'-‘VVR^y.<-z4>)  |ki  <  s?(  )>, 

l-o  27T 


(4.31) 


whore  K|[>l  ( 27^)-'^>',-  |see  ftp  '  h/|  4l'(p/2)|  is  of  order  of  magnitude  unity  for  |>  -  11/3. 

Since  a  j.  „  <|  1.0  per  assumy  ion,  and  we  have  shown  that\VR<  27T,  it  follows  that 
1  in  the  weak  s  .ttering  approximation  where  k<Sj(r,)>17  <  1.  Kvcn 
though  the  weak  scattering  condition  has  been  used  to  derive  the  previous  results,  it 
is  likely  that  only  the  weaker  condition  <<t>\  >  1  is  necessary,  however.  For  turbu¬ 

lence  superimposed  on  a  homogeneous  background  this  has  been  proved  by  Tatarskii 
(1971,  pp  23t>).  One-dimensional,  numerical  computations  of  intensity  scintillations 
superimposed  on  an  inhomogeneous  background  indicate  that  <0\  >  <  1  is  indeed 
the  proper  criterion  also  in  this  case  (French  et  al  1980). 

4  4  Temporal  power  spectra 

4.4.1  lhcorctic.il  phase  and  intensity  power  spectra 

I  he  scintillation  index  provides  only  an  averaged  representation  of  the  effect  of  the 
inhomogeneous  background  on  radio  and  optical  wave  propagation  in  turbulent  plane¬ 
tary  atmospheres.  More  detailed  information  is  contained  in  the  temporal  power 
spectrum,  the  shape  of  which  will  be  shown  here  to  depend  critically  on  both  occul- 
tation  depth  and  geometry  when  coupling  to  the  inhomogeneous  background  is 
accounted  for. 

1  he  phase  autocovariance  (unction  in  the  plane  \  -  1)  of  the  observer  is  given  bv 

B  =  h*  <S ,  (l),y  2 )S ,  ( 1 ) ,y  j  4  Ayo.z,  +  Az,)>  (4.32) 

Geometrical  optics  energy  conservation  implies  that  separations  Ay  2 ,  Az2  at  the 
receiver  project  into  ray  separations  Ay,  -  0vAv2  and  A/.,  =  </>zAz2  at  ray  periapsis. 
Velocity  components  vv.  vz  of  the  spacecraft  are  therefore  related  to  corresponding 
components  v*  and  v*  .it  the  ray  periapsis  according  to 

v*  =  0yvy,  v*  =  <t>v\7  (4.33) 

It  is  customary  to  adopt  the  "frozen-field"  hypothesis,  whereby  turbulent  inhomo¬ 
geneities  are  assumed  to  be  conceded  across  the  line  of  sight  at  constant  velocity, 
which  for  the  present  case  equals  the  relative  velocity  between  that  of  the  rav  at 
periapsis  and  the  atmospheric  bulk  velocity,  due  to  either  planetary  rotation  or  the 
presence  of  atmospheric  winds,  or  both.  We  shall  assume  that  the  components  of  the 
latter  are  small  compared  to  v*.  v*.  In  this  case  ray  separations  Ay,,  Az,  at  the 
periapsis  are  related  to  time  separations  r  according  to  Ay,  =  v*T,  Az,  =  vzr.  Using 
relations  (4.33)  and  drawing  on  the  analysis  of  section  4.1,  we  find  that 

B  ~  k:  fl  <*’s|  (ux)  exp|i(uy0yVy  +  uz0zvz)r|duj_  (4.34a) 

where 

•I’s,  (uLt  -  2jtV‘I>(ul)  cos’  ( P)  (4.34b) 

is  the  spatial  power  spectrum  of  the  scintillations  in  phase  path,  as  given  implicitly  bv 
equation  (4. 1  1 ).  Fourier  transformation  over  r  yields  for  the  temporal  power 
spectrum  of  the  phase  scintillations 

W,  (w)  -  2trk‘  //  ll’S|  (ux)  5[uv0yV  f  u z</>zv7  co|dux  (4.33) 

Using  the  turbulence  power  spectrum  (4.16)  this  integral  may  be  evaluated  by  trans¬ 
forming  to  polar  coordinates  and  performing  the  integration  over  the  angular  variable. 


After  a  suitable  change  of  integration  variable  in  the  remaining  integral  we  find  that 
(Haugslad  1979a) 

■>  ^  \eos:|(-~)2|0yCos:(flj  +  6)  +0/,sin2(dj  +  6!  jj 

U’|  (w)  '  « ,  cj  11  f - — -  -  — - dx  (4.3b) 

1  '  I  |,  MCJ\  )2|P/-,(XJ  1)12 

CJ„ 

We  have  here  defined  0t  -  0,  =■  ros'1  (l/\),  to,,  =  toc  =  2ir1'2  \*/uy  0  , 

6  -  tan'1  (v^Vy),  and  k,  -  | 4tt  1  -  Hp/LM/Rf p-3)/2 ) |k'k:  I.j’Of >/v*.  We  observe  from 
this  result  that  the  shape  of  the  phase  power  spectrum  is  controlled  by  the  four 
parameters  cou,  uJ(.,  0v  and  0Z.  In  view  of  their  definition,  cu0  and  t oc  are  the 
characteristic  frequencies  assoc  iated  with  the  ra\  speed  at  petiapsis  and  the  outer 
scale  and  the  free  space  Fresnel  /.one,  respect ivclv. 

I  (trough  entirely  analogous  procedures  the  radio  and  stellar  intensitv  scintillation 
power  spectra  are  obtained  as 

xsin:  {( -  ~  );  |0vcos2  (dj+6 1  +  fusin'  ( Oj  t  6)|j> 

W;(oj)  -  K:^>  £  7  -  •  '  <K  I  4.4  7 1 

1  *  i  iif  i’ r-ix-  D,: 


x  ’  [0,  cos' (fb  t  6)  *  0  siir'R:  *  6)|" 

W,M  =  1  f  —  -  —  —  |‘[  (  wx  1  tlx  <l.38| 

n  i  1 1  +  ( ):  |l>/'(x:  1)'-  wc 

w.l 

respectively.  For  the  former  result  k :  -  (  Kiir-  l'(p/2)/l  ((p-A)/‘— ))  fek"  I <2/7 >/\ **’ ;  0 ,  , 
02  and  oj„.  ojc  arc  defined  .is  before.  The  stellar  intensity  scintillation  power  spect¬ 
rum  is  subjected  to  the  same  limitations  as  the  expression  (4. lib)  for  the  scintillation 
index.  In  this  case  u,  v*-'as,  and  k  ,  ( 1  (itr'  I ' ( p / 12 )/T((p-3)/L!))  I):l.i’,  <tq >/( as\ *A  I. 

We  note  that  all  of  the  power  spec  tra  calculated  here  are  valid  onlv  for  frequencies 
much  larger  than  the  characteristic  frequency  \"7H  associated  with  either  of  the  scale 
heights  II  or  1 1( . 

Figures  4.3  to  4.9  show  the  phase  and  intensity  power  spectra  as  a  function  of 
normalized  angular  frequency  w/u>(.  or  ui/oj,, .  for  different  values  of  m  v*/v*  and  0. 
assuming  0V  1.  In  all  cases  computed  we  have  assumed  u>(./cl!0  -  3.  For  the  phase 
spectra  coupling  to  the  inhomogeneous  background  is  seen  to  reduce  the  fringes  in 
the  spec  trum  for  w  >  oJ(  in  a  central  occuitation  Int  0)  as  deeper  atmospheric- 
regions  are  probed.  Figure  4.7)  displays  the  change  in  spectral  shape  from  a  central  to 
a  more  nearly  grazing  occuitation,  for  a  fixed  value  0  :  0.1  of  the  mean  signal  inten¬ 
sity. 

by  contrast  to  this  behavior,  the  intensity  power  spectra  display  a  strong  dependence 
on  both  occuitation  geometry  (tn)  and  depth  (0).  According  to  Figure  1.6.  Un¬ 
characteristic  fringes  in  the  power  spectrum  above  the  Fresnel  frequency  u>(.  alto¬ 
gether  disappear  as  deeper  regions  are  probed  in  .1  central  occuitation.  In  a  grazing 
occuitation.  however,  these  fringes  become  increasingly  deeper  and  follow  a  characte¬ 
ristic  decrease  in  scintillation  power  below  u)(.  as  denser  regions  are  probed. 
Figure  4.8  shows  the  progressive-  change  in  spectral  shape  from  .1  c  entral  to  a  more 
nearly  grazing  occuitation  for  a  fixed  occuitation  depth  corresponding  to  0  0.1. 

With  only  minor  modifications  the  changes  in  spec  tral  shape  with  m  and  0  for  a  point 
source  are  reproduced  also  for  an  extended  source  when  the  projected  stellar  radius 
substantially  exceeds  the-  free-spare  Fresnel  zone.  Figure  4.9  shows  the  change  in 
spectral  shape  with  occuitation  depth  for  a  grazing  occuitation.  If  the  projec  ted  stellar 
radius  is  small  compared  to  the-  free-space  Fresnel  zone,  averaging  over  the  stellar  disk 
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Figure  4.3  Phase  scintillation  power  spec¬ 
trum  at  various  depths  0  in  a 
central  oeeultation.  assuming 


Figure  4.4  Phase  scintillation  power  spec¬ 
trum  at  various  depths  0  in  a 
grazing  oeeultation,  assuming 
ojk./co0=3 
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Figure  4.5  Phase  scintillation  /lower  spec 
tram  for  different  occupation 
geometries  and  a  fixed  depth 
Q  0.  /.  assuming  oJ(./ci7()  ~  J 


Figure  4.6  Point  source  intensity  scintilla¬ 
tion  power  spectrum  at  vari¬ 
ous  depths  0  in  a  central  oc¬ 
cupation,  assuming  ojc/oj0 -  3 
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will  yield  no  substantial  effect  on  either 
the  scintillation  index  or  the  power 
spectrum.  In  this  case  power  spectra 
will  correspond  closely  to  the  point 
source  scintillation  spectra  of  Fig¬ 
ures  4.6  -4.8. 


0  =  0.33 


The  classical  Rytov  results  for  the 
phase  and  intensity  power  spectra  cal¬ 
culated  for  turbulence  superimposed  on 
a  homogeneous  background  correspond 
to  the  cases  0  =  1.0  in  Figures  4.3,  4.6 
and  4.9. 


6  =  0.1 


0  =  0.01 


W3(w) 


Figure  4.9  Stellar  scintillation  power  spec¬ 
trum  for  various  occupation 
depths  0  in  a  grazing  occupa¬ 
tion,  assuming  co(/u>0  =3  and 
n\/ah  .0  *  1 


4.4.2  Observed  power  spectra 

Power  spectra  calculated  from  the  Mari¬ 
ner  5  and  10  Venus  encounters  (Woo 
et  al  1974,  Woo  1975a,  Woo  et  al 
1980)  and  from  the  Pioneer-Venus 
entry  probes  (Woo  et  al  1979)  have 
been  reported  in  the  literature.  For  mis¬ 
sions  by  other  planetary  spacecraft  and 
for  stellar  occultations  no  correspond¬ 
ing  results  have  been  published. 

Figures  4.10a-!>,  taken  from  Woo  et  al 
(1974),  show  the  recorded  part  of  the 
Mariner  5  S-band  entrance  and  exit 
occultation  as  a  function  of  universal 
time  and  altitude  above  the  mean 
Venusian  surface.  Regions  marked  A  to 
1)  were  suspected  of  containing  struc¬ 
tures  due  to  turbulence,  and  were 
selected  for  further  analysis.  After  re¬ 
moval  of  a  quadratic  trend  function 
from  the  exit  data,  Woo  et  al  (1974) 
calculated  the  temporal  power  spectra 
and  compared  with  predictions  from  a 
layered  turbulence  model  fitted  to  stan¬ 
dard  weak  scintillation  formulas, 
assuming  a  priori  that  F0  <  ap  ().  In  a 
later  comparison  Woo  ( 1  975a)  claimed 
improved  theoretical  fit  to  the  Mari¬ 
ner  5  spectra  by  instead  assuming 
L0  >  ap  () .  Neither  analyses  accounted 
for  coupling  to  the  inhomogeneous 
background,  however.  As  apparent 
from  the  previous  analysis,  accounting 
for  this  effect  is  instrumental  in  the 
present  case  where  0'1  may  reach  a 
hundred  or  more. 
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Figure  4. 1  On 

1'ime  history  of  signal  intensity  for  the 
recorded  part  of  t he  Mariner  5  Venus  ent¬ 
rance  occultation,  as  adapted  from  Woo  et 
al  (1974) 


Figu  re  4. 1  Oh 

lime  history  of  signal  intensity  for  the 
recorded  part  of  the  Mariner  5  V'enus  exit 
occultation,  as  adapted  from  Woo  et  al 
(1914) 


The  upper  abscissa  uives  the  altiluile  <>!  the  ray  periapsis  above  the  mean  Venusian  surface.  Regions  A 
to  l)  were  suspected  for  their  turbulence  content  and  selected  for  further  examination. 


In  Figures  4.1  1  and  4.12  we  have  overlaid  experimental  spectra  from  regions  A  and  C 
of  the  Mariner  5  exit  occultation  with  theoretical  spectra  calculated  from  equa¬ 
tion  (4.37).  (Kxperimcntal  spectra  from  other  atmospheric  regions  do  not  raise  signifi¬ 
cantly  above  the  noise  level.)  The  relevant  input  parameters  for  this  calculation  were 
obtained  from  Fjcldbo  et  al  (1971),  yielding  aj.  ()  "  1.3-103  m,  vz  =  7.5-1 03  ms'1  and 
v..  "  3-10-  ins’1  for  both  regions  A  and  C.  Since  0V  ?  1.  we  have  v*  s  vv  and 
v*  =  0V,,,  with  0  estimated  from  Figure  4.10b  as  s0.05  for  region  A  and' =0.01  for 
region  C.  Weak  scintillation  theory  is  here  appropriate  since  the  computed  scintillation 
index  for  both  events  is  small  compared  to  unity. 

It  is  apparent  Irom  Figure  4.1  1  that  the  position  of  the  "knee"  in  the  theoretical  and 
experimental  power  spectra  agree  quite  well,  a  fact  strongly  favoring  turbulence  as  the 
cause  ot  the  observed  fluctuations  in  signal  strength.  For  region  A  there  is  a  clear 
discrepancy  between  the  spectral  shapes,  however.  While  the  theoretical  spectrum 
appears  to  match  the  average  high  frequency  slope  of  the  Mariner  5  spectrum,  the 
strong  fringes  in  the  experimental  spectrum  are  not  reproduced  in  the  theoretical 
prediction.  F.vcn  though  both  spectra  display  a  decrease  in  scintillation  power  below 
the  Fresnel  frequency,  they  are  clearly  discrepant  also  in  this  region.  The  reason  for 
this  difference  is  presently  not  clear;  a  contributing  factor  may  be  a  possible  low 
quality  of  the  experimental  spectrum  itself,  however.  Indeed,  the  remarkable  uni- 
formity  of  this  spectrum,  as  compared  to  the  other  Mariner  5  spectrum  in  Figure  4.12 
and  when  contrasted  with  the  90%  confidence  limits  indicated  in  the  figure,  clearly 
cast  doubts  on  the  validity  of  the  spectrum.  Il  may  be  that  the  "zero-fill"  operation 
applied  to  the  original  data  in  order  to  evaluate  the  Fourier  transforms  at  closer¬ 
spaced  intervals  than  would  otherwise  be-  possible  (Woo  etui  1974,  pp  1 700),  has 
introduced  substantial  correlation  between  neighboring  power  estimates  with  a  corres¬ 
ponding  reduction  in  power  variance.  For  a  sharp  spectral  knee  this  correlation  would 
leave  essentially  unchanged  the  position  of  the  knee,  but  would  lend  to  smooth  an\ 
structures,  eg  fringes,  in  the  spectrum.  This  does  not,  however,  explain  the  regular 
fringes  that  still  persist  in  the  experimental  spectrum,  nor  its  strong  decay  at  sub 
Fresnel  frequencies. 
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Figure  4.1 1  Log-amplitude  power  spect¬ 
rum  for  region  A  of  the  Mari¬ 
ner  5  exit  occupation  ( dotted 
curoe),  as  adapted  from  H’oo 
ct  al  (1974) 

Sui>crimpostd  is  the  theoretical  spectrum  (solid 
curse)  based  on  the  actual  occupation  parameters, 
and  displaced  along  the  ordinate  axis  to  vield 
maximum  fit.  The  error  bar  represents  a  VO"!,  i  on 
I'idence  interval. 


Figure  4.12  Log-amplitude  poieer  spect¬ 
rum  for  region  ('.  of  the  Mari¬ 
ner  5  exit  occupation  (dotted 
curve),  as  adapted  from  Woo 
et  al  ( 1974) 

Superimposed  is  the  theoretical  spectrum  Isolid 
curve)  based  on  the  actual  occupation  parameters, 
and  displaced  along  the  ordinate  axis  to  yield 
maximum  fit.  I'hc  error  bar  represents  a  00%  con 
fidence  interval. 


The  experimental  spectrum  shown  in 
Figure  4. 1  2  displays  scatter  among  indi¬ 
vidual  spectral  estimates  that  appears 
more  consistent  with  the  indicated  90% 
confidence  limits.  There  is  an  indica¬ 
tion  that  both  the  predicted  position  of 
the  spectral  knee  and  a  decrease  in 
power  at  lower  frequencies  is  repro¬ 
duced  in  this  event. 

Power  spectra  calculated  by  Woo 
(1975a)  from  the  Mariner  If)  Venus 
occupation  lack  the  low-frequency  part 
below  co(.,  and  also  appear  to  be 
smoothed  at  higher  frequencies  to  an 
extent  that  would  make  a  comparison 
with  the  theoretical  power  spectra  vir¬ 
tually  impossible.  Due  to  the  much 
smaller  turbulent  medium-to-probe  dis¬ 
tances  involved,  the  Pioneer  Venus 
entry  probes  failed  to  detect  any  inten¬ 
sity  scintillations,  establishing  essen¬ 
tially  an  upper  bound  on  the  refracti- 
vitv  fluctuations  in  this  case  (Wo o  et  al 
1979). 


Figure  4. 13  shows  a  power  spectrum 
based  on  the  Mariner  10  Venus  occupa¬ 
tion  profile  as  calculated  by  I.ipa  and 
Tyler  (1979).  The  spectrum  corres¬ 
ponds  in  altitude  approximately  to  the 
region  labeled  A  by  Woo  et  al.  The 
figure  shows  both  the  experimental  and 
theoretical  power  spectra,  together 
with  part  of  the  exit  occupation  inten¬ 
sity  time  profile  and  the  small  segment 
subjected  to  spectral  analysis.  The 
experimental  spectrum  is  calculated  by 
first  fitting  a  polynomial  of  varying 
degree  to  the  intensity  background,  and 
then  subtracting  this  trend  function 
from  the  signal  to  obtain  the  fluctuat¬ 
ing  part.  The  residual  variance  of  the 
signal  relative  to  the  polynomial  trend 
is  shown  as  a  function  of  the  degree  N 
of  the  polynomial  in  Figure  4. 1 4.  Ini¬ 
tially  the  residual  variance  is  seen  to 
decrease  rapidly  with  X,  but  then  be¬ 
coming  virtually  insensitive  to  this  para¬ 
meter  for  N  >  5.  Thus  a  fifth  degree 
polynomial  was  subtracted  from  the 
intensity  time  profile,  and  the  differ¬ 
ence  was  normalized  by  the  trend  func¬ 
tion  and  subjected  to  a  power  spectral 
analysis,  the  result  of  which  is  shown  in 
Figure  4.13.  The  corresponding  theore¬ 
tical  spectrum  has  been  shifted  along 
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Figure  4.  14  Experimental  and  theoretical  power 
spectrum  corresponding  to  the 
apfiroximatelx  Id  s  portion  of  the 
Manner  It)  I'enus  exit  oeeultation 
indicated  in  the  lower  left  part  of 
the  f inure 


the  ordinate  axis  to  yield  maxi 
mum  fit-  Ihe  low-frequency  parts 
are  seen  to  agree  fairly  well.  Also 
the  high  frequency  portions  are 
not  in  obvious  disagreement,  the 
large  scatter  of  the  data  points 
taken  into  account.  The  possible 
discrepancy  between  the  position 
of  the  Fresnel  knees  may,  at  least 
in  part,  he  attributed  to  the  limi- 
ted  statistical  significance  of  the 
experimental  spectrum. 

Preliminary  analysis  of  other  parts 
of  the  Mariner  10  oeeultation  pro¬ 
file  suspected  of  containing  struc¬ 
tures  due  to  turbulence  has  not 
yielded  a  consistent  result.  Most 
of  the  fluctuations  seen  in  the 
deeper  parts  of  the  oeeultation 
profile  turn  out  to  be  noise,  con¬ 
sistent  with  the  findings  of  Woo 
and  colleagues  for  the  Mariner  5 
oeeultation.  At  least  in  one  other 
region  docs  the  power  spectrum 
raise  above  the  background  noise, 
but  yet  does  not  match  its  theore¬ 
tical  prediction.  It  is  possible  that 
this  region  does  not  correspond  to 
well-developed  turbulence,  or  that 
the  statistical  significance  is 
simply  too  low  to  provide  a  mean¬ 
ingful  comparison  with  the  statis¬ 
tical  mean  spectrum  of  the  theore¬ 
tical  analysis. 


.  ,  ,  ,  i  ■  ,  Spectra  of  higher  statistical  signifi- 

indicatcd  in  the  lower  left  part  o,  1  .  ,  .  , 

,  1  cam  e,  essential  to  an  experimental 

the  inure  r  ,  ,  , 

assessment  ol  the  theory  deve¬ 
loped  here,  may  in  principle  be 
obtained  from  the  averaging  of 
power  spectra  from  different  occultations,  or  through  substantial  spatial  averaging  in  a 
single  grazing  oeeultation,  where  the  ray  moves  many  Fresnel  scales  at  essentially 
contain  altitude.  Removal  of  the  background  intensity  from  the  scintillations  poses 
yet  another  fundamental  problem.  When  the  scintillation  spectrum  and  the  trend 
function  have  finite  power  in  a  common  frequency  range,  which  may  normally  be  the 
case,  complete  separation  of  fluctuations  and  trend  is  not  possible  even  in  principle 
(Yaglom  It  is  nevertheless  imperative  that  a  common  and  well-defined  proce¬ 

dure  be  used  when  removal  of  the  trend  is  attempted.  The  technique  suggested  here  is 
both  well-defined  and  statistically  sound  ( Draper  and  Smith  llMiO,  Chatterjec  and  Price 
71. 


Figure  4.14  Rest  variance.  (Jp,  of  intensity  fluctuations  around  a  f>olx noniuil  trend,  as  a 
function  o  f  the  decree  X  *»?  the  txdxnomud 

I  he  lime  vric>  an.tlv/<ii  i  urreNp^mU  to  the  Huinvni  >•!  ihr  humeri  in  Venus  imuli.i 
M«*n  in».lii  .itoi  m  1  mure  l  1  • 


4.5  I  lie  weak  scintillation  approximation  in  planetar\  occullations 


Consistent  with  the  contention  made  m  section  4  d  we  shall  assume  heie  that  tin 

theoretical  results  of  the  previous  sections  are  valid  ’.I  .'  ;  1.  :t  tespei  t  tv  e  ot  the 

magnitude  of  the  phase  fluctuations  at  exit  lr«>m  the  a t m.is'  .her,  below  uc  explore 

to  what  extent  this  weak  scintillation  assumption  max  le  t'  <  sattslied  m  past 

radio  and  stellar  occullations.  examinin''  both  tiiemet  t>  ai  a  i.t  exoenmeutal  evidence. 


He  define  formally  a  critical  oe  collation  depth  a*,  sin  I.  that 
(«?*)>  =  1 


I  he  theoretical  results  developed  here  are  thus  valid  .mU  ; , .  i  .ncultatioti  depths  such 
that  0  0*.  In  terms  of  0*.  the  scintillation  index  at  an  crbitrarv  depth  0  is  from 

equations  1 4.211 )  and  (4.26) 


:0i> 


\V(  1 ,0;j>  I  c0“ 


W(l,0*;p)  0 


l  t 


l 


I 


(4.491 


where  \V  takes  either  of  the  fonns  atven  h\  equations  1 4.  J4h  t  and  l4.2(ib)  for  radio 
and  stellar  o<  cultations,  respectivelv.  In  I  able  4.1  we  have  c  alc  ulated  the  c  ritical 
depth  0*  for  two  past  radio  and  stellar  cue  ill  I  at  ions,  lor  all  cases  we  have  assumed 
'' i' |  > ,'er,  :  1 0  ^  and  I.„  -  11/4.  where  II  ulonu  with  the  free  space  f  resnel  /.one  and 
the  stellar  angular  radius  are  uiven  in  I  able-  4.1.  Since  0 “  1  for  all  four  occ  tiltations, 
equation  (4.49)  simplifies  to 


to  ;or 


(4.40) 


when  0  1.  It  is  apparent  from  I  able  4.1  that  the  weak  scintillation  condition  may 

have  been  violated  in  the  deeper  parts  of  these  occullations  if  turbulence  of  the 
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strength  found  at  the  ~15  km  level 
of  the  earth's  troposphere  is  also  pre¬ 
sent  at  the  higher  altitudes  probed 
on  Venus,  Jupiter  and  Neptune. 

Data  from  the  Mariner  5  Venus 
occultation  have  been  analyzed  by 
Gurvich  (1969),  Golitsyn  and  Gur- 
vich  (1971)  and  by  Woo  et  al  (1974) 
for  the  variance  of  the  log-amplitude 
fluctuations.  Accounting  fora  factor 
of  4  between  this  quantity  and  the 


fable  4.1  Approximate  critical  occultation 
depths  0*  and  minimum  occul¬ 
ted  intensity  $mtn  Jor  two  past 
ratlin  and  stellar  uccultations 
/or  tfn*  stellar  ot  (  filiations  bv  Jupiter  and  \ep 
l line  are  from  laller  rt  til  (1974)  and  Freeman  and 
(1970).  respectively . 


scintillation  index,  Golitsyn  and  Cur- 
vich  (1971)  estimated  the  scintilla¬ 
tion  index  to  be  about  2 - 1  O'2  for 


both  the  night  side  immersion  and 
day  side  emersion.  This  value  is  an 
average  over  altitude  corresponding 
to  an  approximate  pressure  range  of 
1  -a  bar.  Woo  et  al  (1974)  obtain  a 
similar  figure  (— 410'2)  for  the  scin¬ 
tillation  index  in  the  -60  km  altitude  region  of  pronounced  turbulence  identified  in 
their  analysis.  In  his  at  il\ sis  of  turbulence  in  the  Mariner  10  Venus  occultation,  Woo 
(1975a)  found  a  value  =0.1  for  the  scintillation  index  in  the  same  —55  — 70  km 
altitude  region  previously  probed  by  Mariner  5.  While  the  radio  occultation  measure¬ 
ments.  including  that  of  Venera  9  (Timofeeva  et  al  1978).  appear  mutually  consistent, 
estimates  of  the  refractive  index  structure  constant  that  would  yield  much  higher 
values  for  the  scintillation  index  have  been  derived  from  radio  link  analyses  of  the 
Venera  1  8  landers  t Kolosov  et  al  1970,  Yakovlev  et  al  1974).  While  this  discrepancy 
lias  not  yet  been  satisl.u lords  resolved,  the  recent  Pioneer  Veins  entry  probe 
measurements  appear  to  support  the  radio  occultation  data  (Woo  '  1979).  lhe 
spacecraft  measurements  thus  appear  to  predict  a  lower  turbulence  -irength  than 
would  be  predicted  from  extrapolation  of  measurements  in  the  earth's  troposphere. 


No  determination  of  a  scintillation  index  for  the  Pioneer  Jupiter  missions,  for  the 
Voyager  Jupiter  and  Saturn  missions,  or  for  any  of  the  several  past  stellar  occulta- 
tions  have  been  reported.  The  large  excursions  in  intensity  commonly  observed  in 
stellar  occupations  strongly  indicate,  however,  that  non-uniform  atmospheric  struc¬ 
tures.  possibly  turbulence,  would  yield  a  form  al  value  for  the  scintillation  index 
comparable  to  or  larger  than  unity.  It  is  deemed  important,  after  properly  separating 
the  scintillations  from  the  main  structure  of  the  event,  to  calculate  the  scintillation 
index  also  in  such  cases.  Since  numerical  simulations  indicate  that  the  scintillation 
index  saturates  to  unity  also  when  an  inhomogeneous  atmospheric  background  is 
accounted  for  (French  et  al  1980)*,  values  of  <$]>  that  were  larger  than  unity,  for 
example,  would  immediately  rule  out  turbulence  as  the  sole  cause  of  these  intensity 
variations.  Additona!  clues  as  to  the  nature  of  the  underlying  atmospheric  structure 
might  then  be  obtained  from  an  examination  of  the  scintillation  power  spectrum,  as 
discussed  more  fuiiv  in  section  7. 


*  Ret  ent  results  In  this  author,  based  on  .i  slight  ijcncrah/uttou  ol  the  station. irv  |ih.ise  tr<hni,|uc 
developed  in  so  tton  li.  shows  that  the  si  intill.itiun  index  indeed  saturates  to  unity  also  in  the 
present  e  ol  an  inhoinotteneou.s  hat  k^rotind. 
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5  PROFILE  ERRORS  DUE  TO  TURBULENCE 

As  previously  stated,  the  inversion  procedure  to  obtain  profiles  of  refractivity,  tempe¬ 
rature  and  pressure  from  occultation  measurements  assumes  the  atmospheric  refracti¬ 
vity  to  have  a  spherical  distribution.  Turbulence  violates  this  condition,  and  since  the 
perturbations  in  signal  intensity  and  Doppler  frequency  due  to  turbulence  arc  non¬ 
linear  functionals  of  the  generic  refractivity  fluctuations,  such  profiles  will  contain 
both  random  (zero  mean)  and  systematic  (finite  mean)  errors.  We  estimate  in  this 
section  the  magnitude  of  these  errors,  and  shall  find  that  the  random  and  systematic 
errors  are  both  small  in  the  limit  of  weak  scintillations,  consistent  with  the  simple 
argument  given  in  section  2.2. 


5.1  Random  errors  in  atmospheric  profiles  from  radio  or  stellar  intensity  measurements 

Let  =  1  +  v *  be  the  refractive  index  calculated  from  the  formal  inversion  given  by 
equation  (2.1b).  Let  a  ~  a+  p,  where  p/a  '  p/R  -•  1.  lo  lowest  order  in  v  this  equa¬ 
tion  becomes 

v* ( ro )  =  4 *  /  M(a  *  pip1  :dp  (5.1) 

TTj)  i\  p 

We  have  here  used  the  approximation  cosh"1  (a’/a)  '  ( 2 p /" R ) 1  ;,  and  the  argument  a*  p 
of  M  =  cvl)/l I  signifies  that  o  is  the  bending  angle  of  a  ray  with  impact  parameter  a*  p. 
Since  or=  ( 2?rR/I I ) 1  2t'()  for  small  bending  angles  in  a  sphericalK  symmetric  and  iso¬ 
thermal  atmosphere,  direct  substitution  in  equation  (5.1)  yields  i’  _-r0,  as  required  in 
this  case. 


When  turbulence  is  present,  the  signal  intensitv  may  l>e  expanded  in  powers  of  tq  as  in 


0  0,  ♦  0,  ♦  . . .  0,(1  1  <$i  *...)•-  (1  ♦  M)  1 


(5.2) 


The  last  equality  assumes  that  the  signal  intensity  is  still  related  to  the  ray  bending 
angle  as  if  no  turbulence  were  present  in  a  shallow  occultation.  Combination  of  the 
two  previous  equations  yields  for  the  lowest  order  error  e*. ,  relative  to  the  "true” 
refractivity  e0  of  the  ambient  atmosphere 


^  -  1  -  (  -  )‘  v  0'u,l,’p'Mp 

TTl)  K  „  0z(a  *  p) 


(5.3) 


Ihe  argument  a+ p  of  0,  and  0Z  signifies  again  that  these  quantities  are  calculated 
with  respect  to  a  ray  of  impact  parameter  alp.  Clcariv,  as  <0j  >  -*  0  the  refractivity 
error  vanishes  and  v_  -*  tq, . 


Since  is  ,i  stochastic  quantity,  only  its  mean  and  variance  (and  possibly  also 

higher  statistical  moments)  are  of  interest.  1  he  mean  value  is  zero,  since  clearly 
<0i  >  0.  I  h c  fractional  variance  of  this  quantity  is,  from  equation  (5.3) 
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(5.4) 


All  atmospheric  variables  where  .in  argument  has  not  been  assigned  are  referenced  to 
the  rav  impact  jaarameter  a'  ,c  In  establishing  this  equation  we  have  used  the  fact 
that  <0,  (a  +  p,  )0,  (a+p, )  >  <0]  (a)  >exp(-(jt(p,  t-p2))p(p2  p,  ),  where  p  is  the  (nor¬ 

malized)  autocorrelation  function  of  the  (normalized)  intensity  scintillations  at  the 
receiver;  fi" 1  II  and  (3t'  II,  as  before.  In  the  following  we  shall  not  distinguish 
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between  11  and  1  It .  Let 

l  =  /p(Ap)il(Ap)  (5.5) 

i) 

be  the  correlation  length  of  the  intensity  scintillations,  and  let  I-' ( p  |  ,p 2 )  be  a  function 
that  changes  slowly  on  a  length  scale  M,  so  that 

If  ( P ,  .l>a)Pl|>2  Pi)dp|d|>2  -  1/  l-(pi ,pi  )dpi  (5.6) 

on  o 

It  may  be  shown  using  the  definition  (5.5)  together  with  equation  (4.7)  that  I  is  typi- 
iall\  on  the  order  of  aj.  ()  in  a  shallow  occultation.  Therefore,  in  the  range  0<(P.. 
P:>  II  where  equation  (5.4)  derives  its  principal  contribution,  the  co-factor  of  p  in 
the  integrand  of  this  equation  satisfies  the  condition  for  using  relation  (5.6).  For  the 
0  1  region  of  principal  interest  we  thus  obtain 


From  equations  (2.61  ami  (2.7)  it  may  be  shown  that  the  corresponding  fractional 
errors  in  temperature  and  pressure  are  approximately 


1  *i  >  <P*t 


1  A, 

1  «>  > 
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where  and  p*|  are  the  lowest-order  errors  in  temperature  and  pressure,  respecti¬ 
vely.  relative  to  the  "true”  ambient  values  T0  and  p0.  Thus  the  fractional  mean 
square  errors  in  derived  relractivity,  temperature  and  pressure  are  all  small  in  the  limit 
ot  weak  scintillations.  For  stellar  occultations  the  same  result  applies,  except  that  the 
correlation  length  1  may  be  different  from  s  ap  ,,  in  this  case.  If  the  projected  stellar 
radius  as  significantly  exceeds  aj.  „ ,  I  =  as,  and  relations  (5.7)  and  (5.8)  are  still  valid 
if  also  as  •*  11. 

It  is  interesting  to  note  that  if  the  Baum  and  Code  equation  (2.5)  were  used  to 
estimate  the  scale  height  and  hence  the  temperature,  the  corresponding  mean  square 
fractional  error  in  temperature  would  be  approximately  <$]>  for  0  <1.  Thus  the 
integral  inversion  technique  is  less  sensitive  to  intensity  scintillations  than  curve-fitting 
to  the  Baum  and  Code  equation  by  a  factor  1/11  <  1. 


5.2  Random  errors  in  atmospheric  profiles  from  Doppler  measurements 

For  a  spherically  stratified  relractivity  field  the  atmospheric  Doppler  frequency  f  is 
related  to  the  small  bending  angle  a  approximately  as  in 


where  \.  is  the  component  ot  the  spacecraft  velocity  in  the  plane  of  refraction  that  is 
perpendicular  to  the  incoming  ray  asymptote.  When  turbulence  is  present,  the 
Doppler  frequency  may  be  expanded  as  f=  f„  t  f,  +...,  where  f0  =  ©(t'o),  fi  =  0(iq) 
etc.  Since  equation  (5.9)  is  still  used  to  infer  the  bending  angle  also  when  turbulence 
is  present,  the  estimated  angle,  or  -  a  ,  may  be  similarly  expanded  in  powers  of  tq 
according  to  oq  a*o  cVi  +  •  •  •  From  equation  (2.  lb),  and  using  a  similar  expansion 
of  cosh  '(a  /a)  as  in  equation  (5.1),  we  find  that  the  estimated  refractivity  v  is 
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where,  clearly,  v„0  -  Q){v0),  =  0(r’,),etc.  Using  the  standard,  small  angle  expres¬ 

sion  for  a,  equation  (5.10)  is  readily  shown  to  yield  "*o  =  *V  I  he  lowest-order  error 
to  this  ’’true”  refractivity  profile  is  thus  represented  by  equation  (5.11).  In  terms  of 
the  covariance  function  (assuming  local  homogeneity) 


B(Pi  4>2 ) 


<df.  U+  P2_)df,(a^p,J>  =  <(df1(a))>»  p(pj  _pi  } 
dp  2  dp  |  da 


the  mean  square  fractional  error  in  refractivitv  is  in  this  ease 
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(5.12) 


Ihe  first-order  Doppler  perturbation  is  related  to  the  corresponding  perturbations  in 
bending  angle  as  in 


-Xf,  =  vv(3,  -  v7 a- , 


where  J,  and  are  the  first-order  perturbations  in  ray  bending  angle  in  the  xy-  and 
\ /  plane,  respectively. 

Proceeding  from  this  stage  to  evtiluatc  <(df,  dar>  and  p  via  the  statistical  properties 
•  :  p,  and  O;  turns  out  to  he  extremely  difficult,  essentially  for  the  reasons  given  at 
the  end  of  section  4.1.  As  an  alternative  approach  to  a  direct  theoretical  evaluation,  one 
might  estimate  the  right  hand  side  of  equation  (5.12)  from  occultation  data.  In  so 
doing,  care  should  be  taken  to  properly  subtract  from  the  total  atmospheric  Doppler 
the  lowest-order  component,  f0,  due  to  the  ambient  atmosphere,  as  explained  in 
section  4.4.2.  From  a  record  of  t|  versus  impact  parameter  a,  df,  /da  along  with  its 
mean  square  and  correlation  (unction  p  may  he  estimated,  yielding  even¬ 
tually,  through  the  operations  indicated  in  equation  (5.12),  an  estimate  of  the  mean 
square  fractional  error  in  refractivity. 


We  have  shown  previously  that  in  the  evaluation  of  the  mean  square  fluctuations  in 
••tiding  angle,  the  exact  kernel  in  the  field  integral  14.1)  cannot  he  approximated  by 
the  quadratic  I  resncl  kernel.  Also  the  complete  spatial  extent  of  the  turbulent 
medium  will  have  to  be  considered  in  such  a  calculation,  as  explained  in  section  4.1. 
Instead  of  attempting  a  rigorous  calculation,  we  shall  make  the  following  heuristic 
arguments,  which  will  hopefully  yield  the  proper  functional  form  and  order  of  magni¬ 
tude  of  the  right  hand  side  of  equation  (5.12).  It  may  be  argued  on  general  grounds 
that  <a^><<(Ji>  when  0X  <?  !•  We  thus  approximate  -f,  -  (vv/X)/?,,  while  also 
assuming  vv=0(vz).  Due  to  the  assumed  exponential  variation  of" the  rms  turbulent 
refractivity  fluctuations,  the  effective  width  of  the  turbulent  medium  parallel  to  the  limb 
is  iRHj)12.  Perpendicular  to  the  limb  the  effective  width  may  be  shown  to  he  at 
most  on  the  order  of  II  when  0  r  1.  Ihe  smallest  eddy  sizes  capable  of  scattering 
elementary  waves  within  an  aperture  of  these  dimensions  is  therefore  tv  "•  XD/t  and 
kz  "  XD/H  in  the  y-  and  /.-directions,  respectively.  The  correlation  length  of  di3,/da 
can  therefore  he  no  smaller  than  the  least  of  these  dimensions,  if,  (tv.  Also,  when 
•he  turbulent  phase  perturbation  of  the  wave  is  one  radian,  the  rms  turbulent  bending 
angle  can  be  at  most  on  the  order  of  j3 ,  -  X/\v.  Omitting  numerical  factors,  we  thus 
lmd  that  an  approximate  upper  hound  on  the  right  hand  side  of  equation  (5.12)  is  such 
that 
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(5.15) 


Vlopnng  numerical  values  appropriate  for  the  Mariner  10  Venus  and  Pioneer  10 
Jupiter  o<  <  ult.itions.  we  find  that  the  right  hand  side  of  relation  (5.14)  is  approxi¬ 
mately  O  K)  ’  and  4  10'.  respectively.  I  bus  the  refractivity  error  is  small  even  when 
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the  turbulent  phase  perturbation  approaches  the  strong  scattering  regime.  We  also 
observe  from  relation  (5.13)  that  the  refractivity  error  decreases  with  increasing 
planet— spacecraft  distance,  all  other  quantities  remaining  constant.  The  refractivity 
error  incurred  from  the  inversion  of  radio  intensity  measurements  yielded  an  increase  in 
this  error  (proportional  to  D70  for  <t>  <  1;  cf  equations  (4.23)  and  (5.7))  with  pro¬ 
pagation  distance,  and  is  probably  much  larger  than  that  calculated  here  for  Doppler 
measurements  at  any  occultation  depth. 


from  equations  (2.6)  and  (2.7),  and  consistent  with  the  crudeness  of  the  analysis 
presented  here,  we  finally  conclude  that  also  the  fractional  temperature  and  pressure 
errors  are  small,  being  approximately 
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where,  as  before,  T*,  and  p*t  are  the  lowest-order  errors  in  temperature  and  pres¬ 
sure  relative  to  the  ambient  values  T0  and  p0. 


5.3  Systematic  profile  errors 

The  possibility  that  atmospheric  profiles  obtained  from  occultation  measurements 
were  inflicted  by  systematic  errors  in  addition  to  the  random  errors  discussed  here, 
was  first  raised  by  Hubbard  and  Jokipii  (1975).  As  a  result  of  additional  contribu¬ 
tions  and  subsequent  critical  responses  by  other  authors,  certain  aspects  of  the  theory 
of  these  higher-order  effects,  e  g,  their  wavelength  dependence,  are  still  under  debate 
(cf  discussion  and  references  quoted  in  section  1.1).  However,  both  parties  apparently 
agree  that  the  systematic  effects  arc  too  small  to  be  of  much  practical  concern,  and  a 
discussion  of  this  topic  here  mas  therefore  seem  unwarranted.  From  a  theoretical 
point  of  view  these  effects  command  substantial  interest,  however,  and  they  indeed 
appear  to  have  passed  virtually  unnoticed  in  previous  analyses  of  electromagnetic 
wave  propagation  in  randomly  inhomogeneous  media*.  We  shall  therefore  provide 
here  a  brief  discussion  of  these  effects,  developing  first  a  geometrical  optics  result  for 
the  systematic  phase  error  and  proceed  by  generalizing  to  wave  optics.  From  this 
result,  the  form  and  approximate  magnitude  of  the  systematic  errors  in  intensity  and 
bending  angle  will  be  established  along  with  the  resulting  errors  in  the  calculated 
refractivity  profile. 

5.3.1  The  second-order  phase  bias 

Following  the  treatment  of  Eshleman  and  Haugstad  (1977)  and  Haugstad  (1978b), 
consider  a  thin  phase-changing  screen  representing  homogeneous  and  isotropic  turbu¬ 
lence  superimposed  on  a  homogeneous  background,  see  Figure  5.1. 

A  plane  wave  impinging  on  the  screen  as  indicated  is  refracted  through  small  angles 
a,  and  P,  in  the  xz-  and  xy-plancs  of  propagation,  respectively.  The  total  phase 
change  of  the  wave  by  propagating  through  the  screen  to  the  receiver  at  x  =  S/2  +  D  s  I) 
is  kS,  where  in  geometrical  optics  the  optical  path  S  is  simply 
D 

S  =  /  n(x,y(x),  z(x))  ds  (5.15) 


*  Tatarskii  (1971;  pp  2S6)  calculates  the  second-order,  systematic  phase  change  due  to  turbulence, 
but  he  uses  the  result  to  merely  establish  the  conditions  for  convergence  of  the  method  of 
smooth  perturbations.  Keller  (1962)  has  calculated  the  systematic,  second-order  phase  shift  for 
geometrical  optics,  but  his  result  differs  from  that  obtained  here.  Keller  (1977)  has  later  corrcc 
ted  for  an  error  in  his  calculation,  which  brings  his  result  into  exact  agreement  with  that  given 
by  equation  19.22)  here. 
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Figure  5.1  Radio  ivave  propagation  through  homogeneous  and  isotropic  turbulence 
superimposed  on  a  homogeneous  background 

The  spacecraft  is  located  on  the  x-axis  a  distance  D  >  £  from  the  exit  face  of  the  screen. 

Here  n  =  1  +  u  =  1  +  p0  +  p,  is  the  refractive  index,  and 

ds  =  [  1  +  (-{>-r+  (^)2|Wdx  (5.16) 

dx  dx 

is  a  differential  along  the  raypath  represented  by  the  space  curve  v(x),  /.(x). 

The  zero-order  raypath  is  clearly  y0(x)  =  z0(x)  =  0,  since  f0  is  homogeneous.  Denot¬ 
ing  by  Vi(x),  Zi  (x)  the  first-order  perturbations  in  the  raypath  induced  by  the 
turbulence,  vve  find  through  second  order  in  small  quantities  by  expanding  v,  and  ds 
relative  to  the  zero-order  raypath 

S  =  /  (1  +  +  t't  f  t't ..  V|  t  cuz,)(l  4  +  Ja,l)dx  (5.17) 

o  - 

All  quantities  are  here  evaluated  on  the  zero-order  raypath,  i  e,  =  iz,  (x,  0,0), 
fl  y  =  t»i(y  (x,0,0)  etc,  and  subscript  y  or  z  on  c,  denotes  differentiation  with  respect 
to  that  argument.  In  deriving  this  result  we  have  also  used  the  fact  that  0,  =dyt/dx 
and  a!  =  dz,  /dx  for  small  bending  angles  a,  and  0t . 


The  integral  in  equation  (5.17)  should  he  evaluated  subject  to  the  boundary  condi¬ 
tions 


V\  =  t'l.y  =  "i,z  =0  at  x  =  0 
Vi  =  Z|  =  0  at  x  =  1) 

Expanding  formally  S  =  S04-S,  4-S2,  where  S,  =  ©(f,  ),  S2 


(5.18) 

O(t'j),  we  then  obtain 
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(5.19b) 

(5.19c) 
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In  the  last  result  we  have  used  the  fact  that  fa]  dx  ^  l)aj  and/(3^  dxs  D(3f ,  where  Or, 
and  (3,  are  now  the  final,  first-order  bending  angles  in  the  xz-  and  xy  plane,  respecti¬ 
vely.  The  ensemble  averaged  optical  path  length  is  now,  since  <f,>  =  0 

<S>  =  S0+ <S2>  =  S0+ «,  +  C2  (5.20) 

where  we  have  defined,  in  view  of  equation  (5.19c) 

«i  =Il)«a,2>  +  </?,-’>) 

(5.21) 

D  i) 

Kj  =  /  <*'i  v  Vi  >dx  +  /  <iq  /.,>  dx 
o  o 

Integrating  by  parts  the  last  expression  and  applying  the  boundary  conditions  (5.1  8) 
yields 

«2  =  !)(<«, J>  +  <0,:>) 

Thus  there  is  a  net  reduction  in  average  optical  path  by  an  amount 

<S>  S0  =  -  D(<«,2  >  +  <j3,2>)  (5.22) 

relative  to  a  non-turbulent  atmosphere  of  rcfractivity  v 0.  An  alternative  statement  of 
this  result  is  that  the  average  phase  velocity  is  increased  in  the  presence  of  turbulence 
(Haugstad  1976,  F.shlcman  and  Haugstad  1977,  Ilaugstad  1978b).  The  appropriate 
generalization  of  this  result  to  an  inhomogeneous  background  is  (Haugstad  1978b) 

<S2>  =  -  Jj  l>(0y<0r’>  *  0z<a,2>)  (5.23) 

where  a,  and  (3,  are  the  bending  angles  produced  by  turbulence  on  a  homogeneous 
background. 

The  above  derivation  of  <S2>  apparently  allows  the  following  interpretation  to  be 
made.  The  total  change  in  phase  path  has  two  contributions,  Ci  and  82.  The  former 
component  represents  the  increase  in  phase  path  associated  with  the  longer  geometri¬ 
cal  path  to  the  receiver.  This  increase  is  more  than  offset,  however,  by  the  (^-compo¬ 
nent,  which  represents  a  twice  as  large  decrease  in  phase  path  associated  with  the  rays 
"success"  in  traversing  the  refractive  medium  in  such  a  way  as  to  encounter,  on  the 
average,  a  lower  than  average  rcfractivity.  We  observe  from  equation  (5.22)  that  the 
bias  in  phase  path  increases  with  propagation  distance  in  vacuum  between  the  turbu¬ 
lent  medium  and  the  receiver. 

A  wave-optical  calculation  shows  that  the  general  weak  scattering  result  for  the  bias 
in  phase  path  is  (Haugstad  1978c) 


<S2>  =  <X,St>  (5.24) 

where  Xi  and  S,  are  the  first-order,  wave  optical  expressions  for  log-amplitude  and 
phase  path.  Employing  equations  (4.7)  and  (4.8)  and  noting  that  2xi  =  .  we  find 

that 

<S2>  =  trkC /.f  «l>(u_l)sin  (P)duj_  (5.25) 

with  I*  given  as  before.  A  comparison  of  this  result  with  the  geometrical  optics 
expression  (5.23)  shows  that  the  two  results  arc  indeed  identical  in  the  limit  k  -*  A 
further  comparison  of  the  two  results  also  reveals  that  diffraction  strongly  suppresses 
contributions  to  <S2>  from  turbulent  eddy  sizes  smaller  than  the  atmospheric 
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Fresnel  zone.  We  also  find  that  for  Kolmogorov  turbulence  the  bias  in  phase  path 
scales  with  radiation  wavelength  as  A'1*',  implying  that  an  initially  non-dispersive 
medium  becomes  slightly  dispersive  by  the  addition  ol  turbulence  (Haugstad  1978c, d, 
Haugstad  and  Eshleman  1979),  A  result  similar  to  equation  (5.25)  has  also  been 
obtained  by  bee  (1976)  from  an  evaluation  of  the  second-order  scattering  produced 
by  simple  diffraction  grating. 

5.3.2  Systematic  errors  in  atmospheric  profiles  from  Doppler  or  intensity  measure- 


If  the  bias  in  optical  path  changes  spatially,  as  a  result  of  a  spatial  variation  in  either 
t'o  or  <^>,  a  finite  bias  in  both  bending  angle  (and  hence  in  Doppler  frequency)  and 

intensity  mav  result,  as  illustrated  in 
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|vo^>| 
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Figure  5.2  Illustration  of  the  average 
effect  of  turbulence  on  phase 
path,  bending  angle  and  inten¬ 
sity 

The  solid  curve  shows  the  position  of  the  phase 
front  behind  a  homogeneous  ambient  atmosphere 
void  of  turbulence.  The  broken  curve  shows  the 
curved  phase  front  when  turbulence  of  increasing 
strength  in  the  z-direction  is  superimposed. 


intensity  may  result,  as  illustrated  in 
Figure  5.2.  The  turbulence  strength  is 
here  assumed  to  increase  with  atmo¬ 
spheric  depth  so  that  the  average  phase 
advance  is  greater  in  the  denser  region 
of  the  atmosphere  (cf  equations  (5.22) 
or  (5.25)).  It  is  seen  from  the  figure 
that  a  linear  spatial  variation  of 
<v\>  gives  rise  to  a  bias  in  ray  bend¬ 
ing  angle,  whereas  a  non-vanishing 
second  derivative  of  this  quantity  yields 
a  net  focusing  or  dcfocusing  of  the  rays 
and  hence  a  bias  in  signal  intensity. 

Mathematically  the  bias  in  bending 
angle  is  related  to  that  in  phase  path  as 


<cv  >  =  al.ST_>  ,<e2>  = 

9  v  ‘  9z 


where  a2  and  )32  are  the  second-order 
components  of  the  ray  bending  angle  in 
the  xz-  and  xy  planes  of  propagation, 
respectively.  The  exact  relation  be¬ 
tween  the  intensity  bias  and  that  in 
phase  path  is  still  a  matter  of  contro¬ 
versy,  even  for  the  geometrical  optics 


case  (Eshleman  and  Haugstad  1978, 
Haugstad  1978a,  Eshleman  and  Haugstad  1979,  Hubbard  1979).  However,  for  shallow 
occultations  several  lines  of  reasoning  indicate  that  the  intensity  bias  must  be  of  the 
form 


<02  > 


pa<a2> 

0/  /  -r  dx 


(5.27) 


in  both  geometrical  optics  and  wave  optics,  since,  for  reasons  of  energy  conservation, 
a  finite  value  of  <02>  is  conditional  on  a  systematic  spatial  variation  of  the  angular 
bias".  As  the  arguments  over  the  precise  form  of  both  the  angular  bias  and  the 
intensity  bias,  along  with  the  question  of  their  wavelength  dependence,  have  become 
extremely  technical  in  nature,  we  shall  make  no  attempt  here  to  provide  a  balanced 
expose  of  these  issues.  Rather,  since  the  resulting  effect  of  these  biases  on  derived 
profiles  is  agreed  by  both  parties  to  be  very  small,  we  shall  instead  use  simple  order 
of  magnitude  arguments  to  merely  establish  the  approximate  magnitude  of  these 


Relation  (5.27)  also  assumes  that  <j32>  =  0.  This  will  always  be  true  if  both  f0  and  <fj>  have 
strictly  spherical  symmetry,  which  will  be  assumed  here. 


f 

* 
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Considering  the  nominal  occultation  level  (<t>z  -  1/2)  and  disregarding  numerical  fac¬ 
tors,  relations  (5.26)  and  (5.27)  yield 


<a2>  ~  <S2>/H 
<02  >  ~  <S2>/H2 


(5.28) 


where  02  s  02/0z  is  the  fractional  bias  in  intensity,  and  H  is  a  typical  vertical  scale  of 
variation  of  either  <v\>x'7  or  v0.  By  virtue  of  the  assumed  spherical  symmetry  of  these 
quantities,  <02>  =  0.  In  equation  (5.26)  the  z-derivative  at  the  receiver  is  related  to  a 
similar  derivative  in  the  atmosphere  according  to  3/9z(receiver)  =  0zd/9z(atmosphere), 
so  that  relations  (5.28),  being  representative  of  (j>z  —  1/2,  actually  overestimates  the 
biases  in  deeper  atmospheric  regions. 


Refractivity  profiles  obtained  from  radio  Doppler  or  intensity  measurements  will  con¬ 
tain  a  bias,  <r.2>,  that  is  related  to  the  biases  (5.28)  and  to  other  second-order 
effects.  This  refractivity  bias  may  be  obtained  from  equation  (2.1b)  by  formal  expan¬ 
sion  of  v  in  powers  of  iq  and  collecting  all  terms  of  second  order,  yielding 


<v*2> 


_  I 

Doppler  2 

=  I 

Intensity  2 


,  .  9  ?  d<a2>  ,,2  , 

J~i rp  p 

,  Di( 2/  =  7«,2m<«,> 

0 


(5.29a) 

(5.30a) 


A  closer  examination  of  equation  (5.29a),  using  the  previously  derived  approximate 
upper  bound  on  <t>^,>/ ^o>  reveals  that  the  refractivity  bias  is  probably  dominated  by 
the  second  term  in  this  case,  ie  by  the  angular  bias.  Bv  contrast,  the  right  hand  side 
of  equation  (5.30a)  is  dominated  by  the  second  term  in  the  integrand,  yielding 
approximately  for  the  two  cases 


Doppler  ~  (  b'  *  7  1>‘  2  dP 


R  o  dp 


^  2  ^  Intensity  ( 


i 


2RD2 


<0,>p!  2  dp 


(5.29b) 

(5.30b) 


Explicit  evaluation  of  these  relations,  using  the  first  of  relations  (5.28)  along  with  an 
explicit  expression  for  <S2>  obtained  from  equation  (5.25),  yields  the  following 
approximate  results  for  the  fractional,  absolute  bias  in  calculated  refractivity 


^  1^*  ^  HfiLp  j  JF,o^P  <'t'i  ^  (5.30) 

^0  Doppler  Raf*0  L0  r'o 

(<!!jitlj>)  ~<0j>  (5.31) 

P0  Intensity 

Although  approximate  only,  these  results  are  significant  in  several  respects.  The 
former  result  shows  that  the  fractional  refractivity  bias  is  always  small  in  the  weak 
scattering  limit.  For  the  Mariner  10  Venus  and  the  Pioner  10  Jupiter  occultations,  for 
example,  we  estimate  typical  values  '10'6  — 10'5  if  <v\>lv\  =  10‘5  is  a  representative 
turbulence  strength.  We  also  observe  that  the  refractivity  bias  decreases  slightly  with 
propagation  distance  (as  D'1*  for  Kolmogorov  turbulence)  and  scales  with  radio  wave¬ 
length  in  the  same  way  as  the  bias  in  phase  path,  i  e,  as  X',<6  for  Kolmogorov 
turbulence.  By  contrast,  the  fractional  systematic  error  in  refractivity  profiles  derived 
from  intensity  measurements  is  controlled  by  the  entirely  different  quantity  <0  2>, 
which  may  approach  unity  in  the  present  theory.  However,  as  long  as  the  scintilla¬ 
tions  are  weak  also  the  fractional  refractivity  bias  will  be  small,  although  being  much 
larger  than  the  corresponding  error  in  profiles  derived  from  Doppler  measurements. 
Thus  the  fundamentally  different  character  and  magnitude  found  for  the  random 
refractivity  error  for  the  two  cases  also  carry  over  to  the  refractivity  bias,  giving  the 
radio  Doppler  technique  a  further  advantage  over  the  less  accurate  intensity  measure¬ 
ments  as  a  tool  for  remote  probing  of  planetary  atmospheres. 
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6  TURBULENCE  IN  DEEP  RADIO  OCCULTA  HONS 
6.1  General  considerations 

Most  previous  spacecraft  occultations  by  planets  have  involved  fairly  shallow  penetra¬ 
tion  of  the  probing  radio  waves  into  their  atmospheres.  Increased  sensitivity  of  the 
radio  links  of  current  (Voyager  1  and  2)  and  future  spacecraft  will  allow  deeper 
atmospheric  regions  to  be  probed  bv  the  radio  occultation  method,  following 
Eshleman  ct  al  ( 1979a),  we  shall  reserve  the  term  "deep  occultation"  to  denote  situa¬ 
tions  where  atmospheric  refraction  deflects  the  ray  Ivy  about  a  tenth  of  a  planetary 
radius  or  more,  as  viewed  in  the  plane  of  the  sky.  Under  such  circumstances  the 
probing  wave  may  experience  significant  focusing  due  to  the  curved  planetary  limb. 
As  explained  in  sections  3.2  anil  3.3,  in  the  central  regions  behind  the  planet  the 
focusing  may  become  comparable  to  and  eventually  substantially  exceed  the  initial 
defocusing  caused  bv  differential  refraction  in  the  plane  of  propagation.  Correspond¬ 
ing  optical  enhancements  have  been  observed  at  earth  during  the  central  occultation 
of  e-Geminorum  by  Mars  (Elliot  ct  al  1977b). 

For  a  general  refractivitv  field  the  focus,  or  more  properly  the  caustic  surface,  is  the 
locus  of  intersection  of  infinitesimally  close  normals  (rays)  to  the  equiref ractivity 
contours.  If  these  contours  are  everywhere  parallel  to  the  planetary  limb,  this  locus  is 
also  the  evolute  of  the  planetary  limb.  I  bis  circumstance  has  given  rise  to  the  term 
"evolute  flash"  (Eshleman  ct  al  1979a).  which  refers  to  the  brief  but  very  strong  signal 
enhancement  associated  with  a  crossing  ot  the  evolute.  For  an  oblate  planet  characte¬ 
rized  bv  two  orthogonal  radii  of  curvature,  the  locus  resembles  a  four-cusp  cylinder, 
extending  along  a  line  through  the  planetary  center  from  a  certain  minimum  distance 
behind  the  planet  to  infinity.  For  a  strictly  spherical  refractivitv  field  the  evolute 
degenerates  into  a  single  line. 

The  basic  theory  of  stub  evolute  flashes  has  been  worked  out  in  two  recent  publica¬ 
tions.  These  papers  have  either  emphasized  the  potential  information  about  gravitatio¬ 
nal  moments,  planetary  rotation  and  zonal  winds  that  can  be  obtained  from  observa¬ 
tion  of  the  evolute  flash  (Eshleman  ct  al  1979a).  or  the  possibility  of  using  the  gravita¬ 
tional  focus  of  the  sun  for  eavesdropping  and  communication  over  interstellar  dis¬ 
tances  (Eshleman  1979).  fhe  theoretical  predictions  of  these  papers  were  based, 
however,  on  the  assumption  of  a  spherical  or  near-spherical  distribution  of  refracti¬ 
vitv-.  Any  departure  from  such  a  degree  of  symmetry,  as  caused  for  example  by  a 
more  complex  global  distribution  of  refractivitv  or  by  smaller-scale  turbulence,  will  in 
principle  affect  the  predictions  from  this  theory.  The  extent  to  which  this  occurs  is 
investigated  here  by  evaluating  the  general  results  developed  in  section  4  for  the 
scintillation  index  and  intensity  power  spectrum  close  to  the  focal  line,  following 
essentially  the  development  of  Haugstad  (1981).  Fhe  general  results  of  this  analysis 
will  then  be  applied  to  the  actual  conditions  of  the  recent  Voyager  1  Jupiter  focal 
evolute  crossing,  and  also  to  examine  the  effect  of  solar  corona  plasma  irregularities 
on  the  performance  of  the  gravitational  lens  of  the  sun. 

For  the  reasons  given  at  the  end  of  section  4.1.  no  attempt  will  be  made  to  evaluate 
the  fluctuations  in  Doppler  frequency  near  the  fonts. 


6.2  Scintillation  index  and  power  spectrum  close  to  the  focus 

fhe  weak  scattering  scintillation  index  associated  with  the  near-limb  ray  is,  from 
equations  (  4. 10)  and  (4.23) 


H7rk  k  *  //'!>( 


"r''"':",){12,1(asux)|:|dux 


(6.1) 
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where  the  upper  form  is  for  a  coherent  (radio)  source,  and  the  lower  form  for  an 
extended  (stellar)  source  of  projected  angular  radius  ag  above  the  planetary  limb.  For 
the  far-limb  ray,  the  quantities  0V  and  0Z  contained  in  the  definition  of  P  remain  of 
the  same  form,  except  that  0V  -+'|0V|,  since  0V  is  negative  (aD  >  R)  for  the  ray  over 
the  far  limb. 


Partly  because  of  the  ray-optical  nature  of  0 y  and  0Z,  and  in  part  because  of  approxi¬ 
mations  made  in  the  derivation  of  equation  (6.1)  itself,  this  result  is  only  valid 
beyond  a  certain  minimum  distance  from  the  focus.  The  former  constraint  coincides 
with  the  distance  w  s  \I3/R  given  in  relation  (3.21)  as  the  width  of  the  focal  maxi¬ 
mum.  I  he  more  constraining  condition  associated  with  the  neglected  terms  in  the 
stationary  phase  expansion  was  established  in  section  4. l  as  0V  =  (R/rx)  <  (RH) 1/2 / 
ay  (>,  where  rx  is  the  distance  from  the  spacecraft  to  the  tocai  line  in  the  plane  of 
the  sky.  Thus  equation  (6. 1 )  should  not  be  used  close  to  or  within  a  distance  r±  = 
ay  (,(R/II)12  from  the  focal  line. 

Considering  in  the  sequel  the  coherent  source  form  of  equation  (6.1 ),  the  form  of  the 
scintillation  index  in  the  focal  region  is  found  to  depend  strongly  on  which  of  the 
two  inequalities 


aK.h  & 


that  are  satisfied.  When  ay  o  ay  ^  <  Lo.  the  appropriate  form  of  equation  (6.1)  is 
given  by  equations  (4.23a,b).  Evaluation  of  these  equations  yields  for  0Z  ^  1 


(6.2) 


where  K,  -  (tr1  *  ( 2tr| :H»  -  Isccjrp/4I2»’- 1 P2  ((p-1  )/2)|  /  |F((p-3)/2)r(p-l)|.  This  result  is 
valid  for  3  <  p  <  6,  and  for  depths  of  oecultation  such  that  0Z/0V  —  Hr^R-r^R  <  1, 
while  simultaneously  ay  (} ( R/rx) 1  2  <  l,0.  Thus  focusing  around  the  curved  limb  causes 
the  scintillation  index  to  increase  as  rj/''’  (assuming  p  =  11/3)  as  the  focus  is 
approached. 


When  the  oecultation  has  proceeded  to  the  level  where  ay  y,  >  L0,  equation  (6.1)  may 
he  evaluated  subject  to  the  stronger  condition  a,.  h  >  I.0,  to  yield  asymptotically 

<01  >  >  K:l-9-V  <iq>  (6.3) 

A* 

where  K:  Hitr1!  (p/2)/  (p-2)I  ( 3/2 )T  ((p-3)/2)].  According  to  this  result  there  is  no 
explicit  ^-dependence  ol  the  scintillation  index  in  this  regime,  the  only  possible 
change  in  <0j  >  being  due  to  the  very  slow  change  in  <rq>  towards  the  focus.  But 
close  to  the  locus,  where  ay  h  **  I,0  and  0Z .3:  Il/R,  the  change  in  altitude  of  the  rav 
periapsis  as  rA  -*  0  is  a  small  fraction  of  a  scale  height  only.  Thus  <tq>,  and  hence 
also  <01  >.  is  very  nearly  constant  in  this  region.  This  situation  is  depicted  schemati- 
(  ally  in  Figure  6.1 . 


1  he  immediate  reason  for  the  plateau  in  Figure  6.1  is  the  assumption  imbedded  in  the 
turbulence  power  spectrum,  equation  (4.1 6),  that  is  itself  constant  for  u  LJ,1 .  A 
constant  power  spectrum  in  this  region  is  merely  a  convenient  way  of  imposing  a 
finite  variance  on  the  fluctuations  in  refractivity,  however,  since  cannot  increase  as 
u"  all  the  way  towards  u  0  if  p>  1  and  <Tj  >  is  finite.  In  actual  practice  no 
well-defined  outer  scale  exists,  and  any  structure  in  ‘1*  in  the  range  u<I-o  will,  of 
necessity , also alfect  <01  >  dose  to  the  focus. 
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Figure  6.1  Illustration  of  the  qualitative  variation  in  scintillation  index  with  fractional 
occultation  depth,  rJR,  for  a  fixed  planet  to  spacecraft  distance 
R  is  the  planetary  radius,  and  rx  is  the  distance  from  the  spacecraft  to  the  focal  line  in 
the  plane  of  the  sky. 


The  point  source  intensity  scintillation  power  spectrum  is,  from  equation  (4.37) 


2 

W(cj)  =  k2oj  Z 

i  =  l 


xsin2{(^-)  [<>ycos2(0i+  S)  +  0zsin2(0j  +  6)|j- 
/ - - - - _dx  (6.4) 

i  1i  +  (^)2]p/2(x2_1),,2 


As  before,  oa  =  2ir  1/2v*/a p  0  >s  the  Fresnel  frequency,  co0  =  v*/L0,  and  5  =  tan'1  (v*/ 
v*.),  where  vy  = 0yvy,  v£  =  0zv_,  are  the  ray  velocity  components  at  the  ray  periapsis 
corresponding  to  the  velocity  components  vy  and  vz  of  the  spacecraft  in  directions 
parallel  and  perpendicular  to  the  limb,  respectively.  For  the  current  conditions  it  is 
possible  to  derive  an  approximate  analytical  form  of  equation  (6.4).  To  this  end  wc 
first  observe  that  normally  v£/v*  <  (vz/vv)-(10Hrj/R2 )  <  1  in  deep  occultations,  so  that 
6  a  0.  characteristic  of  a  grazing  occultation.  Defining  oa'c  =  coc0y'2,  and  noting  from 
the  definition  of  coc  that  co't.  in  deep  occultations,  we  find  that 


W(cj)  «  co  sin2  ( -^-)  •  /  — - -  (6.5) 

«c  1  (1  +  (^)2]p/2(x2  -  l)1  2 

<x)n 

Examination  of  the  integral  shows  that 


\V(co)  a 


’sin2  ( -—-)  for 
wc 

co1_psin2  (-^4) 

CO  2 


CO  <  co0 


for  co  >  co0 


(6.6) 


Note  that,  while  initially  co0/co'c  ~  (aj.  0/L0)  ^  1  per  assumption,  in  very  deep  occul¬ 
tations  this  ratio  is  reversed,  since  there  co0/co^.  "  (ap /l-0 )•( R/r_,_) 1/3  >  1. 


The  relations  (6.6)  imply  that  in  very  deep  occultations  the  scintillation  power 
spectrum  has  the  form  depicted  qualitatively  in  Figure  6.2,  with  the  bulk  scintillation 
power  confined  within  a  band  coj- ~  co  £  co0.  As  the  focus  is  approached,  both  co^. 
and  co0  increase.  The  bandwidth  also  increases,  however,  since  co^.  ®  (R/rx)*  2  while 
co0  ac  R/r±  and  therefore  increases  faster  than  toc  by  the  factor  (R/rx),/2.  The  increase 
in  co0  towards  the  focus  by  the  factor  R/rx  is  caused  alone  by  the  similar  increase  in 
ray  speed  at  the  periapsis.  The  slower  increase  of  co .  is  due,  however,  to  the  compet¬ 
ing  expansion  of  the  horizontal  Fresnel  scale  by  a  factor  (R/rx)12,  which,  if  acting 
alone  would  reduce  cot  by  a  similar  amount,  and  the  increased  ray  speed  which  is 
proportional  to  R/rx. 
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figure  (i.2  Qualitative  sketch  of  the  point  source  intensity  scintillation  power  spect¬ 
rum  in  eery  deep  occultations 

The  characteristic  frequencies  COt  and  to 0  arc  defined  in  the  text. 


In  the  temporal  domain,  the  bandlimitcd  power  spectrum  described  above  corresponds 
to  typical  scintillation  time  scales  less  than,  or  at  most  equal  to,  that  associated  with 
the  low  frequency  boundary  of  the  pass-band,  that  is  .  Thus  cj^.  is  also  a 

characteristic  temporal  frequency  in  very  deep  radio  occultations. 

by  way  of  illustration,  it  is  readily  proved  using  the  previous  relations  that  one  scale 
height  from  the  focus,  where  the  background  intensity  has  risen  to  the  unocculted 
lc  el  (0=1)  according  to  equation  (3.1 1),  the  lower  boundary  of  the  scintillation 
pass  band  is  at  about  100  II/.,  if  vv  =  1  kms'1 .  a j.-  ()  =  1  km  and  R  =  10<I  km  arc  repre¬ 
sentative  values.  Such  rapid  scintillations  may  severely  impair  the  ability  to  track  the 
radio  signal  through  an  evolute  Hash,  in  particular  if  the  scintillations  are  also  strong 

( <0 i >  2  i). 


6.3  Applications 

6.3.1  1  he  Voyager  1  Jupiter  encounter 

Ihe  first  positive  identification  of  the  refractive  focus  of  a  planetary  atmosphere  is 
due  to  Elliot  et  al  ( 1977b),  who  observed  the  "central  flash"  of  e  Geminorum  as  it  was 
occulted  by  Mars  on  8  April  1976.  By  contrast,  the  first  attempted  radio  transmission 
from  the  focus  occurred  by  the  recent  passage  of  Voyager  1  through  the  focal  evolute 
of  Jupiter.  According  to  Martin  et  a  I  (1980),  the  radio  signal  displayed  several  peaks 
in  intensity  near  the  expected  position  of  the  evolute  crossing,  but  none  could  be 
identified  as  the  evolute  flash.  The  absence  of  a  detectable  evolute  flash  led  Martin 
and  colleagues  to  establish  a  lower  bound  on  the  microwave  absorption,  most  likely 
caused  by  Nil,,  near  the  5 4.3  bar  level  of  the  ray  periapsis  when  crossing  of  the 
evolute  occurred.  We  estimate  below  the  weak  scattering  scintillation  index  during 
this  evolute  crossing,  and  argue  that  its  high  value  may  have  contributed  to  the 
unsuccessful  attempt  to  identify  the  flash  in  this  case. 

A  formal  calculation  of  the  scintillation  index  near  focus  requires  an  estimate  of 
'  t'j  >.  We  will  assume  that 

(6.7) 
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i 
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where  n0  is  the  ambient  refract ivity  at  the  rav  periapsis  during  the  evolute  crossing. 
1  he  value  in  relation  (ft. 7)  corresponds  roughlv  to  turludence  o)  the  strength  found 
immediately  above  the  earth’s  boundary  layer  (Button  cl  til  1972).  if  also  the  outer 
scale  and  the  ambient  temperature  at  the  ray  periapsis  are  equated  to  1  km  and 
200  K,  respectively.  The  ambient  refraetivity  is  related  to  its  value  at  standard  tempe¬ 
rature  and  pressure  (S  I  P)  by 

_  =  IT  [o 

*'e).STP  1*0  1 

where  (p0,!'0)  refer  to  the  standard  condition  (SIP),  fsing  <'(,.s  1 1'  0.000  I  36  as  lor 

a  pure  H: -atmosphere,  we  find  that  vu  "  S-IO'4,  which  yields  the  value  '6  -  10'* 
for  0]>.  Putting  t  =  (RH)‘*  and  11  7  30  km,  we  obtain  for  the  S-band  (X  0.1:1  m) 
oecultation  using  equation  (6.31  and  the  previous  estimates 

<0j  >  a  6  •  10'  (6.8 1 

Hus  high  value  of  the  weak  scattering  scintillation  index  suggests  that  the  scintilla¬ 
tions  have  saturated  to  unity  well  before  the  focus  is  reached.  If  real,  such  strong 
scintillations  would  significantly  change  the  maximum  intensity  relative  to  the  ideal 
prediction  for  an  oblate  planet  (c  /  equation  la  of  Kshlcman  c  l  til  1070a).  further¬ 
more,  the  matched  filter  technique  employed  bv  Martin  <•/  ul  llOSO)  to  detect  the 
flash,  requires  signal  coherence  over  a  certain  minimum  time.  In  view  of  the-  discus¬ 
sion  in  section  6.2  it  is  clear,  however,  that  tins  provision  cannot  be  met  close  to  the 
focus  where  the  scintillation  time  scale  approaches  zero.  1  hits  scintillations  ntav  have 
impaired  the  mere  detection  ol  the  evolute  Hash,  above  and  beyond  the  fundamental 
difficulty  imposed  bv  the  very  strong  Ml,  absorption  <>i  the  S-bancl  radio  signal 
along  the  atmospheric  portion  of  the  lavpath. 

6.3.2  The  gravitational  lens  of  the  sun 

Kshlcman  (1979)  gives  the  maximum  lot  used  mteuMiv  lot  this  case  as 


0 


max 


2ff  kr„ 


1 6.9 1 


where  r,<  is  the  gravitational  radius  ol  the  sun  c  t _’uv>  me  I  his  tcsult  is  derived  for  a 
plane,  coherent  wave  < > t  wavenumber  k  delict  ted  thimigh  an  angle'  o  2r,.  a  ( a  is  the 
impact  parameter)  bv  the  solar  gravitational  held.  \nv  sol.u  oblateness  is  too  small  to 
significantly  degrade  the  maximum  mtensific  aiion  given  bv  equation  Stull 

degradation  max  be  caused,  however.  1  > v  the  sol.u  corona  plasma,  which  imposes 
additional  bending  on  the  wave  uteri v  from  the  lotus  bv  an  amount  depending  on  the 
local  plasma  gradient  and  the  radio  wavelength 


Inference  from  the  continuum  of  scale  si/es  louud  :n  the  solar  vv  mil  plasma  (eg, 
Jokipii  19711),  suggests  that  also  the  solar  corona  is  inhomogeneous  on  a  continuum 
of  scale  sizes  ranging  from  the  Dcbvc  length  to  those  associated  with  the  global 
dimensions  of  the  corona. 


I  he  approximate  effec  t  of  large  st  ale  coronal  plasma  stun  lures  on  ntav  be 

appreciated  bv  finding  the  condition  under  which  the  additional  phase  shift  imposed 
by  the  plasma  is  less  than  X/2.  I  his  is  the  condition  t hot  the  I  resnel  zone  In¬ 
coherent  and  the  signal  intensity  being  esscnltullv  unaffected  bv  the  plasma  for 
larger  phase  shifts  the  Kresnc-I  zone  breaks  up  into  two  or  more  non  e  ounc  e  ted  pails, 
with  a  corresponding  drop  in  maximum  intensilv  below  the  value  given  1  > v  c-cpta 
lion  lb-9l.  We  thus  require 

k:-  >  •»:  I  It,  10) 

where  <Sj  >  is  the  mean  septare  c  hange  in  phase  path  of  the  radio  wave  imposed  In 
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the  coronal  plasma.  An  additional  constraint  to  ensure  the  validity  of  equation  (6-9) 
is  that  <0]>  <1.  In  very  deep  occultations,  where  aj.  h  S>  L0,  this  condition  is,  from 
equation  (6.3 1 

<0f>  '  k,  <tq>  <1  (6.1 1) 

X' 

It  is  important  to  observe  that  as  Ion”  as  a j.  h  >  L0,  the  two  previous  conditions  are 
essentially  identical  regardless  of  the  magnitude  of  L0.  Indeed,  from  an  explicit  evalu¬ 
ation  of  <S^>  it  follows  that  in  this  case 

<01  >  --  2  k:<S}>  (6.12) 

It  is  thus  immaterial  which  of  conditions  (6. 1  0)  or  (6.11)  that  arc  imposed  to  ensure 
the  validity  of  equation  (6.9).  A  result  similar  to  equation  (6. 1 2)  has  been  obtained 
by  Salpeter  (1967)  for  a  one-dimensional  turbulence  model.  Equation  (6.12) 
extends  this  result  to  three-dimensional  turbulence  on  an  inhomogeneous  back¬ 
ground  when  the  leading  dimension  of  the  atmospheric  Fresnel  zone  substantially 
exceeds  the  outer  scale  of  turbulence. 

Assuming  I.0  '  R  ( R  is  the  solar  radius),  so  that  the  rms  phase  change  is  comparable 
to  that  imposed  by  the  ui'cra ge  corona,  we  calculate  below  the  coronal  change  in 
phase  path  corresponding  to  the  modified  Baumbach-AUen  plasma  model  proposed  by 
lyler  rt  nl  ( 1977),  which  accounts  for  latitudinal  variations  and  the  effect  of  the  solar 
wind  plasma.  Denoting  by  p  and  0  the  heliocentric  distance  (in  units  of  the  solar 
radius)  and  solar  latitude,  respectively,  this  model  asserts  that 

\  I  (  +  1-5’,)1014  +-  ~  —  10"  |  (cos20  +  ,4Tsin20)1'2  (6.13) 

e  p16  ph  p-  od 


Figure  (i.j  I  (dues  of  the  phase  shift  «  induced  by  the  solar  corona  plasma,  for  various 
fractional  heliocentric  distances  p  and  wavelengths  X 
I  hr  corresponding  fractional  distance  D  I)/Dmjn  to  the  focus  is  also  indicated. 


r>6 


where  N'c  is  the  density  of  electrons  per  in3 .  The  associated  refractivity  is 

v  = -4.48- 10-'6  X2  Nc  (6.14) 

with  X  in  units  of  meters.  Using  equation  (6.14)  anil  the  previous  equations  we  have 
calculated  in  Figure  6.3  the  ratio  K  =  < S 2  > 1/2  / X  for  different  values  of  p  and  X.  The 
very  large  values  attained  by  «  for  even  large  heliocentric  distances  (p  =  10)  and  small 
wavelengths  (X  a  1  mm),  strongly  suggest  that  coronal  plasma  irregularities  both 
reduce  the  coherent  size  of  the  Fresnel  zone  (and  thus  0,nax)  and  cause  severe  signal 
fading  as  the  focus  is  approached. 

According  to  Figure  6.3,  near-optimum  conditions  can  be  restored  either  by  moving 
the  ray  periapsis  to  larger  heliocentric  distances,  or  by  reducing  the  wavelength  below 
about  10‘4  m.  The  former  alternative  is  not  particularly  attractive,  since  the  distance 
to  the  focus  behind  the  sun  increases  as  p1 ,  and  is  already  10()Dmjn  at  p  =  10 
(l)min~  350  AU  is  the  distance  to  the  nearest  focus  corresponding  to  p=  1).  The 
apparent  solution  is  therefore  to  use  very  short  wavelengths,  provided  present  techno¬ 
logical  difficulties  of  sub-millimeter  communication  can  be  overcome. 


7  DIAGNOSTICS  OF  SMALL-SCALE  ATMOSPHERIC  STRUCTURES 
7.1  General  considerations 

turbulence  in  the  tree  atmosphere  owes  its  existence  to,  (i)  thermal  instabilities,  (ii) 
wind  shear,  or  (iii)  the  breaking  of  large-scale,  vertically  propagating  waves,  such  as 
internal  gravity  waves.  1'hese  are  physically  distinct  processes  that  mav  lead  to  differ¬ 
ent  spatial  distributions  and  characteristics  of  the  turbulence.  For  example,  a  convec¬ 
tive  instability  may  be  more  likely  to  occur  near  the  equatorial  plane  on  the  dav  side 
where  the  solar  heating  rate  is  a  maximum.  On  the  other  hand,  wind  shear  and 
the  associated  turbulence  may  be  sharply  confined  in  altitude,  but  may  not  be  other¬ 
wise  restricted  to  particular  planetographic  locations.  A  common  feature  of  well- 
developed  turbulence,  regardless  of  the  characteristics  of  the  source,  is  the  universal 
form  of  the  turbulence  power  spectrum  in  the  inertial  subrange  where  LJ,1  <  u  <  5J,1. 

1  he  magnitude  of  the  outer  scale,  being  a  characteristic  dimension  of  the  input 
region,  depends  directly  on  the  specifics  of  the  source,  however. 

In  the  turbulence  model  (4.16),  the  two  quantities  <rq>  and  I.0  are  of  principal 
interest.  One  might  also  want  to  estimate  the  slope  parameter  p  from  the  occultation 
data,  as  has  been  done  by  Woo  ct  al  (1976,  1977).  It  should  be  emphasized,  however, 
that  a  statistically  significant  departure  of  p  from  the  universal  value  p=  11/3  in  the 
inertial  subrange,  would  imply  the  existence  of  turbulence  that  is  not  well-developed. 
1  hat  is,  it  is  of  an  intermittent  and  nonstationary  nature.  The  positive  identification 
of  such  intermittent  sources  of  turbulence  would  in  itself  be  important,  but  the 
particular  spectral  shape  recorded  would  be  directly  linked  to  the  specific  event  and 
would  not  bear  any  universal  significance.  The  quantities  <vf>  and  L„  in  well- 
developed  turbulence  would  be  of  much  greater  interest,  since  the  atmospheric  pro¬ 
cesses  generating  the  turbulence  will  be  of  a  more  stationary  character  in  this  case. 
I  he  specifics  of  such  processes,  such  as  the  turbulence  strength  and  characteristic 
dimensions  of  the  source,  may  provide  important  observational  constraints  on  theore¬ 
tical  circulation  models  of  planetary  atmospheres,  as  suggested  initially  bv  Woo  ct  ill 
(1974). 

In  order  to  estimate  <v*>  and  L0  from  occultation  data,  one  should  in  principle 
consider  quantities  that  exhibit  the  strongest  possible  dependence  on  the  magnitude  of 
these  parameters.  For  the  former,  the  scintillation  index  and  the  associated  power  spec- 
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trum  is  such  a  quantity,  while  for  the  outer  stale  the  phase  stimulations  is  the 
natural  quantity  to  examine.  Also  other  signal  characteristics  may  be  used,  hut  will 
not  he  considered  here.* 

We  discuss  helow  specifically  how  the  characteristics  of  small-scale  atmospheric  struc¬ 
tures  can  he  obtained  from  occultation  measurements.  It  is  deemed  important  to 
explore  both  the  limitations  and  capabilities  of  the  occultation  technique  in  this 
regard,  since  once  this  potential  is  realized,  such  experiments  may  he  planned  to  take 
1  idlest  possible  advantage  of  their  intrinsic  capabilities.  Parts  of  the  discussion  here 
closely  parallel  that  of  llaugstad  (1979a). 


7.2  Estimation  of  the  turbulence  strength 

Before  any  "vviggly"  occultation  curve  is  analv/.ed  lor  its  expected  turbulence  con¬ 
tent.  it  should  be  ascertained  that  the  recorded  small-scale  structures  are  in  fact  due 
to  turbulence  and  not  to  other  atmospheric  structures  (e  g.  "layers"),  or  to  numerical 
artifacts  introduced  during  the  data  recording  or  the  subsequent  analysis.  The  best 
way  to  make  such  a  check,  would  probably  be  by  examining  the  intensity  power 
spectrum  and  the  scintillation  index.  Assuming  that  the  intensity  background  has 
been  properly  removed,  the  weak  scintillation  power  spectrum  should  conform,  with¬ 
in  proper  limits  ol  statistical  significance,  with  one  of  the  characteristic  shapes  dis¬ 
played  in  f  igures  4.5  to  4.9,  for  various  geometries  and  depths  of  occultation.  As  an 
alternative,  one  might  identify  the  Fresnel  "knee"  at  u>c  in  the  intensity  power  spectrum, 
and  determine  if  it  scales  with  wavelength  as  A1  2  for  radio  occultations.  Within  wide 
bounds  on  the  spectral  slope  parameter  p  (U  <  p<  (i).  this  scaling  is  independent  of 
p.  In  stellar  occultation  power  spectra,  the  position  and  shape  of  the  corresponding 
"knee"  depend  on  the  relative  si/e  of  the  projected  stellar  radius  and  the  frec-space 
Fresnel  /.one.  11  as  a,  ().  the  "knee"  is  determined  by  as  and  is  therefore  essentially 
independent  ot  wavelength.  For  as  a(,  ().  stellar  scintillation  power  spectra  would  in 
all  respects  resemble  closelv  those  of  a  coherent  source. 

A  final  check  on  the  expected  authenticity  of  the  turbulence  is  provided  bv  the 
scintillation  index,  which  scales  with  wavelength  as  A'h  for  values  small  compared  to 
unit  v . 

Once  well-developed  turbulence  has  been  reliably  identified,  its  strength  and  spatial 
variation  mav  in  principle  be  obtained  from  the  observations  via  the  theoretical  rela¬ 
tionships  developed  in  section  4.4.  Finis,  for  shallow  radio  occultations  and  <  1 . 
equation  I  l.L’aal  yields 


<rf>  j  1.,;  ’  11,  1 2  K  1 2  l)'**'  X76  <<p,2 >  (7.1) 

where  £  is  a  numerical  constant.  It  an  estimate  of  l.„  has  been  obtained,  as  discussed 
in  the  next  section,  equation  ( 7. 1 )  provides  an  estimate  of  <iq2>  in  terms  of  the 
scintillation  index  it  also  llt  is  known.  A  natural  first  choice  of  Ut,  consistent  with 
the  atmospheric  model  adopted  here,  would  be  11,^11.  A  more  realistic  estimate  can 
be  obtained  from  the  spatial  variation  of  <$f>  itself  if  this  quantity  has  been  com¬ 
puted  over  a  few  turbulence  scale  heights  or  more.  Indeed,  if  good  quality  scintilla¬ 
tion  data  are  available  over  a  range  in  altitude,  the  actual  variation  of  the  scintillation 
index  may  be  used  to  improve  the  exponential  turbulence  model  suggested  here,  and 
thereby  yield  a  more  accurate  estimate  also  of  <»;>■ 


Wottrf  al  ( 1977)  have  suggested  using  the  mutual  i  nherence  function  and  also  the  quantity  <X|S,^ 
(Woo  I97:»l»l.  which  is  proportional  to  <S27>.  Ihe  latter  quantity  does  not  appear  to  offer  any 
advantage*  over  the  s<  m filiation  index.  In  the  weak  scattering  regime,  incorporation  of  the  inhomo¬ 
geneous  bac  kground,  not  considered  bv  Woo  rt  al  (1977),  shows  that  the  mutual  coherence 
func  tion  is  essentially  identical  to  the  phase  scintillations  in  terms  of  its  potential  usefulness  for 
determining  or  1,0. 
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For  stellar  occultations  the  analogous  relationship  to  equation  ( 7. 1 )  may  he  obtained 
from  equation  (4.25).  As  previously  noted,  however,  the  explicit  form  of  this  relation 
ship  depends  critically  on  the  ratio  as/.i|  ().  No  general  relationship  ol  the  lorm  (7.1) 
ran  therefore  he  given  in  this  case,  except  when  the  ratio  is  either  very  small  or  very 
large  compared  to  unity.  In  the  former  case  equation  (7. 1 )  clearly  applies,  while  in 
the  latter  case  the  appropriate  relationship  is  obtained  from  equation  (4.2(>a)  as 


<*T> 


T}  I -o  ' 


R1 2  D'2  <cy  > 


(7.2) 


where  r;  is  a  numerical  constant.  Again  knowledge  of  1,„  and  II,  yields  <c2  >  in  tenns 
of  the  scintillation  index  for  the  high  altitude  regions  probed  in  stellar  occultation 
experiments. 

Instead  of  <»>2>  one  might  want  to  estimate  the  refractive  index  structure  constant, 
C,2.  Since 


Dn(Ar)  =•  C,2|Ari2-'  ,  V„  «  Ar-i  l.„ 


where  I)n  is  the  refractive  index  structure  function,  it  follows  from  the  definitions  of 
I)„  and  I,0  that 
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From  the  previous  results  it  therefore  follows  that 
'kr1I;'  2R•|'2l)';ih^7,' 


,Ksas'-’lli,»R-'-l>-^ 


(7.4) 


The  upper  form  with  kr  —  2f  is  for  radio  occultations,  and  the  lower  form  with 
xs^'2ri  for  stellar  occultations.  The  unknown  value  of  l.„  is  here  absorbed  into  the 
structure  constant,  flic  temperature  structure  constant  corresponding  to  C2  for  as¬ 
sumed  adiabatic  conditions  is  given  bv  equation  (4.21).  For  radio  occultations  we  find 
that 

/x  , 

Ci  "  2f(7  1 ) -  1 .7 1  Ij '  ’  R  1  '  I )  ^ X ' (7.3) 

1  r- 

i) 

and  with  a  corresponding  result  for  stellar  occultations  when  the  conditions  leading  to 
equation  (7.2)  are  satisfied. 

Thus  equations  I  7. 1 )  and  (7.2)  vield  the  variance  of  the  refractive  index  fluctuations 
in  terms  of  the  scintillation  index  once  l.„  and  1 1 1  have  been  stipulated.  From  the 
relationships  (7.4)  and  (7.5),  rough  estimates  of  the  appropriate  structure  constants 
may  also  be  obtained  from  a  measurement  ol  the  scintillation  index,  using  prior 
knowledge  of  II,  only. 

From  the  analysis  of  section  b  it  may  be  shown  that  formal  relationships  closely 
similar  to  those  developed  here  also  hold  for  deep  occultations. 


59 


7.3  Estimation  of  the  outer  scale 


As  we  have  seen,  an  approximate  measure  of  L0  is  a  prerequisite  for  estimating  <«'2> 
from  the  scintillation  data.  However,  by  virtue  of  its  character  as  a  typical  dimension 
of  the  macroscopic  flow  generating  the  turbulence,  knowledge  of  the  magnitude  of 
I.0  is  also  important  in  its  own  right. 


1'he  natural  quantity  to  examine  in  this  case  is  the  phase  scintillation  power  spec¬ 
trum,  as  given  in  analytic  form  anti  computed  for  different  occultation  geometries 
and  depths  in  section  4.4.  The  relevant  feature  of  the  spectra  is  the  "knee”  at 
oj  s  -  v^/L.o.  Since  v*  is  known,  the  position  of  the  ’’knee”  is  directly  related  to 
the  magnitude  of  1,0.  It  is  moreover  significant  that  the  position  of  the  "knee”  is 
independent  of  both  occultation  depth  and  geometry  for  shallow  occupations.  In 
actual  practice  the  turbulence  power  spectrum  will  not  have  the  idealized  form  as¬ 
sumed  here,  with  constant  power  for  u  I/J.  A  significant  departure ,  one  way  or  the 
other,  from  the  average  high  frequency  slope  of  the  phase  spectrum  is  probably  the 
feature  that  in  practice  signifies  the  wavenumber  location  of  the  input  region. 


In  deep  occupations  the  estimation  of  I.0  is  less  easy.  1  he  situation  is  simplest  for 
moderately  deep  occupations,  where  ay  q  <  ayj,  <  I.0.  In  this  case  the  spectral  shape 
corresponds  closely  to  that  of  Figure  4.4,  with  the  transition  to  small  fringes  at  cj  =  u>c 
approaching  oo0  as  ay  h  increases  towards  L„.  At  still  greater  depths  of  penetration  of 
the  ray,  where  ay  h>  I.0,  the  phase  power  spectrum  assumes  the  asymptotic  form 
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The  power  spectrum  thus  has  the  form  displayed  qualitatively  in  Figure  7.1.  Identifi¬ 
cation  of  I.,,  is  seen  to  depend  on  the  possibility  of  detecting  the  onset  of  the  co’8  3 
decrease  in  scintillation  power  above  >co0.  The  large  fringes  in  the  power  spectrum  in 
this  region  may  hamper  this  identification,  a  problem  likely  to  be  compounded  by  the 
expected  rapid  imrease  of  oj„  itself  near  the  focus.  On  the  other  hand,  if  the  transi¬ 
tion  in  spectral  shape  at  u)  "  has  been  identified,  determination  of  1,0  is  readily 
accomplished  by  the  usual  relationship  ton  =  v*/ 1,0. 


Figure  7.1  Qualitative  sketch  of  the  phase  scintillation  power  spectrum  in  very,'  deep 
occultatio  ns 

I  he  i  h.ir.n  tenstu  frc,jucn<  ics  u)(  and  OJ0  arc  defined  in  the  text. 
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7.4  Planetary  rotation  and  atmospheric  winds 

We  discuss  here  the  possibility  of  extracting  from  the  scintillation  measurements 
information  about  planetary  rotation  and  atmospheric  winds.  In  the  earth’s  atmo¬ 
sphere  this  technique  has  been  used  extensively,  having  in  fact  reached  the  level  of 
sophistication  where  a  horizontal  wind  profile  has  been  obtained  (Lee  1974).  Scintil¬ 
lation  measurements  have  also  been  used  to  derive  bulk  flow  velocities  of  the  solar 
corona  and  solar  wind  (<  g.  Woo  1978). 

1  lie  theoretical  feasibility  of  such  measurements  rests  in  all  cases  on  the  identification 
of  characteristic  features  in  the  scintillation  power  spectrum  whose  location  or  shape 
depends  on  the  atmospheric  velocity  across  the  line  of  sight.  Considering  the  fixed- 
axes  coordinate  system  of  Figure  2.1,  the  ray  velocity  components  at  periapsis  relative 
to  the  atmospheric  motion  are 


\*  -  Vv 


rt  vs 
vz  z 


(7.7) 


where  a  superscript  s  has  now  been  added  to  the  spacecraft  velocity  components  to 
distinguish  from  the  corresponding  components  of  atmospheric  motion.  v‘‘.  and  vi. 
The  latter  may  be  dominated  by  atmospheric  rotation  or  by  winds,  depending  on 
planetary  rotation  rate  and  occultation  geometry. 

While  in  the  previous  sections  we  have  assumed  both  va  and  v‘‘  to  be  negligible,  we 
here  ask  for  effects  that  can  be  attributed  to  the  finite  values  of  these  quantities.  For 
a  spherically  symmetric  atmosphere  the  rotational  part  of  v^'  is  zero.  Since  vertical 
wind  velocities  are  normally  very  small  compared  to  their  horizontal  counterparts,  we 
shall  assume  in  the  sequel  that  va  1  0.  In  this  case 


f-n<Vy  vv'):  >  0-v;.:r 
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(7.8b) 


where  the  relative  ray  speed  vK  and  obliquity  factor  m  determine  the  location  and 
shape,  respectively  ,  of  the  intensity  power  spectrum,  as  depicted  in  Figures  4.6  4.9 

lor  shallow  occultations.  These  spectra  are  also  valid  in  deep  occultations,  provided 
aK  !■<>•  as  explained  in  section  6.2.  In  very  deep  occultations,  where  a^  ^  9>  L0, 
the  form  of  the  intensity  power  spectrum  must  in  general  be  established  from  the 
general  result  (4.37).  If  v'  is  such  that  itnl  1.  the  power  spectrum  has  the  charac¬ 
teristic  shape  display  ed  in  Figure  6.2,  with  the  low  frequency  boundary  related  to  the 
characteristic  frequency  u/  -  27T1  2v*/aj.  h,  and  hence  to  v*  and  va. 

It  is  clear  from  e(]uations  ( 7.8a)  and  (7.8b)  that  in  order  to  relate  v  a  to  the  location 
of  the  spectral  "knee"  at  co  =  ojc  (u>{  for  deep  occultations),  or  to  the  spectral  shape, 
it  is  necessary  that  iva|  be  at  least  comparable  to  the  component  0vv‘,  due  to  the 
spacecraft  alone.  This'  fact  accentuates  in  the  present  context  the' importance  of 
near-central  (/ e,  deep)  occultations.  Indeed,  if  over  some  range  of  the  occultation 
0V |vyl  ^  Ivyi,  then  both  v*  and  m  depend  crucially  on  v a  whenever  0zlv* .1  <  :va|, 
which  is  likely  to  occur  at  some  stage  of  the  occultation'.  In  this  regime,  both  the 
shape  and  location  of  the  power  spectrum  yield  indirect  information  on  the  magni¬ 
tude  of  v-y.  In  some  cases  the  magnitude  of  v“  may  be  such  as  to  dominate  v*  and  m 
even  for  shallow  occultations.  For  Jupiter,  for  example,  the  rotational  component 
alone  is  about  13  kins'*  for  a  meridional  occultation  over  the  equator.  For  an  occul¬ 
tation  over  the  polar  regions,  the  rotational  component  is  very  small  and  va  may  here 
be  dominated  by  the  wind  component. 


It  is  apparent  from  the  sketehy  discussion  above  that  information  on  wind  speeds  and 
atmospheric  rotation  rates  are  indeed  within  the  capability  of  the  occultation  tech¬ 
nique.  However,  since  such  measurements  depend  critically  on  details  of  the  occulta¬ 
tion  geometry,  it  is  important  to  fully  examine  these  capabilities  in  the  planning  stage 
of  any  spacecraft  mission. 

7.5  Practical  considerations 

In  the  discussion  of  the  previous  sections  we  have  only  briefly  alluded  to  the  several 
practical  problems  involved  in  the  estimation  of  turbulence  parameters  and  atmo¬ 
spheric  motion  from  occultation  measurements. 

The  major  practical  obstacle  in  the  estimation  of  <fj>  appears  to  be  the  limited 
statistical  significance  of  the  scintillation  index.  This  problem  is  most  acute  for  a 
central  or  near-central  occultation,  since  a  meaningful  evaluation  of  <$}>  must 
employ  scintillation  data  covering  a  range  in  altitude  of  less  than  a  turbulence  scale 
height.  In  a  near-central  occultation.  the  number  of  independent  ’’ripples”  in  the 
intensity-time  record  is  at  most  \0  ^  Ht /a k.o .  since  ap  ()  is  an  approximate  lower 
bound  on  the  correlation  length  of  the  intensity  fluctuations  in  this  case.  Since  N0  is 
typically  on  the  order  of  10.  it  is  clear  that  the  statistical  significance  of  <0j>,  and 
hence  of  <1-4  >.  may  indeed  be  very  poor  in  this  case.  The  problem  of  low  statistical 
significance  is  compounded  by  the  principal  difficulty  of  subtracting  the  intensity 
background  from  the  fluctuation  part  of  the  signal,  as  discussed  more  fully  in  sec¬ 
tion  4.4.2.  Because  also  1  l/.i j.  0  =  CO(IO),  the  intensity  trend  function  may  have  signifi¬ 
cant  power  at  spatial  frequencies  — a .  implying  that  irremovable  contamination  of 
the  scintillation  index  and  power  spectrum  from  the  ambient  background  may  be  the 
normal  situation  in  near-central  occultations.  Despite  this  fundamental  difficulty,  it  is 
important  that  an  efficient  and  common  procedure  be  used  to  separate  the  intensity 
background  from  the  fluctuation  part  in  order  to  minimize  this  problem,  as  pre¬ 
viously  argued. 

1  he  situation  is  entirely  different  for  shallow  occultations.  In  this  case  the  ray  has 
only  a  very  slow  vertical  motion,  but  covers  a  large  horizontal  distance  at  essentially 
constant  altitude.  Indeed,  in  the  central  portion  of  the  occultation  the  ray  moves  a 
horizontal  distance  -A’  S  (KH,)1*  before  the  altitude  of  the  ray  periapsis  is  changed 
by  one  turbulence  scale  height.  In  this  case  N ,  Q/af.o  *%>  N„,  so  that  considerable 
spatial  averaging  and  corresponding  high  statistical  significance  of  <<?;>  and  </-q>  is 
possible. 

I  he  practical  estimation  of  l.„  is  probably  most  hampered  by  the  inherent  lack  of  a 
well-defined  value  of  this  quantity.  Indeed,  the  transition  from  spatial  scales  in  the 
inertial  subrange  to  the  input-range  will  for  all  practical  cases  be  gradual,  and  any 
quantity  or  explicit  result  that  depend  on  the  precise  value  of  L0  should  therefore  be 
handled  with  corresponding  caution.  Apart  from  this  inherent  imprecision  of  L0 
itself,  the  statistical  argument  presented  above  apply  for  similar  reasons  also  to  the 
estimation  of  this  quantity. 

I  he  estimation  of  atmospheric  bulk  flow  and  turbulence  parameters  in  deep  occulta¬ 
tions  is  affected  by  statistical  limitations  in  the  same  way  as  described  above  for 
shallow  occultations.  Extremely  dose  to  the  focus,  where  both  co0  and  coj.  change 
rapidly  due  to  the  rapid  change  in  tf>v,  the  effective  length  of  the  scintillation  time 
scries,  over  which  cq  may  be  considered  stationary,  depends  critically  on  the  various 
parameters  involved  and  must  be  assessed  for  each  particular  case. 

We  have  described  above  some  important  constraints  of  statistical  nature  that  are 
likely  to  affect  the  practical  estimation  ol  turbulence  parameters  and  atmospheric 
motions  from  occultation  measurements.  Neither  the  potential  nor  the  statistical  limi¬ 
tations  of  this  technique  appear  to  have  been  fully  appreciated  in  the  past.  It  is 


iU'cmcil  important  that  both  aspects  Ik*  addressed  ai  an  earliest  possible  sluice  in  ilu* 
pl.mninu  dI  |d.mct.trv  IK  1  >\  missions,  so  .is  in  realize  In  t-.it  h  fast-  the  maximum 
ji  ilt  nu.il  ol  ilu-  occult. itiem  lechnicpir. 


S  DISCISSION  AND  I  l  RUH  R  IM.RSIM  O  l  i\  I  S 

ll  is  i  n>  >  I  ,ii  >  It-  I  .it  t  ill. it  wliilt-  t  •  1 1  ft  Is  ol  turlmleiu  f  on  taebo  and  optical  wave  ptopa- 
;.iljon  m  ilif  t'.iiili  s  ai mosjilitTf  and  ionosphere  ha\t-  attracted  i  nnsidrraldt-  ailt-nlion 
o\ri  i  'it-  ji.tst  .It  t  .nit's.  onl\  v tux  rcccntlv  lias  (Ins  issue  Keen  pursued  svstcm.it  if.tllv  in 
at  i  nliui  ion  nn  isuivnienls  ol  pl.mfi.iiv  .itmi>s|ilnTes.  I  Ins  l.u  k  of  interest  is  naturalK 
explained  t •  > i  i.ulio  m  t  ullations  I > \  the  short  hision  ol  this  kind  ol  measurement s. 
ami  also  l>\  the  l.u  t  that  the  lirst  analvscs  ol  such  ohscrv  at  ions  (/ e.  the  Mariner 
missions  to  Mars  ami  Venusl  produced  atmospheric  politics  that  were  laruclv  accept- 
i *  It  .  I  he  [•  »\s  inleiest  is  less  cxpln  aide  Ku  stellat  t  u  tult.it  ions,  the  milch  longer 
lust.ux  ot  w  h i<  i i  is  marretl  In  "spilt \"  Imninosiiv  ptoliles  indicative  ol  Uirhulenee  or 
othei  small  si  ale  atmospheric  strut  tines. 

I  he  pit-sent  an.iKsis  attempts  to  catcuM'i/c  am!  an.iK  /e  the  various  aspects  ol  the 
overall  problem  |>oscd  In  turbulence  in  plain  tan  .unios|)heies  subjt-t  t  t-t  I  to  radio  or 
stellat  otiiillatioii  measurements.  I  Inis  the  prim.irv  task  has  been  to  extend  existin'; 
theorv  id  elec  1 1  oma  ;net  it  wave  prop. i -.at  ion  m  a  t. union)  medium  to  acoiiinl  lot  the 
inhomogeneous  ambient  at  nn  >sj  diei  e  upon  which  the  iiiihuleme  is  siiperiurposed. 
VV  hilt-  i  n-  ;i  i  a.i  b  It  in  terrestrial  prop.i.aiiou  c  xpet  munis  nnoKiu;  shot  l  r,i\  paths,  in 
i  ad  to  and  si  ellar  m  til  It  at  ion  ex  pci  line  nt  s  oni  mu.  h  longer  jiat  I  is  sin  h  eoupltn-;  to  the 
ambient  atmosphere  tn  .omes  veiv  impott.nti.  I  lie  theoretic  al  results  derived  here  lot 

the  si  nit  ilia!  ion  index  a  tie.  with  similar  results  obtained  smiultaneouslv  lllubbard 

c  1 1  *  7  S  ii  oi  late:  ilti.o  .  •  I  USUI  !«\  nine  ant  Inns.  In  the  lorin  expressed  here 

the  t  o  * .  1 1  cite.;  I.|  tin  c  . .  u  p  I :  n  on  '  he  s.  nil  dial  ion  index  is  i  cadih  assessed.  He  hat  e 

I*  mild  tliat  III  shallow  ...  iihat  |o-,s  ,  niiplui  :  to  the  nnbient  I  >.n  k  ;t  mi  ml  reduces  the 
sc  mi  ill. i ;  ions  :  >  \  at  in. .si  a  In  t.n  .1  two  i  point  soilin',  while  lor  an  extended 

sled  u  soma  e  the-  fee  1 1  n  t  ion  varies  bet  vv  e.-n  this  limit  when  tile-  projec  ted  stellar  t  ad  ills 
.  .imputed  t  the  ta.lins  ol  the  lice  space  I  resiiel  /one.  to  a  maximum  teduc 
o;;  e.  appi ..  xiniat  elv  a  la.  lot  ..!  v.iivvlnn  this  i.itio  is  l.u-;e.  I  liese  lairlv  moderate 
c  han  -  tie  .  onti.cstc'd  I )  v  vet'  .east  ini  i.i  I  c  hail  ns  in  the  shape  ol  the  inleusitv  power 
spec  t : . ; : :  l .  which  is  now  i  sti.-n  turn  lion  ol  both  >>c  c  ult.it  ion  de|>th  and  ;cmnctrv. 

I  liese  .  om  hisioti'.  as  initliei  dc  i  uled  in  the  appropriate  so  lions,  express  a  major 

t  lie.  >t.  ..  ll  : '  'Cl’.t  ..  I  the  :u  esc  III  w  "I  I- 


i  id,  sis  also  show  ,  ’he  ; !  n  a  net  all  I  s,  in  t  illal  mn  I  umulus  del  iv  ed  het  e  should 
'Ken  a  n  i!. .  .s  :n  pi.  s'  nn  l.u  ■  1  Is’  v  t.n  approui  I)  vv  lienev  c  r  dll  I  e  lent  ial  rcti  ac  I  ion 
s  'he  iveia  .  i :  1 1 .  •  1 1  s  1 1  x  in  a  siniilii.ml  utioiltil  \n.dv  sis  ol  ladto  and  optical 
1 1 1 .  nis  . !  c  i  :  r  i  ■_  .  >c  c  u'.i  at  i.  >:i  obsei  v  at  tons  should  then  tote  alwavs  proceed  Irom 
'  ict  all  '  c  . :  I." mulls  deiiv.,!  here  but  also  otliei.  teitestii.il  applic  at  mns  mav  be 
lot  tin's.-  i .  suits,  x.  1 1  edit  e  to  sal.  'hie  <  ■•iiiniunic  ation.  where  the  connnuilica 
III:  is  III  I  .a  .  I  ■;  III  .  i  lev  the  e. Hill's  1 1  111.  >s;  illt'l  e  ol  Ionosphere.  111. IV  involve  Slllll 
. !d  e:it  id  i  el  i  ’ :  ■ .  •:  l’i .  ., .  a  it  ton  ..  I  lie  e  art  h'  at  mospliere.  vv  here  the  t  av 

is  i.jous  .  >  1  thdinii  mveisioii  mav  also  v  leld  substantial  values  ol  \|  nD'll. 
h  hi  ;h  tin  el!.,  live  distam.  D  between  tin  medium  and  the  transmitter  or 

■  is  sin  ill  in  this  case  li  is  possihl,  .  Im.illv.  that  a  version  <>l  the  present 
a!i,.;i  lice..:',  mav  I'll. I  application  in  acoustic  lot  elastic  I  wave  |>ropa  ;at  ion  ill 
r:  -  .:  mo -jihete.  in  tin  ocean  oi  in  the  earth's  interim,  whenever  siillicieni 

.  touts  due  to  mom  ilu  s  in  c  iihei  tc-mpei  nine  or  densitv  aie  present. 

ditilu  it  ic  >n  1 1  id  i  valuation  ol  the  sec  oml  older,  svsteinatic  prop.iu.it  ion  ellects 

■  list it 1 1 1 c-  a  piuicip.il  icsiilt.  thou;h  ol  ni.unlv  theoretical  interest,  ('.outran  to 
u  id  o  '  ill  ol  lllibb.ud  and  J  ■>  k  1 1 '  i  i  ll‘)7"i),  who  concluded  that  the  aver.i"e 


phase  was  retarded  hv  turhulem e  relative  l < >  a  mm  turbulent  innlmm  willi  similar 
average  ret  rue  t  iv  it  v,  we  have  found  (rum  a  smijile  \;eomel  rieal  optics  .inulv  sis  llial  till’ 
average  phase  is  m  lari  (ulruiu  t'<l  in  i  lie  presence  ol  /.ei'o-iiiean  turbulence.  I.xtcnsion 
ot  this  result  to  wave optics  has  Inrtlierniore  demonstrated  a  weak  wavelength  depen 
deuce  id  the  Inas  m  phase  path  Im  phase  speed),  \  leldiii1;  the  important  result  that  an 
initialK  mm  dispersive  inedmin  heroines  slightly  dispersive  hv  the  addition  o|  turhu 
Icme,  When  a  svsteinatic  spatial  variation  in  turbulence  stiem;th  hi  in  the  ambient 
rclructivitv  is  present,  a  bias  in  rav  bcnclmv;  .male  and  intensity  mav  also  result. 

While  these  sec  oml  order  propagation  elite  ts  were  lotind  to  v  ield  only  msianit  it  ant 
errors  in  occullation  experiments,  in  other,  more  controlled  experiments  rc<|uirina 
extreme  at  t mac  v  ,  suc  h  eilects  mav  have  to  he  ac  counted  lor.  A  potential  c  andidate 
mav  he  relativistic-  timedelav  experiments,  m  which  the  i.njiath  is  intercepted  hv  a 
turbulent  neutral  or  loni/eel  medium. 

It  is  a  principal  conclusion  limn  the  present  work  that  the  propagation  effects  dis¬ 
cussed  ahove  product-  onlv  veiv  maramal  eriors  in  proliles  ol  temperature  and  pres- 
sure  derived  limn  m  t  iiil.il ion  measurements,  ptovided  that  the  weak  sc  intillation 
assumption  lundamental  to  these  calculations  is  not  violated.  In  view  ol  the  lar-ze 
•  ■riot’s  m  the  Ihonei  I11  and  I  I  |  •  1 1  >  n  c  -  r  proliles.  as  caused  In  the  tna-tnil  it  at  ion  effect 
discussed  in  section  this  conclusion  is  in  no  wav  obvious.  I  he  result  is  however, 

lullv  consistent  with  i.idio  .  >c  c  ult  at  imi  expel  iiueiil',  when-  individual  I  -p  proliles 
limn  the  same  atm<"p!ieie  ! •  * . 1 1 1  exhibit  subsi.uiti.il  mutual  agreement,  and  are  also 
hioadlv  consistent  with  other  observations  where  comparisons  can  In  made-  Stellar 
oi  c  ult  at  ion  pi  o  I  lies  are  t  \  pit  icd  hv  a  m tic  It  I  -  idler  iuteiisit  v  flue  I  uat  ion  lev  el.  probably 
violatin',  the  weak  scintillation  assumption,  hut  are  also  more  susceptible  to  instru¬ 
mental  etrots.  It  is  clear  that  errors  ol  this  tvpe  plaved  a  major  role  in  the  several 
independent  recordings  o|  the  ,cSc  mpn  oec  ult  a  ( ions  hv  |upitei,  which  displaved  sub 
si.inn.il  dis. i -n  a  1 1 1  c  ■  1 1 1  i  see.  Ibinten  and  N'evc-rka  1  *  <  7  *  • :  1  Inis  stellar  oc  c  tiltaliou  pro 
tiles,  thoiidi  clearlv  containin':  si/.,hle  errors,  do  not  dec  ess.u  ilv  contradict  the'  con¬ 
clusion  that  the  part  ot  these-  errors  that  can  be  attributed  to  turbulence  is  small 
w  henev ei  the  weak  sc  intillation  assumption  ts  satisfied. 

I  he  natural  and  important  extension  •> I  the  themv  developed  in  this  work  will  he  to 
tin  or:  oi  ate  stum;  scatterm;  elicits.  I  Ins  extension  should  also  include  an  evaluation, 
hv  t  lieoi  et  n  a  I  or  nuinc  rn  a!  anal  v  sis.  ot  the  mean  s.  uiaie  flue  t  uat  ions  in  bend  i  tv;  ar.  ;le. 
essential  to  a  more  ti  aiimis  assessment  ol  luihuleme  elicits  ill  radio  Doppler  lliea- 
sii  i  cm  c  - 1 1 1 '  dium,  icotli  weak  and  stroll:  sc  if  tit  in-:  conditions.  1  hc-se  tasks,  prohahlv 
in-in  :  ol  a  inline:  a  ai  nat  me.  indeed  constitute  tin  i in  tin  ipa!  remain m  ;  problems  in  the 
thc-dt',  o|  -i.iihulc-nii  eilects  m  tadio  and  'lell.u  occultation  expet  itiients.  I  he  rapid 
•n.  tc-.isc-  ■  ■:  o..  . .  I : i •  i i  1  1 1 1  •  - 1  i!..iex  well  iiinos'.i'iern  de  pth  in  i  idi<>  m  (  ult.iiions.  a|on: 

will  'ill  .  OIIUII- .  i*  ist-l  V  c-i  I  “spike  sten  IIHC-"  o|  slejl.il  •••  i  llll.ltlotl  lldlt  i  llIVC-s. 
Ill  ike  ll  lie  S II  .  c :  >  1c  to  milslli  tins.  oils  VI  o|ol|s|v  Ol  the  two  I  asks,  a  iuuuetlc.il 

i  ah  tila t  'a  >n  ■  >!  ’ '  - e  n a  i  , ua i  e  t  In.  t  na 1 1>  c ns  , u  :  > •  ■  1 1 <  i in  ,  a  u  . a  o  j  n .  >i  m  .  >S v  t  in  simp .  ■- 1 

I  ■  e  •  ■ 1  -o  ■ :  i  1 1 .  i :  s,  i : : et in  nu  o  'a  I .  evaluation  o!  a  double  in!  -  .1  a  I  is  i  e.  p  i :  i  •  a  i . 

while  !o!  i  in  '■  !  w  o  .  I  line : '  si.  u  i.c  1  'iiode!  a  Imu  tuple  I  n  t  >  a.h  o  involved.  1 V  x 

i  a-iti.ia.  tile  a  mi  ll. a!  Imi  index  involves  a  hoi!  tuple  liltc  aal  even  III  the  ope 
■  imn  usi'-!iai  c  i .  *  - .  whil'  -i  two  dimetl-sfc  ui'  ev  llll.ltlotl  o|  an  ei  hi  t :  i ;  >  1 1  ■  mte  ■  i"  a  1  Is 
c  dad  !  oi  \  - 1 1 !  1 1  *  i  c  t  eipio.nii  mvolvm;  less  computation,  will  m  to  calculate  a 

s ,  ■  i  ’  i !  !  e  d ,  i  ’ :  •  • :  c  ii:-o|viii  .  e  I  -  nit  ’  uple  1 1 1 1  ■  ■  ■  t  a  I  ill  a  two  dmn  usi.  mi  I  model  i  ol  tin 

mtensitv  i  oi  h  s:  i.  .mini  in  a  spatial  u  ali/ation  o|  tin  tiirhulcin  e.  and  then  lorm  the 

:m  is  s.jii'i.  i !  r  h,  in'  i  -isii  v  I  i  .'in  m  ensem  Me-  e  it  sii.  Ii  .  .de  111  ct  i.-lls 


I  mui  ;h  e  ve  n  lies  i.iliei  appio.i.  Ii  mav  lead  In  piohihitivc  e  omputatlmi  times,  a  cllli'i  1 
in.ilvric.il  i  ali  iilal  imi  mav  In  no  c.isiei.  Indeed,  even  ordmarv  siion  ;  scintillation 
t  hco  '  I  v  his  '  ii  1 1  v  Te'ee  ill  I  \  heen  deveioped  t«'  a  staoe  where  cert. mi  Imiiliut:  expressions 
lot  the  sc  inlill.il  Ion  index  an  av.iilalde  (see.  i  1  ante  l'17o,  1  re  and  Jokipli 

1'*.  >a.lii.  Jud  un •;  from  the  sulisi.uin.il  added  complexwv  hv  m 1 1 odiic  in  ;  the  mlionto 
a  nemis  .imlnenl  .ilmosphen  ut  l  he  vveak  w  .itterni .;  case,  a  similar  approach  lit  the 
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strong  scintillation  regime  may  be  beyond  fcasibilit'  .  Ii  is  possible,  however,  that  a 
more  heuristic  approach,  eg  of  the  sort  suggested  In  Salpetcr  (1967)  for  turbulence 
on  a  homogeneous  background,  can  be  used  to  treat  strong  scintillations  with  the 
ambient  background  included.  Together  with  selected  numerical  computations,  first  in 
a  one-dimensional  model,  a  combined  approach  of  lids  «ui  may  be  the  only  reali¬ 
zable  path  towards  a  better  understanding  of  the  various  aspect >  of  strong  scattering 
effects  in  occultation  experiments*. 

An  extension  of  the  present  theory  to  strong  scattering  should  also  encompass  the 
systematic  effects.  In  strong  scattering  they  are  likely  to  involve  all  powers  of  ,  and 
should  therefore  exhibit  entirely  different  characteristics  in  this  regime.  Remarkably, 
the  intensity  bias  may  be  the  easiest  quantity  to  generalize  to  strong  scattering.  This 
contention  rests  of  the  result  of  a  preliminary  analysis  (Haugstad  1979b)  which  indi¬ 
cates  that  the  mutual  coherence  function,  which  reduces  to  the  average  intensity 
when  the  two  field  points  of  its  argument  coalesce,  may  be  extended  in  analytical 
form  to  strong  scattering  and  an  inhomogeneous  background.  No  similar  indications 
have  been  found  for  the  phase  and  bending  angle.  These  latter  quantities  are  more 
easily  obtained  from  numerical  calculations,  since  they  both  involve  only  double 
integrals  for  even  a  two-dimensional  propagation  model.  While  a  numerical  evaluation 
may  not  be  satisfactory  for  assessing  the  qualitative  aspects  of  the  systematic  effects 
in  strong  scattering,  this  approach  is  fully  acceptable  for  assessing  their  approximate 
quantitative  impact  on  occultation  experiments. 

The  fundamental  problem  and  question  motivating  the  present  analysis  has  been  to 
what  extent  turbulence  in  planetary  atmospheres  limit  the  accuracy  of  profiles  of 
refractivity,  temperature  and  pressure  obtained  from  occultation  measurements.  We 
have  contended,  either  explicitly  from  rigorous  analysis,  or  from  approximate  or 
heuristic  reasoning,  that  in  the  limit  of  weak  scattering  this  influence  is  very  small. 
An  answer  to  the  basic  question  should  be  pursued  also  for  strong  scattering  condi¬ 
tions,  and  we  have  indicated  how  this  might  be  done.  The  fact  that  turbulence 
effects  arc  in  all  probability  small,  at  least  in  weak  scattering,  implies  an  added 
potential  of  the  occultation  technique  to  also  yield  information  on  small-scale  atmo¬ 
spheric  structures.  Thus,  even  though  the  effects  of  turbulence  on  derived  profiles  arc 
small,  its  various  imprints  on  the  occultation  data  may  be  extracted  and  subjected  to 
a  separate  analysis.  Subject  to  the  practical  limitations  discussed  in  section  7.5,  this 
analysis  would  yield  estimates  of  the  turbulence  strength  and  the  outer  scale  in  both 
shallow  and  deep  occultations.  In  very  deep  occultations,  where  atmospheric  rotation 
or  velocities  intrinsic  to  the  atmosphere  dominate  over  the  ray  speed,  information 
also  on  the  two  former  quantities  are  within  the  capabilities  of  the  occultation  tech¬ 
nique.  Since  the  formal,  weak  scattering  scintillation  index  is  typically  very  high  in 
such  experiments,  as  examplified  by  the  estimate  for  the  Voyager  1  Jupiter  evolute 
crossing,  an  extension  to  strong  scattering  conditions  of  the  theory  developed  here  is 
a  prerequisite  for  measurements  of  this  kind.  If  this  has  been  successfully  done,  and  if 
other  sources  of  error  can  be  controlled  and  kept  on  a  low  level,  careful  planning  and 
interpretation  of  occultation  observations  may  provide  a  very  important  source  of 
information  also  on  atmospheric  circulation  and  dynamics. 


*  Recent  results  by  this  author,  based  on  a  slight  generalization  of  the  stationary  phase  technique 
developed  in  section  :l.  docs  show,  however,  that  some  analytical,  strong  scattering  results  for  a 
homogeneous  background  (c  j,  the  mutual  coherence  function  and  the  asymptotic  scintillation 
index  anil  power  spectrum)  may  be  simply  modified  to  accomodate  the  effect  of  an  inhomogeneous 
background. 
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